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1.  SYNOPSIS  OF  CONTRACT  ACTIVITIES 


During  the  past  five  contract  years,  we  have  completed  all  the  milestones  outlined  in  our 
original  program.  These  include  the  study  of  (i)  the  generation  of  ion  conics  of  ionospheric 
origin,  (ii)  nonlinear  wave  propagations  resulting  from  auroral  electron  beam-plasma 
interactions,  (iii)  turbulent  heating  of  charged  particles  in  the  ionosphere  and  magnetosphere  due 
to  wave-particle  interactions,  and  (iv)  the  evolution  of  counterstreaming  electrons  in  the  auroral 
zone. 


During  this  contract  period,  we  published  a  total  of  30  sc,>.  ntific  papers  addressing  the 
various  states  of  plasma  turbulence  and  associated  stochastic  heating  processes  of  charged 
panicles  in  the  geoplasma  environment.  Members  of  the  MIT  group  engaged  in  the  research  of 
ionospheric-magnetospheric  processes  have  been  invited  by  numerous  intemational/national 
conferences  and  scientific  institutions  to  deliver  invited  and  review  lectures  on  the  scientific 
findings  contained  in  these  papers. 

The  close  proximity  between  MIT  and  the  Geophysics  Directorate,  Phillips  Laboratory, 
USAF  has  allowed  direct  personal  exchanges  and  collaborations  in  several  interesting  research 
interactions.  These  included  visits  by  and  discussions  with  Drs.  J.R.  Jasperse,  J.M.  Retterer,  M. 
Heinemann,  Su.  Basu  and  others.  In  addition,  our  research  efforts  have  been  enhanced 
considerably  by  periodic  visits  of  a  number  of  established  scholars  across  the  nation  and  from 
abroad.  Our  contributions  in  the  study  of  the  turbulence  and  headng  processes  in  the  ionosphere 
and  magnetosphere  were  particularly  influenced  by  Professor  Paul  Kintner  of  the  Cornell 
University,  Drs.  William  Peterson  and  Dave  Klumpar  of  the  Lockheed  Palo  Alto  Research 
Laboratory,  Dr.  Mats  Andre  of  the  Swedish  Space  Institute,  Dr.  Hannu  Koskinen  of  the  Finish 
Meteorlogical  Institute,  Dr.  Dave  Winningham  of  the  Southwest  Research  Institute,  Dr.  Roger 
Amoldy  of  the  University  of  New  Hampshire,  Dr.  Noah  Hershkowitz  of  the  University  of 


Wisconsin,  Dr.  C.T.  Dum  of  the  Max-Planck  Institute  in  Extraterrestrial  Physics,  and  Dr.  Don 
Gumett  of  the  University  of  Iowa. 

Five  scientific  symposia  on  the  "Physics  of  Space  Plasmas",  four  Cambridge  Workshops 
on  various  topics  of  geoplasma  physics,  and  a  Chapman  Conference  on  Ion  Acceleration  in  the 
Ionosphere  and  Magnetosphere  were  organized  during  the  period.  These  activities  resulted  in  the 
publication  of  six  books  based  on  the  scientific  presentations  at  these  conferences. 

The  content  of  this  report  is  organized  as  follows.  In  Sections  2-5,  we  discuss  in  detail  the 
major  contributions  of  our  research  activities.  In  Section  6,  a  complete  listing  of  the  scientific 
publications  of  the  MIT  group  during  the  past  five  years  is  provided.  A  listing  of  the  titles  oe  thr 
books  of  the  published  proceedings  is  given  in  Section  7. 
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2.  PARTICLE  ACCELERATION  BY 

ELECTROMAGNETIC  ION  CYCLOTRON  TURBULENCE 


G.  B.  Crew  and  Tom  Chang 
MIT  Center  for  Space  Research,  Cambridge,  MA  02139 

ABSTRACT 

Low  frequency  electromagnetic  turbulence  is  proving  to  be  an  impor¬ 
tant  source  of  energy  for  the  acceleration  of  ions  in  various  regions  of  the 
Earth's  magnetosphere.  In  particular  it  has  been  shown  to  account  for  some 
of  the  energetic  oxygen  conics  found  in  the  auroral  regions,  and  a  convincing 
case  is  being  built  for  its  role  in  the  cusp/cleft  region  of  the  magnetosphere. 
The  transfer  of  energy  from  the  waves  to  the  particles  is  efficiently  accom¬ 
plished  through  ion  cyclotron  resonance  with  the  left-hand  polarized  com¬ 
ponent  of  the  turbulence,  and  the  result  of  the  interaction  is  a  heating  of  the 
particle  distribution.  In  this  tutorial  review,  we  shall  present  a  general 
theoretical  treatment  of  ion  cyclotron  resonance  heating  in  a  weakly  inhomo¬ 
geneous  magnetic  geometry  and  then  proceed  to  examine  the  formation  of 
auroral  ion  conics  in  somewhat  greater  detail.  For  the  auroral  case,  the  pro¬ 
perties  of  the  electric  field  spectral  density  and  the  Earth’s  dipolar  magnetic 
field  allow  the  introduction  of  a  similarity  transformation  which  results  in  a 
considerable  simplification  of  the  analysis  for  the  altitude-asymptotic  form  of 
the  conic  distribution.  The  merit  of  this  approach  is  that  it  makes  it  possible 
to  directly  compare  the  theory  with  observations,  and  the  agreement  is  found 
to  be  excellent. 


I.  INTRODUCTION 

This  is  a  tutorial  review  of  the  theoretical  progress  that  has  been  made 
towards  an  understanding  of  the  energetic  ion  populations  found  within  the 
Earth’s  ionosphere-magnetosphere  system  which  are  generically  referred  to  as 
ion  conics.  In  fact,  there  is  such  a  wide  variety  of  ion  conic  observations  as 
well  as  viable  theoretical  mechanisms  to  explain  them  that  it  is  not  possible 
to  do  them  justice  within  the  bounds  of  an  hour’s  talk  or  an  article  of  this 
length.  Thus  we  shall  here  restrict  our  attention  to  a  class  of  these 
phenomena  where  the  energization  in  accomplished  through  an  ion  cyclotron 
resonance  interaction  with  electromagnetic  plasma  turbulence. 

We  are  concerned,  then,  with  a  resonant  wave-particle  interaction 
which  results  in  a  real  heating  of  the  ion  population,  which  we  term  an  ion 
cyclotron  resonance  heated  conic,  or  simply  ICRH  conic  for  short.  Such  dis¬ 
tributions  are  termed  conics  because  of  the  distribution  of  ions  in  their  velo¬ 
city  space:  loosely  speaking,  they  lie  on  a  cone  whose  axis  is  aligned  with 
the  mag.  -tic  *eld  direction.  Of  course,  there  are  a  wide  variety  of  ways  in 
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which  the  ions  may  be  so  distributed,  accounting  for  the  veritable  zoo  or 
conic  observations.  By  the  same  token  any  theory  which  purports  to  explain 
the  formation  of  conics  will  produce  some  characteristic  form  for  the  velocit> 
distribution.  The  point  of  the  term  ICRH  conic  is  that  these  are  events  where 
theory  and  observation  appear  to  be  talking  about  the  same  thing,  namely,  a 
conic  due  to  ion  cyclotron  resonance  heating,  and  the  term  is  used  to  distin¬ 
guish  such  events  from  others  which  must  have  some  other  explanation. 

Conics  were  first  reported  [lj  more  than  a  decade  ago,  and  have  since 
been  found  over  a  wide  range  of  altitudes  within  the  Earth’s  auroral  zone 
(see,  eg.  the  observational  reviews  by  Klumpar  [2]  and  Burch  [3]).  Conics 
are  important  because  they  provide  a  means  for  ionospheric  ions  to  escape  the 
Earth's  gravitational  field  and  consequently  constitute  a  major  component  of 
an  ionospheric  contribution  to  magnetospheric  plasma  [4-8]. 

They  are  also  interesting  to  study  because  the  ion  energization  can  be 
over  several  orders  of  magnitude  and  generally  appears  to  be  accomplished 
through  wave-particle  interaction  or  other  bona  fide  plasma  processes  (see, 
e.g.  the  theoretical  reviews  by  Lysak  [9]  and  Chang  et  al.  [10]).  Aside  from 
the  mechanisms  to  be  discussed  in  this  paper,  considerable  attention  has  been 
directed  both  towards  wave-particle  interaction  processes  as  well  as  other 
possibilities.  The  former  include  interactions  with  current-driven  electrostatic 
ion  cyclotron  waves  [11-16],  and  lower  hybrid  waves  [17-21]  which  are 
excited  by  the  energetic  precipitating  auroral  electron  beam.  Another  possi¬ 
bility  that  has  received  some  attention  is  energization  through  double  layers 
[22-24]. 

The  present  work  on  energization  through  interaction  with  electromag¬ 
netic  low  frequency  turbulence  was  motivated  by  a  conic  event  [25]  which 
did  not  seem  to  have  been  caused  by  any  of  these  mechanisms.  That  is,  the 
required  waves  were  not  observed,  and  the  conic  was  well  equatorward  of  the 
auroral  electron  beam/discrete  auroral  (or  boundary  plasma  sheet)  region 
where  they  would  be  expected.  Rather,  the  conic  was  observed  throughout  a 
broad  region  above  the  diffuse  aurora  which  maps  to  the  central  plasma  sheet 
(CPS)  of  the  magnetosphere  [25].  On  the  other  hand,  there  was  an  intense 
background  of  low  frequency  electromagnetic  turbulence  as  is  commonly 
observed  [26],  and  Chang  et  al.  [27]  considered  the  possibility  that  ions  might 
be  energized  in  much  the  same  fashion  as  heating  in  the  ion  cyclotron  range 
of  frequencies  (ICRF)  has  been  used  for  years  in  the  fusion  community. 
(ICRH  has  also  been  shown  to  be  an  important  mechanism  in  the  Earth’s  ring 
current;  see  e.g.,  the  review  by  Gendrin  [28].)  That  is.  some  remote  agency 
would  be  responsible  for  the  generation  of  electromagnetic  waves  which 
would  propagate  to  lower  altitudes  and  result  in  the  observed  ion  energiza¬ 
tion.  Indeed,  not  only  did  the  initial  estimates  demonstrate  the  viability  of 
this  scenario,  but  a  subsequent  numerical  simulation  [29]  found  that  the 
mechanism  produced  a  distribution  of  ions  remarkably  like  the  observed  con¬ 
ics. 
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Moreover,  the  heating  mechanism  works  for  other  events  as  well. 
Additional  auroral  events  [30]  have  also  been  shown  [31]  to  be  explicable  as 
ICRH  conics.  There  are  also  conic  events  in  the  cusp/cleft  where  ICRF  heat¬ 
ing  is  probably  responsible  for  the  observed  energization  of  ions  [32-34],  In 
view  of  these  exciting  developments,  this  seems  an  appropriate  time  to  review 
our  understanding  of  the  subject  and  look  to  future  discoveries. 

We  begin  in  the  next  section  with  a  primer  on  the  conic  formation  pro¬ 
cess  in  general  and  describe  the  basic  components  which  make  an  ICRH 
conic  in  particular.  Section  in  describes  a  number  of  theoretical  techniques 
that  have  been  brought  to  bear  on  this  problem  and  describes  the  type  of 
results  obtained.  These  results  are  then  compaied  with  the  experimental 
observations,  leading  into  a  brief  discussion  which  concludes  the  paper. 

II.  CONIC  FORMATION  PRIMER 

It  is  clear  that  there  are  many  ways  to  make  a  conic.  To  discuss  the 
conic  formation  process  in  its  most  general  terms,  we  offer  Fig.  1,  which 
schematically  indicates  the  various  components  which  may  play  a  role.  At 
the  bottom  of  the  diagram  we  have  indicated  a  source  population  of  cold  ions. 
These  are  the  raw  materials  for  the  conic:  typically  an  ambient,  cold,  isotro¬ 
pic  Maxwellian  plasma  with  temperatures  of  at  most  a  few  eV.  Under  the 
influence  of  various  agents,  this  population  will  evolve  through  space  and 
time  to  eventually  be  observed  as  an  energetic  conic  population  broadly 
characterized  by  an  energy  of  tens  to  thousands  of  eV,  a  perpendicular  tem¬ 
perature  in  excess  of  parallel  temperature,  and  of  course  a  conic-like  pitch 
angle  distribution. 

This  evolution  is  controlled  by  two  types  of  interactions:  eg., 
wave-panicle  interactions  which  are  responsible  for  the  energization  or  heat¬ 
ing  of  the  ion  population,  and  adiabatic  interactions  which  rearrange  the  ions 
in  their  phase  space  in  various  ways,  but  neither  add  nor  remove  energy.  We 
have  already  mentioned  several  wave-particle  interaction  processes  in  the 
Introduction;  the  adiabatic  forces  include  effects  due  to  the  geomagnetic  field 
and  the  large-scale  DC  electric  fields  found  in  the  magnetosphere.  The  rela¬ 
tive  contributions  of  these  two  types  of  interactions,  together  with  the  under¬ 
lying  geometry,  determine  the  morphology  of  the  conic  which  is  ultimately 
observed.  The  wide  variety  of  conic  observations  is  largely  due  to  the  details 
lying  behind  Fig.  1. 

Aside  from  questions  concerning  these  two  classes  of  interactions  are 
questions  that  are  in  some  sense  secondary  concerning  the  origin  of  the  wave 
turbulence  or  the  adiabatic  forces.  For  the  former,  this  can  include  the 
identification  of  a  free  energy  source,  the  manner  of  its  transformation  into 
wave  activity,  and  its  possible  propagation  to  the  region  of  the  wave-particle 
interacts  There  can  also  be  troublesome  questions  of  self-consistency 
whrr  one  extends  the  analysis  beyond  the  simple  linear  instability  phase. 
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Observed  Distribution  (Conics):  T||  <  Tj_  —  10  -  1000  eV 


Wave-Particle 


Interactions 


Evolution 
of  Ions 


Adiabatic 


Interactions 


ICRH,  LH, ... 


Source  of  Turbulence? 


Origin  of  Fields? 


Source  Distribution  (Maxwellians):  Tn  -  Tx<  1  eV 


Figure  1.  A  schematic  diagram  of  the  conic  formation  process 
whereby  low  energy  ionospheric  plasma  is  energized  through 
wave  particle  interactions  to  magnetospheric  energetics.  The 
variety  of  conic  formed  is  determined  by  both  the  type  of  wave 
turbulence  participating,  as  well  as  the  various  adiabatic  forces 
which  also  may  play  a  role. 
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Some  progress  has  been  made  on  this  question  for  the  ICRH  conics  a 
plausible  origin  for  the  wave  turbulence  has  been  identified  in  the  anisotropy 
of  the  panicle  distributions  in  the  ring  current  [35,  36],  The  waves  then  pro¬ 
pagate  to  lower  altitudes  where  they  interact  with  the  ions  to  form  conics. 

The  adiabatic  forces  are  usually  determined  by  the  placement  of  the 
conic  formation  process  within  some  larger  setting:  the  Eanh’s  dipolar  mag¬ 
netic  geometry  and  the  large-scale  potentials  of  the  auroral  acceleration  region 
and  the  polar  cap  are  determined  by  processes  which  are  indifferent  to  conics 
In  either  case  these  secondary'  questions  go  considerably  beyond  the  scope  of 
this  tutorial,  so  we  shall  not  discuss  them  further. 


A.  Basic  Physical  Ingredients 


/.  Heating  Rate.  The  first  process  to  consider  is  the  process  responsi¬ 
ble  for  the  energization  and  heating  of  the  conic.  For  ICRH  conics  this  is,  of 
course,  the  ion  cyclotron  resonance  heating  mechanism  which  we  illustrate 
schematically  in  Fig.  2.  Each  ion  of  charge  q  and  mass  m,  finds  itself  adnft 
in  a  sea  of  low  frequency  electromagnetic  turbulence.  In  general,  this 
includes  activity  spanning  a  broad  frequency  range  which  at  each  frequency 
may  be  decomposed  into  three  separate  polarizations  (left  hand,  right  hand, 
and  parallel  to  B).  To  most  of  this  activity,  the  ion  is  oblivious.  However, 
the  ion  is  sensitive  to  components  in  the  left-hand  polarization  at  frequencies 
near  to  its  cyclotron  frequency.  To  see  this,  we  need  only  integrate  the  left- 
hand  components  over  some  band  Aco  in  the  vicinity  of  the  ion  cyclotron  fre¬ 
quency  Q  (3  qB/mtc)  to  obtain  an  electric  field  E/_(r)  of  some  magnitude 
and  polarization,  EL(Cl)(&bi)iri,  which  for  some  brief  time  Ar.  rotates  at  the 
same  rate  and  in  the  same  sense  as  the  ion.  This  results  in  a  secular  change 
Av!  in  the  ion’s  velocity  vector  v,  given  by 

m,  AVi  =  qEL(t)At  ,  (1) 


and  a  consequent  change  in  perpendicular  energy 

AWX  =  ym,(vi  +Avi)  -  yrn.lvj2  =  m,  vpAv,  +  y  m,  |Avi|2  (2) 


Of  course  for  a  population  of  ions,  or  alternatively  for  a  random  sequence  of 
such  interactions,  the  term  linear  in  Avx  averages  out;  we  may  use  the  last 
term  to  construct  an  effective  heating  rate: 


W 


.  rts 


AW; 

A/ 


<72|EL(fl)l2 

2m, 


AcoA  t 


(3) 


The  interaction  time  Ar  is  limited  by  the  coherence  time  of  the  electric  field 
spectrum;  roughly  speaking,  AwA t  =  2n.  Thus  the  heating  rate  (3)  is  directly 
proporti  *1  m  the  electric  field  spectral  energy  density  evaluated  at  the  local 
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Figure  2.  A  schematic  of  the  ion  cyclotron  resonance  heating 
mechanism.  The  left-hand  side  of  the  diagram  indicates  a 
power-law  spectu^m  (log  scales)  of  electric  field  turbulence,  a 
band  near  the  cyclotron  frequency  fl  (shaded),  and  a  small  frac¬ 
tion  with  which  the  ions  are  resonant  (shaded  darker).  This  com¬ 
ponent  appears  to  each  ion  as  an  electric  field  E L(t)  which  coro¬ 
tates  with  each  ion  and  gradually  increases  its  energy. 

cyclotron  frequency,  Cl. 

2.  Diffusion.  The  preceding  discussion  of  ICRH  is  rather  heuristic 
(qualitatively  correct;  quantitatively  inadequate);  fortunately,  a  quantitatively 
correct  description  is  readily  available  from  the  standard  kinetic  theory  of 
wave-particle  interactions.  In  this  context,  the  interaction  is  completely 
described  by  a  diffusion  tensor,  whose  general  form  for  interaction  with  elec¬ 
tromagnetic  waves  has  been  worked  out  in  some  detail  [37,  38].  In  the  space 
physics  applications  which  interest  us,  we  do  not  require  such  complexity 
(and  indeed,  the  observational  data  set  is  not  currendy  sufficiently  complete  to 
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merit  it).  For  long  wavelengths  (4^,  *,v9  «  Q),  a  much  simpler  calculation 
[39]  results  in  a  diffusion  operator  acting  on  the  distribution  function  /  of  the 
form 


where 


D[f] 


r  \ 


(4) 


Dl  =  (^2/2m(2)  Jo  dtu  |£l|2((i))  8(ci)-^i,v„-Q)  (5) 

is  the  diffusion  coefficient  for  ions  of  charge  q  and  mass  ml  expressed  in 
terms  of  the  electric  field  spectral  density  in  left-hand  polarized  waves,  |£^  |2. 
Here  the  phase  space  velocity  variable  has  been  decomposed  into  components 
V||  and  parallel  and  perpendicular  to  B,  respectively,  and  we  have  retained 
only  the  most  important  term  in  expression  (4). 

The  Doppler  shift  in  Eq.  (5)  may  usefully  be  neglected:  we  have 
already  taken  Jt»vM  «  Q  and  for  minority  species  ions  there  is  typically  little 
observed  structure  in  the  electric  field  spectral  density  near  Q.  Thus  we  may 
take 

D,  =  (nq2/2m,2)\EL\2(Cl)  .  (6) 

This  results  in  the  considerable  simplification  of  having  a  velocity- 
independent  diffusion  coefficient. 

There  are  a  few  other  comments  which  should  perhaps  be  made  at  this 
point  concerning  the  applicability  of  formula  (6).  One  is  that  this  diffusion 
coefficient  depends  on  the  spectral  density  in  the  proper  (i.e.,  left-hand)  polar¬ 
ization,  rather  than  the  total  observed  turbulence.  As  we  shall  see  below,  we 
must  make  some  accommodation  for  the  fact  that  the  polarization  is  typically 
not  determined  by  the  observations.  On  the  other  hand,  the  main  study  of 
this  turbulence  [26]  did  find  evidence  for  a  characterization  of  the  turbulence 
as  at  least  partially  electromagnetic,  and  of  indefinite  (though  at  times 
predominantly  right-hand)  polarization,  implying  that  some  \EL\2  is  present  as 
some  fraction  of  the  observed  total. 

A  second  point  is  that  formula  (6)  represents  an  approximation  required 
in  part  by  our  ignorance  of  many  details  of  the  actual  wave  spectrum.  In  par¬ 
ticular,  we  are  usually  almost  totally  ignorant  of  the  k-spectrum  of  the  tur¬ 
bulence.  For  the  plasma  modes  under  consideration,  inclusion  of  such  infor¬ 
mation  should  it  become  available,  would  most  probably  result  in  corrections 
of  at  most  order  unity.  That  is  the  main  physics  and  effects  we  discuss 
would  still  transpire,  and  quantitative  results  should  still  be  valid  to  something 
like  a  factor  of  two. 

Finally,  it  is  our  intention  in  this  work  to  apply  this  heating  mechanism 
only  to  minority  species  ions.  For  the  high-altitude  auroral  and  cusp  regions, 
this  me?nv  j  and  possibly  He+,  but  probably  not  H  ions.  The  main  reason 
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is  that  cold  plasma  wave  theory  does  not  allow  waves  at  the  majority  specie?, 
cyclotron  frequency  to  propagate  in  a  left-hand  polarization,  for  precisely  the 
same  reasons  that  we  propose  such  modes  for  ion  heating.  That  is,  such 
modw  do  not  propagate  because  there  are  plenty  of  ions  to  resonate  with 
them  and  sap  them  of  energy.  The  situation  is  improved  somewhat  for 
minority  species  ions  because  they  represent  less  of  a  drain.  In  particular,  the 
‘‘stop  gap”  in  propagation  is  sufficiently  narrow  that  the  Doppler  shift  men¬ 
tioned  above  shifts  the  resonance  condition  to  neighboring,  allowed  wave 
modes.  Of  course,  cold  plasma  theory  is  only  an  approximation,  as  is  the 
infinite  and  homogeneous  geometry  it  which  it  is  derived.  The  inclusion  of 
realistic  effects  change  the  picture  somewhat.  And  of  course,  the  observa¬ 
tions  typically  show  a  structureless  spectrum  near  the  n  inority  cyclotron  fre¬ 
quency,  arguing  that  we  are  reasonably  safe  with  formula  (6). 

3.  Mirror  Force.  At  the  same  time  that  the  ions  are  being  heated,  they 
are  also  subjected  to  a  variety  of  adiabatic  (i.e.,  independent  of  heating) 
forces;  the  most  important  of  these  is  usually  the  mirror  force  which  arises 
from  the  Earth's  geomagnetic  field.  One  way  to  understand  this  force  is  to 
consider  that  in  the  absence  of  heating,  each  ion  is  constrained  in  its  motion 
to  preserve  both  its  energy  /n,  (V|2  -¥v})H  and  its  magnetic  moment  m,  vf/B 
Thus  as  an  ion  drifts  upwards  toward  regions  of  decreasing  magnetic  field  its 
vi  must  decrease  with  B .  However,  it  can  only  preserve  its  energy  with  a 
corresponding  increase  in  v„.  The  result  is  a  conversion  of  perpendicular 
energy  into  parallel  energy.  In  the  absence  of  heating,  this  would  result  in  an 
upward  moving,  field-aligned  population  at  high  altitudes.  However,  if  there 
is  heating,  the  perpendicular  energy  will  be  replenished  as  it  is  converted  into 
parallel  energy — some  form  of  conic  distribution  of  ions  then  results  as  a 
compromise  between  these  two  effects.  This  is  illustrated  schematically  in 
Fig.  3  which  shows  a  portion  of  the  ion  velocity  space.  The  effects  of  heat¬ 
ing  (horizontal  arrows)  and  the  mirror  force  (curved  arrows)  combine  to  pro¬ 
duce  a  conic  (closed  contours). 

Aside  from  it  role  in  shaping  the  form  of  the  ion  distribution,  the  mirror 
force  plays  a  secondary  role  which  is  crucial  to  the  degree  of  energization 
which  may  take  place.  Since  the  mirror  force  accelerates  ions  along  the  mag¬ 
netic  field,  the  ions  are  induced  to  travel  along  the  field  line  to  higher  alti¬ 
tudes,  and  most  importantly,  lower  gyrofrequency.  This  means  that  cyclotron 
resonance  successively  brings  the  ions  into  interaction  with  waves  of  different 
frequencies;  since  each  interaction  is  limited  to  some  brief  time  At,  the  ions 
do  not  appreciably  deplete  any  of  the  waves  they  encounter.  That  is,  the 
feedback  of  the  ions  on  the  waves  is  an  insignificant  part  of  the  story,  and  the 
process  does  not  saturate  in  the  usual  sense  of  the  term.  The  limits  to  ion 
energization  are  determined  instead  by  geometry  and  the  frequency  extent  of 
the  spectrum. 

4.  Electric  Fields.  A  second  adiabatic  force  which  can  play  a  role  in 
the  conic  formation  process  is  large-scale  electric  fields.  For  example,  in  the 
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Vi 


force  in  shaping  the  conic  distribution  (closed  contours).  Ions 
diffuse  to  larger  vL  under  the  influence  of  the  heating  mechanism 
(horizontal  arrows).  As  the  ions  drift  to  greater  altitudes  the  mir¬ 
ror  force  directs  ions  to  larger  v,  at  the  expense  of  vL  in  order  to 
conserve  energy;  the  ions  move  on  curves  of  constant  energy 
(curved  arrows  on  dotted  line). 

cusp  there  is  a  significant  cross  polar  cap  potential  leading  to  a  dawn-dusk 
electric  field.  In  combination  with  the  Earth’s  magnetic  field  this  leads  to  a 
bulk  plasma  drift,  v£  i  ExB,  which  we  must  include  in  order  to  properly 
follow  the  ions.  No  direct  energization  of  ions  results  from  this  electric  field, 
however. 

Another  possibility  is  the  field-aligned  potential  drop  such  as  is  found 
in  the  discrete  aurora.  This  results  in  an  electric  field  E\  whose  sense  is  to 
accelerate  ions  upward.  Indeed,  upward-directed  ion  beams  rather  than  conics 
are  indeed  observed  in  the  discrete  aurora.  Aside  from  the  discrete  aurora, 
however,  evidence  for  such  electric  fields  is  not  as  conclusive,  and  as  we 
shall  see  below,  ICRH  conics  can  be  adequately  formed  with  Et  a  0. 

B.  Kinetic  equation 

The  ideas  of  the  preceding  section  can  be  made  considerably  more  pre¬ 
cise  by  developing  a  kinetic  model  for  the  conic  evolutionary  process  [29, 
40-42].  Specifically,  we  are  interested  in  obtaining  the  distribution  function 
/(r,  r,  v)  n  the  usual  six-dimensional  phase  space  of  r,  v.  In  principle,  /  is 
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given  as  a  solution  of  a  kinetic  equation  of  the  general  form 

+  V-IV 1  =  Ol/l 

given  suitable  initial  and  boundary  conditions.  Here  Tj/  ]  is  the  phase  space 
flux  of  ions  due  to  adiabatic  forces  (convection,  mirror  force,  etc.),  and  D  [/  ] 
is  the  diffusion  operator  describing  the  heating  process  as  discussed  above 
Without  its  right-hand  side,  Eq.  (7)  makes  a  statement  of  conservation  of  ions 
in  phase  space.  With  the  inclusion  of  D,  this  statement  is  modified  to 
include  the  wave-induced  diffusion  of  ions. 

In  practice,  we  seldom  have  enough  information  to  merit  a  solution 
containing  the  complete  dynamics  of  the  ion  evolution  in  its  magnetospheric 
setting.  Thus  it  behooves  us  to  consider  the  process  as  more-or-less  station¬ 
ary,  and  neglect  the  time  variable  t.  Thus  we  are  more  concerned  with  the 
boundary  conditions  on  Eq.  (7)  rather  than  initial  conditions. 

In  addition,  typically  one  or  more  spatial  coordinates  are  ignorable,  and 
the  remaining  coordinate  system  is  best  aligned  with  respect  to  the  magnetic 
field  (e.g.,  use  /,  the  arc-length  along,  and  some  variable  A  across  the  mag¬ 
netic  field  line).  Additionally,  both  the  evolution  of  the  distribution  and  the 
measurement  of  the  complete  distribution  take  place  on  time  scales  longer 
than  the  gyroperiod,  so  it  is  most  convenient  to  average  Eq.  (7)  over  this 
fastest  timescale  and  work  with  a  gyrophase  average  distribution 
/(/.A,  v,,  V]^).  The  resulting  equation  for  /  is  still  rather  complicated,  and  it 
is  best  to  write  it  down  as  it  applies  to  cusp  and  the  auroral  zone  separately. 

1.  Cusp  Conics.  We  consider  the  cusp  case  first,  because  it  can  be 
made  sufficiently  simple  that  it  admits  an  exact  solution  of  the  kinetic  equa¬ 
tion.  The  typical  observational  geometry  is  illustrated  in  Fig.  4.  In  events  of 
this  type,  the  satellite  crosses  from  a  relatively  quiet  region  where  a  cold  ion 
population  is  .observed  into  a  region  of  intense,  low  frequency  turbulence 
where  a  heated  population  is  found.  At  the  same  time,  there  is  a  substantial 
E  x  B  drift  of  ions  from  the  quiet  region  into  the  turbulent  region.  In  these 
cases,  the  use  of  the  word  “conic”  is  something  of  a  misnomer  to  the  purist, 
since  the  dominant  pitch  angle  of  the  distribution  is  near  90°,  suggesting  that 
the  minor  force  does  not  play  a  major  role,  although  there  is  component  of 
drift  along  B . 

As  a  simple  example  to  analyze,  we  can  neglect  everything  but  the 
E  x  B  drift  and  the  heating.  In  this  limit,  the  kinetic  equation  (7)  becomes 

(8) 

where  vc  is  the  E  x  B  drift  speed,  and  DL  was  given  above  as  expression  (6). 
The  remaining  spatial  variable  A  is  designed  to  follow  the  motion  of  the  ions 
as  they  drift  across  the  magnetic  field  in  the  direction  of  increasing  invariant 
latitude. 
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Figure  4.  A  cartoon  illustrating  the  observational  geometry  typi¬ 
cal  of  ion  conic  events  in  the  cusp.  Here  the  satellite  drifts 
across  the  magnetic  field  from  a  relatively  quiet  region  into  a  re¬ 
gion  of  intense  electromagnetic  turbulence.  The  ions,  drifting  at 
some  angle  to  the  magnetic  held  (dashed  lines),  are  successively 
encountered  at  points  (a)  unheated,  (b)  partially  heated,  and  (c) 
significantly  heated;  contour  plot  sketches  of  the  distribution 
function  /  (A)  in  velocity  space  (vr  and  vy  are  parallel  and  per¬ 
pendicular  to  B)  at  the  corresponding  coordinates  Aa,  Ab,  and 
Ac  appear  in  the  lower  part  of  the  diagram. 


If  we  introduce  an  effective 
a  simple  diffusion  equation: 
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‘heating  time”  t(A)  =  DxA/v£,  we  obtain 


i  a 


if l  dvL 


K. 

av! 


'  9  ! 


The  ‘‘initial  condition”  is  in  fact  a  boundary  condition  at  the  edge  of  the  tur¬ 
bulent  region  which  we  identify  by  A  =  0  for  convenience.  For  simplicity, 
we  require  the  ions  to  start  at  this  point  with  a  cold,  source  Maxwellian  distri¬ 
bution,  fs,  characterized  by  a  thermal  velocity  vt.  It  is  then  straightforward 
to  obtain  the  distribution  /  (A,  vn,  vi )  at  all  values  of  A: 


/  (A,  v  (l,  vA)  = 


;c3/2vt[vf2  +  4r(A)) 


exp 


[v,2  +  4r(A)J  V[ 


(10) 


This  solution  is  sketched  in  the  lower  portion  of  Fig.  4  and  is  simply  a 
Maxwellian  population  with  a  temperature  increasing  linearly  with  A. 

In  practice,  this  simplified  model  is  only  valid  near  A  =  0;  away  from 
the  boundary,  other  effects  begin  to  play  a  significant  role.  We  shall  not  dis¬ 
cuss  the  cusp  conic  events  in  any  greater  detail  in  the  paper,  largely  since  the 
topic  has  evolved  considerably  since  the  workshop  [43,  44], 

2.  Auroral  Conics.  The  kinetic  equation  in  the  main  auroral  zone  is 
somewhat  more  complicated,  mostly  because  the  dominant  drift  of  ions  is  up 
along  the  magnetic  field  rather  than  across  it.  The  mirror  force  now  plays  a 
major  role,  consequently  terms  reflecting  its  effects  must  be  retained.  On  the 
other  hand,  we  can  neglect  terms  arising  from  electric  fields  as  these  are  now 
secondary  effects.  The  resulting  equation  may  be  written  in  the  conservative 
form  [31] 
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This  is  a  convective-diffusion  equation  for  the  density  of  ions  per  unit 
length  of  flux  tube,  f  IB,  in  the  coordinate  space  of  (/,  v„,  vj).  That  is,  the 
divergence  of  the  flow  of  ions  in  phase  space  due  to  adiabatic  effects  is  bal¬ 
anced  by  the  divergence  of  the  flow  induced  by  the  wave-particle  interac¬ 
tions.  The  components  of  the  former  are  the  field-aligned  convection  (v#)  and 
the  parallel  and  perpendicular  accelerations  of  the  mirror  force;  the  diffusive 
flow  is  proportional  to  the  velocity- space  gradient  of  f  IB  with  coefficient  0*. 

An  additional  complication  implicit  in  Eq.  (11)  is  the  dependence  of  B 
and  Di  on  /  in  a  non-crivial  fashion.  To  some  extent  we  can  simplify  this  by 
scaling  the  magnetic  field  with  altitude  roughly  as  /~3.  On  the  other  hand,  0L 


depends  on  /  directly  through  the  special  density  | EL  |2(co)  found  at  each  alti¬ 
tude,  and  indirectly  though  the  value  of  the  cyclotron  frequency  Q  at  /.  a 
simple  way  to  deal  with  this  is  first,  to  neglect  the  direct  dependence  on  /  and 
use  the  spectrum  observed  at  the  satellite  for  all  altitudes,  and  second,  to 
approximate  the  observed  spectrum  as  a  power  law  with  some  spectral  index 
a: 

\El  |2(a»  oc  oTa  .  (12) 

Then  since  the  cyclotron  frequency  varies  with  the  magnetic  field 

Cl(l)  s  ( q/mlC)B(l )  ~  (l/l0)-\  (13) 

where  we  have  used  a  convenient  reference  location  la  (e.g.  at  the  satellite)  to 
express  the  scaling,  the  diffusion  coefficient  varies  with  altitude  as 

01  =  D0(l/l0?a.  (14) 

This  model  can  be  made  slightly  more  general  than  these  assumptions  allow, 
provided  we  are  willing  to  relax  the  meaning  of  a  from  being  the  observed 
spectral  index  to  being  the  spectral  index  of  the  resonant  turbulence  (which 
spans  the  altitude  range). 

An  additional  complication  is  that  the  polarization  of  the  turbulence  is 
usually  not  determined  by  the  observations.  Typically,  what  is  observed  is  a 
spectrum  |£x  |2(co)  based  on  the  response  of  one  antenna  to  the  total  (all  polar¬ 
izations)  turbulence.  Thus  we  introduce  an  observational  form  factor  rj  which 
relates  the  turbulence  reported  to  the  fraction  in  the  left-hand  polarization: 

\el\2  =  nl£,l2. 

The  situation  is  illustrated  schematically  in  Fig.  5.  Here  we  illustrate 
the  evolution  of  the  ions  from  low  altitudes  (bottom  of  the  figure)  to  high 
altitudes  within  some  flux  tube.  The  boundary  condition  applied  to  Eq.  (11) 
is  that  the  low  altitude  population  be  characterized  by  ionospheric 
parameters — typically  a  cold  Maxwellian.  At  any  altitude  / ,  the  ions  resonate 
with  the  turbulence  at  their  local  gyrofrequency,  gain  a  bit  of  energy,  convert 
some  perpendicular  energy  into  parallel  energy,  and  drift  to  higher  altitudes. 
This  sequence  is  repeated  at  successively  higher  altitudes  as  the  ions  evolve 
into  an  energetic  conic.  We  shall  examine  this  scenario  in  greater  detail  in 
the  next  section  as  we  discuss  different  techniques  for  extracting  information 
from  Eq.  (11). 


ID.  SOLUTION  TECHNIQUES 

In  this  section  we  shall  discuss  some  of  the  different  approaches  which 
have  succeeded  in  contributing  to  our  understanding  of  ICRH  conics.  Each 
approach  has  its  respective  merits  and  limitations.  In  concert,  however,  we 
obtain  a  wealth  of  information. 
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Figure  5.  A  cartoon  of  the  evolution  of  ICRH  conics  confined  to 
an  auroral  flux  tube.  Different  portions  of  the  electric  field  spec¬ 
tral  density  (shaded  bands  in  plot  at  the  left)  interact  with  the 
ions  at  successive  altitudes  (wave-particle  interaction  cartoon  at 
right)  to  produce  a  gradual  evolution  of  the  ion  distribution  func¬ 
tion  from  cold  distribudons  (at  the  center  bonom)  to  hot,  energet¬ 
ic  distributions  (at  the  center  top).  Following  an  initial  relaxation 
(not  shown,  but  somewhere  below  the  figure),  the  form  of  the 
conic  remains  unchanged. 
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A.  Numerical  Solution 


A  fairly  direct  approach  is  look  for  a  numerical  solution  of  the  problem, 
and  this  is  readily  obtained  using  the  Monte  Carlo  solution  technique  [29], 
which  was  first  applied  to  the  formation  of  lower  hybrid  conics  [45],  The 
idea  here  is  to  follow  some  number  of  “test”  panicles  through  their  journey 
up  the  flux  tube.  The  action  of  the  wave  turbulence  is  mimicked  by  the 
application  of  a  sequence  of  random  velocity  kicks  Av  which  result  in  velo¬ 
city  space  diffusion.  An  observed  spectrum  of  turbulence  may  be  invoked  to 
determine  the  strength  of  the  diffusion  via  Eq.  (6);  the  principle  uncenainties 
here  are  the  fractional  polarization  parameter  rj  and  the  exact  form  of  the  ini¬ 
tial  distribution. 

This  technique  was  applied  to  the  conic  event  observed  on  the  Dynam¬ 
ics  Explorer  1  (DE-1)  satellite  with  the  High  Altitude  Plasma  Instrument 
(HAPI)  reported  in  Ref.  [25];  the  results  are  shown  in  Fig.  6  [29],  The 
simulation  used  a  spectral  index  a  =  1.7  and  intensity  |£z|2  =  1.2xl(T* 
V2/m2Hz  at  fci  (=fl/2jt)  =  5.6  Hz,  derived  from  the  simultaneous  DE-1 
Plasma  Wave  Instrument  (PWI)  measurements.  The  mean  particle  calcula¬ 
tions  [27]  indicated  a  value  of  r\  -  1/8  was  appropriate.  The  temperature  of 
the  initial  distribution  was  taken  -  0.2  eV,  which  is  appropriate  for  0+  ions 
at  1.2  Re  geocentric.  No  field- aligned  potential  drop  was  required  to  produce 
the  net  upward  drift  of  the  conic. 

One  of  the  interesting  results  of  the  simulation  (aside  from  the  excellent 
agreement  with  the  observations)  is  the  insensitivity  of  the  result  to  the  initial 
distribution.  Basically,  the  cooperative  action  of  the  diffusion  and  mirror 
forces  is  sufficiently  strong  to  rapidly  rearrange  the  initial  form  of  /  into  that 
of  a  conic.  Indeed,  in  the  simulation,  the  ion  distribution  was  recognizably  a 
conic  after  only  a  several  hundred  km  increase  in  altitude.  Thereafter,  the 
form  of  the  distribution  remained  relatively  fixed  while  becoming  more  ener¬ 
getic  with  increasing  altitude.  This  sort  of  evolution  was  included  in  the 
sketch  of  the  distributions  in  Fig.  5,  is  termed  self-similar,  and  is  discussed  in 
the  next  subsection. 


B.  Similarity  Regime 

The  essential  point  we  shall  exploit  here  is  the  fact  that  the  high  alti¬ 
tude  conic  observations  are  made  in  what  can  be  termed  a  similarity  regime 
where  the  conic  evolution  is  self-similar.  That  is,  the  observations  are  made 
above  the  altitudes  where  the  distribution  is  in  the  process  of  relaxing  onto 
the  conic  shape  and  still  retains  some  of  its  initial  character.  This  gives  us 
the  freedom  to  ignore  details  of  the  initial  (unknown)  distribution  fs  because 
they  are  irrelevant.  The  situation  is  entirely  analogous  to  the  process  of  strik¬ 
ing  a  tuning  fork  to  produce  a  tone.  After  a  brief  interval,  one  obtains  the 
same  tone  ndependent  of  how  the  tuning  fork  was  struck. 
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Figure  6.  The  top  panel  is  the  result  of  a  Monte  Carlo  calcula¬ 
tion  for  the  oxygen-dominated  conic  which  was  observed  on  the 
DE-1  satellite  with  the  HAPI  instrument  and  is  displayed  in  the 
bottom  panel.  The  contours  in  both  panels  are  uniformly  spaced 
with  an  increment  of  0.4  in  the  logarithm  of  the  phase  space  den¬ 
sity  (from  Ref.  [29]). 

1.  Distribution  Function.  The  physical  basis  for  the  similarity  regime  is 
that  at  each  altitude  throughout  the  evolution  process  the  physical  conditions 
are  in  some  sense  the  same — the  evolution  of  the  distribution  at  any  altitude 
is  similar  to  the  evolution  at  other  altitudes.  Formally,  some  of  the  physical 
quantities  (5  and  D{)  vary  as  powers  of  (///„),  and  it  is  possible  to  scale  the 
rest  in  the  same  fashion,  leaving  the  kinetic  equation  (11)  invariant.  One  way 
to  do  this  is  to  introduce  new  velocity  variables  x  and  y  for  the  perpendicular 
and  parallel  components  of  velocity,  respectively  [41]: 
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*(vit/)  3  (Vi /v„  )(///„  r° . 


(15aj 


y(vt,l)  s  (v#  /va  )(///„ )~°  , 


(15b; 


and  ask  that  the  distribution  function  depend  on  velocity  and  l  only  through 
these  variables,  except  perhaps  for  some  overall  factor  of  (l/l0 ).  The  parame¬ 
ters  O  and  v0  give  us  two  degrees  of  freedom  with  which  to  make  this 
approach  work.  In  fact,  if  we  take 


o 


3a  +  1 
3 


(16a) 


=  (D0L) 


\n 


(16b) 


we  find  that  the  shape  of  the  conic  is  then  given  by  a  dimensionless  distribu¬ 
tion  F(x,  y)  related  to  /  (/ ,  v«,  vL )  by 

/(/,vfcVl)  =  n((/<,)v„'3F(Jt(vi,/),y(vll,/))(///(,r{4<,+3).  (17) 

Here  the  first  two  factors  carry  the  dimensions  of  the  distribution:  «,  (/„)  is 
the  density  of  ions  at  the  observation  point,  va  plays  the  role  of  a  characteris¬ 
tic  thermal  velocity,  and  the  final  factor  ensures  that  ions  are  conserved  in  the 
flux  tube.  At  any  given  altitude,  the  distribution  F(x,y)  is  the  same  as 
/(/,  vii,  vL )  with  the  velocity  scales  relabeled  according  to  Eqs.  (15). 

The  equation  which  determines  F  may  be  usefully  cast  as  a 
convective-diffusion  equation  for  a  density  N(x,  y)  =  x  F(x,  y): 

V-(uN)  =  ^  ,  (18) 

9x 


where  there  is  an  imposed  flow  field 


J_ 

x 


*y  . 


(19a) 


This  flow  results  from  the  combined  effects  of  the  magnetic  mirror  force,  the 
mean  particle  heating,  and  the  scaling  of  velocities.  The  random  component 
of  the  diffusion  process  is  described  by  the  right  hand  side  of  Eq.  (18). 

An  important  point  to  note  concerning  these  equations  is  that  they  are 
parameterized  only  by  o.  The  parameter  va  enters  only  into  the  determina¬ 
tion  of  the  velocity  scale  and  the  dimensional  factors.  The  situation  is  analo¬ 
gous  to  the  characterization  of  a  bi- Maxwellian  distribution  by  a  thermal 
velocity  (vf )  and  an  anisotropy  factor  (A  *  T/T 0.  The  energy  of  the  distri¬ 
bution  is  characterized  by  vf,  but  its  shape  is  completely  described  by  the 
anisotropy  tctor  A . 
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The  Monte  Carlo  approach  may  be  applied  effectively  to  the  solution  of 
this  equation  [41}.  The  central  idea  is  to  treat  Eq.  (18)  as  the  time-asymptotic 
limit  of  a  Fokker-Planck  equation.  One  may  then  launch  a  single  particle 
(with  essentially  arbitrary  initial  conditions)  and  follow  its  evolution  accord¬ 
ing  to  Langevin  equations.  The  distribution  N(x,  y)  is  then  available  as  the 
time  integral  of  the  particle  distribution  divided  by  the  duration  of  the  simula¬ 
tion.  From  N(x,  y)  we  may  construct  F(x,  y)  through  division  by  x .  This 
procedure  is  considerably  more  efficient  than  the  more  general  Monte  Carlo 
procedure  discussed  above,  which  must  retain  all  of  the  details  of  the  altitudi¬ 
nal  evolution. 

This  procedure  is  illustrated  for  a  =  2.53  in  Fig.  7.  Several  streamlines 
of  the  flow  field  u(x,  y)  are  indicated  in  the  left-hand  panel  of  Fig.  7.  The 
first  thing  to  note  is  that  the  flow  converges  on  a  fixed  point  where  u  =  0 
This  occurs  at  the  point 
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The  convergence  that  this  flow  pattern  exhibits  is  then  balanced  by  the 


Figure  7.  The  left-hand  panel  shows  the  distribution  N(x,  y)  for 
a  =  2.53  obtained  by  integrating  the  Langevin  equations  for  a 
panicle  in  a  flow  field  u(x,y),  indicated  by  the  dotted  trajec¬ 
tories.  The  distribution  F(x,y)  obtained  from  N(x,y )  appears 
in  the  right-hand  panel.  The  contours  in  these  plots  are  placed  at 
half  decades  (from  Ref.  [41]). 
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diffusive  action  of  the  right-hand  side  of  Eq.  (18).  Thus  as  the  density  tends 
to  pile  up  at  the  fixed  point,  the  action  of  this  term  is  to  spread  out  the  pile 
into  something  of  a  mound.  Although  the  diffusion  acts  only  in  the  x  direc¬ 
tion,  the  flow  pattern  produces  a  spread,  albeit  smaller,  in  the  y  direction 
This  is  seen  in  the  form  of  the  function  N(x,  y)  also  included  in  the  left-hand 
panel  of  Fig.  7.  Division  by  x  then  skews  the  result  to  obtain  the  distribution 
F(x,  y).  One  can  determine  that  N(x~* 0,  y)  goes  linearly  to  zero,  with  the 
result  that  F  maximizes  at  x  =  0.  The  resulting  distribution  F  (x ,  y )  appears 
in  the  right-hand  panel  of  Fig.  7. 

We  note  that  this  distribution  has  all  of  the  qualitative  features  of  the 
observed  conic.  One  feature  of  special  note  is  the  absense  of  ions  at  small  v , 
'i.e.  we  find  F  (x ,  y -*0)  -*  0  rapidly).  This  is  the  feature  that  invites  expli¬ 
cation  in  terms  of  a  field-aligned  potential  [30],  but  here  we  clearly  see  how 
the  feature  arises  without  it:  the  flow  u(jc  ,  y )  sweeps  particles  away  from  the 
origin  toward  the  fixed  point,  and  there  is  no  way  for  the  ions  to  diffuse  back 
to  the  origin. 

2.  Moments  Further  progress  can  be  made  with  the  theory  in  the  simi¬ 
larity  regime.  So  far,  we  have  concentrated  on  methods  (albeit  numerical) 
which  may  be  used  to  obtain  the  distribution  function.  However,  this  func¬ 
tion  often  contains  more  information  that  we  need.  In  many  situations,  one  is 
quite  happy  to  make  due  with  less  information,  and  this  is  often  expressed  in 
terms  of  some  fluid  moments  of  the  distribution: 

<  vf  vT>  3  2  rt  J0**VidV|  |o"dv„  vf  vy/a.VfcVi)  .  (21) 

Analogously,  we  may  use  the  scalings  (15)  to  obtain  dimensionless  variants 
<xnym>.  Typically,  we  are  interested  in  a  few  of  the  lowest  order 
moments,  as  these  tell  us  the  density,  momentum,  energy,  and  energy  flux  of 
the  conic. 

The  usual  road  to  the  calculation  of  moments  is  to  construct  fluid  equa¬ 
tions  from  the  kinetic  equation  (11).  Unfortunately  for  our  case,  all  this  does 
is  to  generate  an  infinite  hierarchy  of  equations  with  no  obvious  “equation  of 
state’’  which  may  be  used  to  correctly  truncate  the  hierarchy.  Applying  the 
same  procedure  to  the  velocity-scaled  equation  for  F  is  not  of  much  more 
help  as  the  result  is  a  three  term  recursion  relation  among  the  moments  of  the 
distribution.  Re-expressed  in  terms  of  the  physical  velocities,  this  relation  is 

3m  <  vf+2  vf"1  >  -  [/i(3+2o)+2mo]  <  vfvy*1  > 

(22) 

+  2n2v03<vf-2vr>  =  0  . 

Unfortunately,  the  solution  to  this  relation  has  infinitely  many  free  parameters 
(e.g.,  <  vf  >  and  <  vf  >  with  integers  n  >  0  even  and  m  >  0).  Neverthe¬ 
less,  it  turns  out  that  it  is  possible  to  obtain  these  moments  of  the  distribution 
directly  via  a  path  integral  technique  [4^j. 


21  . 


The  details  would  take  us  considerably  beyond  the  scope  of  this  paper 
However,  the  basic  idea  is  that  a  probability  density  functional  P  is  available 
to  assign  relative  probabilities  to  any  path  'P(r)  the  ion  might  take  in  obeying 
the  Langevin  equations  for  a  single  particle  which  may  be  derived  from  the 
kinetic  equation  (18).  An  example  of  such  a  path  is  the  one  actually  taken  by 
the  single  particle  we  followed  in  constructing  Fig.  7.  There  are  infinitely 
many  such  paths,  and  their  relative  probability  ¥(;)  can  be  expressed  as  a 
functional  of  position  (.t ,  y )  and  conjugate  momenta  (pz ,  py)  components  of 
the  path  'F(r), 

PpF]  =  exp {i  jdt[px(ux-x)  +  py(uy-y)  +  ip?x]  }  .  (23) 

Arbitrary  moments  M  C¥)  of  the  distribution  may  then  be  constructed  as 
weighted  path  integrals: 

fD*F  M{'¥)P['¥] 

<M0F)>  =  1 — -  .  (24) 

Jd*F  PpF] 

It  turns  out  that  evaluation  of  these  moments  is  practical  via  a  Feynman 
diagram  perturbation  expansion  technique.  The  unperturbed  system 
corresponds  to  an  ion  oscillating  about  the  fixed  point  of  the  flow  u.  The 
effective  perturbation  parameter  is  the  relaxation  rime  for  the  ion  to  return  to 
this  fixed  point;  this  time  decreases  with  increasing  o.  It  is  then  practical  to 
obtain  arbitrary  moments  of  the  conic  by  algebraic  means. 

Results  of  the  calculation  for  the  parallel  (n,  v02<x°y2>)  and  perpen¬ 
dicular  (ntvj<xz y°>)  energy  densities,  and  the  parallel  (n,vo3<x0y3>)  and 
perpendicular  (n, v03< z*y *>)  energy  flux  densities  are  shown  in  Fig.  8. 
While  the  parallel  quantities  vary  with  o  by  over  an  order  of  magnitude  for 
the  range  plotted,  the  perpendicular  quantities  are  relatively  insensitive  to  a. 
This  is  due  to  the  fact  that  smaller  values  of  a  correspond  to  effectively 
larger  parallel  drifts,  while  at  the  same  time  the  relative  strength  of  the 
diffusion  is  roughly  independent  of  o.  The  solid  points  (•)  correspond  to 
numerical  solutions  of  the  kinetic  equation  (18).  The  zeroth  and  first  order  of 
the  perturbation  expansion  are  indicated  by  dashed  ( — )  and  solid  ( — )  lines, 
respectively. 

A  number  of  features  are  worthy  of  note.  One  is  the  fact  that  the 
expansion  improves  with  increasing  c.  This  is  due  to  the  fact  that  at  large  c 
particles  are  more  rapidly  swept  into  the  sink  of  the  flow  u.  The  result  is  a 
comparatively  tighter  conic.  Conversely,  as  o  decreases,  the  turbulence  is 
able  to  better  scatter  the  ions  against  this  flow,  resulting  in  a  greater  dispersal 
of  ions  in  phase  space — a  distribution  that  might  appear  to  be  a  “heated 
beam”  rather  than  a  conic. 

A  final  point  to  note  is  that  the  results  are  not  terribly  sensitive  to  a. 
One  could  interpret  this  to  mean  that  formation  of  a  conic  is  relatively 
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Figure  8.*  Results  of  the  path  integral  calculation  for  the  parallel 
(n,  v„2<  x  V>)  and  perpendicular  x2y0>)  energy  densi¬ 
ties,  and  the  parallel  and  perpendicular 

(/i,v03<xy>)  energy  flux  densities.  The  solid  points  (•) 
correspond  to  numerical  similarity  regime  solutions.  The  dashed 
( — )  and  solid  ( — )  lines  are  the  result  of  the  zeroth  and  first  ord¬ 
er  of  perturbation  theory,  respectively  (from  Ref.  [31]). 

insensitive  to  the  spectral  index  of  the  applied  turbulence.  Alternatively,  a 
determination  of  the  spectral  index  from  the  moments  of  the  conic  distribu¬ 
tion  is  not  likely  to  be  an  accurate  procedure.  On  the  other  hand,  the 
moments  nevertheless  contain  substantial  information  concerning  the  particle 
distribution  and  indeed  could  be  used  to  categorize  particle  observations,  e.g., 
as  an  extension  of  the  current  use  of  densities,  drifts  and  temperatures. 


IV.  THEORY  MEETS  OBSERVATIONS 


With  the  preceding  theoretical  discussion  behind  us,  we  are  now  in  a 
position  to  consider  how  well  the  theory  meets  the  observations.  In  fact,  ue 
shall  see  that  the  agreement  is  rather  good.  To  date,  three  conic  observations 
have  been  given  a  detailed  examination  within  the  context  of  this  theory  [29. 
31].  An  overview  of  these  events  appears  in  Table  I.  These  events  were 
observed  with  the  DE-1  satellite  within  the  (diffuse)  auroral  zone.  A  gross 
characterization  of  each  event  may  be  made  with  the  energies  a  m,v,2/2. 
Comparable  values  for  each  event  are  also  obtained  from  the  wave  data  using 
the  heuristic  theory  [27]. 


Table  I.  Three  0+  ICRH  conic  events  observed  with  the 
Dynamics  Explorer  1  satellite  (from  Ref.  [31]). 


81288$ 

81318t 

82061$ 

Year 

1981 

1981 

1982 

Day 

Oct.  15 

Nov.  14 

Mar.  2 

UT 

2046 

2346 

1743 

h  min 

lo 

4.28 

2.0 

4.52 

Re 

IL 

73.7 

60.0 

66.7 

0 

MLT 

811 

1911 

000 

h  min 

Kp§ 

3- 

5 

5+ 

fa 

0.67 

5.6 

0.43 

Hz 

|£,I2(/c.) 

-  1.5xlO"5 

-  1.2X10-6 

-8.8X10"6 

V2m'2Hz~l 

-0.3 

-  10 

-  1.4 

cm"3 

-  1020 

-40 

-240 

eV 

-850 

-30 

-  190 

eV 

$Cf.  Ref.  [30];  tCf.  Ref.  [25];  §from  Refs.  [46-48]. 


Here  we  shall  review  only  the  first  event  (day  81288)  in  detail  [31]. 
This  event  has  received  a  fair  amount  of  attention  in  the  literature,  beginning 
with  its  initial  publication  [30],  where  the  form  of  the  distribution  was  taken 
to  be  suggestive  of  a  two- stage  acceleration  process.  The  bulk  parallel  drift 
of  the  conic  was  supposed  to  be  the  result  of  passing  a  heated  distribution 
through  a  field-aligned  potential  drop.  Subsequent  explanations  for  the  distri¬ 
bution  were  offered  [49,  50]  in  which  the  necessary  acceleration  arose  out  of 
the  transport  of  ions  across  the  polar  cap.  However,  as  we  have  seen  in  the 
HAPI  event  (Fig.  6,  and  Ref.  [29]),  this  apparent  parallel  acceleration  arises 
naturally  out  of  the  wave-particle  interaction  process,  and  indeed  no  parallel 
electric  field  is  required.  Indeed,  Temerin  [51]  showed  that  this  effective 
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parallel  acceleration  arises  from  any  diffusive  acceleration  process  occurring 
over  a  wide  altitude  range  in  the  Earth’s  magnetic  geometry  and  was  the  first 
to  apply  such  an  argument  to  the  81288  event.  Indeed,  preliminary  work  [52] 
demonstrated  that  the  ICRH  mechanism  was  viable.  We  review  here  some  of 
the  subsequent  work  [31]  which  supports  the  view  that  the  event  is  most 
probably  an  ICRH  conic. 


A.  Observations 


I.  Wave  Data.  As  we  have  seen  in  the  theoretical  analysis,  the  charac¬ 
teristics  of  ICRH  conics  are  completely  determined  by  the  two  parameters  o 
and  v0.  We  therefore  need  to  extract  these  parameters  from  the  observational 
data. 

Owing  to  the  high  altitude  of  the  observations,  the  oxygen  cyclotron 
frequency  fell  below  the  normal  frequency  range  of  the  Low  Frequency 
Correlator  (LFC)  of  the  Plasma  Wave  Instrument  (PWI)  [53]  and  it  was 
necessary  to  obtain  the  electric  field  spectral  density  from  the  DC  electric 
field  portion  of  the  instrument.  The  spectral  information  could  then  be 
obtained  from  a  Fourier  transform  of  each  spin’s  worth  (6  s)  of  data — in 
practice  we  obtained  smoother  spectra  using  the  maximum  entropy  method 
(MEM)  which  has  been  shown  to  do  a  better  job  representing  spectra  of  the 
power-law  type  [54], 

A  detailed  examination  of  the  DC  wave  forms  [31]  also  suggested  that 
the  left  and  right  wave  polarizations  were  dominant  over  the  parallel  polariza¬ 
tion.  Thus  in  using  the  approximation  \EL  |2  =  t\\Ex  |2  we  note  that  unpolarized 
turbulence  (\EL  |2  -  IE#  |2)  corresponds  to  r|  -  2. 

The  electric  field  data  for  the  event  is  presented  in  rig.  9.  Here  the 
spectra  is  plotted  in  a  time  series  format,  using  a  combination 
contour/halftone  representation.  What  is  immediately  obvious  from  the  figure 
is  that  the  peak  of  wave  activity  occurs  near  the  center  of  the  figure  (-  UT 
2046:40  =  UT  74800  s).  This  is  in  some  respects  analogous  to  the  sequence 
observed  in  the  energy-time  spectra  of  the  EICS  data  (left-hand  side  of  the 
top  panel  Plate  l  of  Ref.  [30]). 

With  these  spectra  in  hand,  it  is  rather  straightforward  to  compute  the 
parameters  a  and  v0.  In  fact,  c  is  determined  from  the  spectral  index 
(modulo  uncertainties  in  the  model)  by  Eq.  (16a),  and  v0  is  determined  by  the 
electric  field  spectral  density  at  the  resonant  frequency  via  Eqs.  (6)  and  (16b); 
in  SI  units  we  have 


v 


O 


38.52 


lo 

1/3 

THEZ|2(/C() 

1  Rg 
i  c  j 

.  10-6  VV^Hz'1 . 

1/3 

(km/s) 


(25) 


where  fci  *fl(/0  )/2x  is  the  local  local  cyclotron  frequency  in  Hz. 
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Figure  9.  Electric  field  spectral  density  based  on  the  DE-1  PWl 
DC  data  plotted  as  a  function  of  UT  on  day  81288.  A 
contour/halftone  format  is  used,  where  the  contours  are  logarithm¬ 
ically  spaced  by  decades,  and  the  darker  regions  correspond  to 
the  greater  spectral  density.  The  vertical  lines  near  UT  74900  s 
result  from  a  data  gap  (from  Ref.  [31]). 

2.  Particle  Data.  The  81288  conic  is  sufficiently  energetic  to  be  well 
resolved  by  the  Energetic  Ion  Composition  Spectrometer  (EICS)  [55]  instru¬ 
ment  in  the  “FAST”  mode  in  which  it  was  operating.  In  order  to  construct  a 
distribution  function,  each  individual  measurement  (tested  for  significance 


against  the  noise  background)  was  used  to  provide  an  estimate  of  fQ  within 
ranges  of  pitch  angle  and  velocity  as  determined  by  the  angular  and  energv 
acceptances  of  the  instrument  together  with  the  duration  of  the  measurement. 
Over  time,  values  of  f  0  become  available  for  a  set  of  points  distributed  on  a 
plane  in  velocity  space  determined  by  the  satellite  rotation  axis.  Assuming  a 
reasonably  continuous  distribution,  values  anywhere  on  the  plane  may  be 
obtained  through  interpolation  or  simply  reference  to  the  nearest  measured 
value.  In  the  FAST  mode,  alternate  6  s  satellite  spin  periods  are  devoted  to 
different  ion  species  and  different  (interleaved)  energy  acceptance  ranges. 
Thus,  a  partial  energy  scan  is  available  after  6  s,  while  1 8  s  are  required  for  a 
complete  energy  scan  of  one  species. 

However,  a  transformation  from  this  plane  of  measurement  to  the  space 
of  parallel  and  perpendicular  velocity  components  is  still  required.  Given  the 
relative  drift  between  the  satellite  frame  and  the  plasma  frame  (where  it  is 
assumed  to  be  gyrotropic),  the  required  transformation  is  easily  constructed. 
For  the  events  considered,  the  plane  of  measurement  nearly  includes  the 
field-aligned  direction,  so  each  half  of  this  plane  spans  most  of  the  (v^  ,v  „ ) 
space.  There  are  thus  a  variety  of  ways  to  use  the  data  to  construct 
/0(vi»vn):  the  differences  between  them  may  be  taken  as  an  uncertainty  of 
the  determination  of  fa,  given  the  statistics  of  the  measurements  and  other 
assumptions  about  the  plasma  conditions  which  may  be  only  approximately 
satisfied. 


Table  n.  Moment  relations  that  determine  the  parameters  a  and  va . 
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We  determine  o  and  v0  from  the  EICS  observation  of  fa  using  the 
recursion  relation  (22)  as  illustrated  in  Table  II.  In  particular,  we  note  that 
each  instance  of  the  recursion  relation  (22)  presents  us  with  a  linear  relation 
in  these  parameters.  Any  two  of  these  relations  may  be  solved  to  yield  their 
values,  but  we  would  not  expect  different  pairs  to  result  in  the  same  values 
unless  f0  is  truly  an  ICRH  conic  and  all  the  approximations  are  applicable. 
In  practice,  we  obtain  o  and  va  from  a  least  squares  fit  to  the  three  relations 
given  in  Table  II.  Some  appreciation  of  the  relative  error  involved  may  be 
obtained  by  considering  the  spread  of  the  solutions  to  these  relations  taken 
pairwise. 
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B.  Comparison 


1.  Distribution  Functions.  We  first  address  the  extent  to  which  the 
conic  observed  is  consistent  with  the  form  predicted  by  theory.  We  would 
expect  to  see  the  most  energetic,  and  therefore  best  resolved,  conic  near  the 
peak  of  the  event  determined  from  the  wave  activity  (UT  74800  s  from  Fig. 
9).  In  order  to  reduce  the  level  of  noise  in  the  data,  an  average  of  three  com¬ 
plete  energy  scans  (UT  74794-74836  s)  is  presented  in  Fig.  10. 

For  comparison  with  the  numerical  solutions  for  F  (* ,  y ),  we  compute 
the  moments  of  this  distribution,  finding  a  density  n,  (/0)  =  0.3  cm-3,  v0  ~ 
162±5  km/s,  and  a  =  2.76±0.2.  Thus  (in  round  figures)  we  apply  a  velocity 
scale  va  =  160  km/s  to  the  numerical  solution  obtained  with  a  -  8/3  and 
present  the  result  in  Fig.  11.  Allowing  for  the  artifacts  of  the  data  acquisition 
and  reduction  (e.g.,  the  “boxes”  in  Fig.  9  which  correspond  to  the  individual 
data  points),  the  agreement  between  Figs.  9  and  10  is  remarkable:  both  dis¬ 
tributions  peak  near  v*  -  80  km/s  as  — *  0,  and  the  V|  extent  of  the  con¬ 
tours  is  comparable.  There  are  several  measurements  in  the  observational 
plot  (e.g.,  for  vx  -  220-300  km/s)  which  have  no  counterpart  in  the  theoreti¬ 
cal  distribution;  these  are  probably  pan  of  an  unrelated  oxygen  population. 

2.  Moments.  We  next  turn  to  the  time  series  of  the  observations.  This 
is  most  conveniently  presented  in  terms  of  va  and  o  which  may  be  deter¬ 
mined  separately  from  each  instrument,  and  which  arc  presented  in  Figs.  12 
and  13.  Indeed,  an  important  reason  for  using  these  parameters  is  that  the 
observations  of  each  instrument  may  be  independently  reduced  to  these 
parameters,  thereby  testing  the  theory  which  says  they  should  produce  similar 
results. 

The  determination  of  v0  and  o  with  the  EICS  observations  can  be  made 
on  a  spin  by  spin  basis.  Two  six-second  spin  periods  separated  by  a  spin 
period  (or  a  total  of  18  seconds)  are  required  to  construct  a  complete  O 
energy  distribution  which  includes  left  and  right  halves;  thus  Fig.  12 
expresses  the  results  obtained  by  various  ways  in  order  to  express  some  of 
the  uncertainties  inherent  in  this  calculation.  In  particular,  there  are  dotted 
and  dashed  lines,  which  correspond  to  the  right  and  left  halves  of  the  conic, 
respectively,  and  these  traces  are  doubled  because  there  are  two  ways  to  pair 
adjacent  spins.  The  sum  of  the  left  and  right  halves  is  plotted  as  points 
(squares  and  triangles  denoting  the  pairing  of  spins),  with  vertical  error  bars 
ba&d  on  the  uncertainty  in  the  determination  of  v0  and  a  as  described  above 
and  horizontal  error  bars  indicating  the  interval  of  measurement  (18  s). 

An  immediate  impression  from  Fig.  12  is  that  the  determination  of  c  is 
not  nearly  as  precise  as  that  of  v0 .  This  shows  up  in  the  error  bars  as  well  as 
the  general  scatter  of  the  data  points.  The  primary  explanation  for  this  is  the 
fact  that  although  <3  determines  the  shape  of  the  conic,  its  control  is  some¬ 
what  subtle.  While  there  is  considerable  qualitative  difference  between  a 
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Figure  10.  Conic  distribution  f0(vltv,)  constructed  from  four 
spins  of  DE-1  EICS  data  near  UT  7480C  s  on  day  81288.  Con¬ 
tours  are  placed  at  half  decades,  and  the  darker  regions  denote  the 
greatest  phase  space  density  (from  Ref.  (31)). 

conic  with  o  =  0  and  0=1,  there  is  little  difference  between  o  =  2  and 
o  =  5.  An  excellent  way  to  see  this  is  via  the  path  integral  results  (Fig.  8) 
which  show  a  general  insensitivity  of  the  moments  to  o  values  beyond  2. 

Another  observation  to  be  made  from  Fig.  12  is  that  the  error  bars  for 
times  before  UT  74780  s  are  qualitatively  larger  than  those  after  this  time. 
One  may  .nterpret  this  as  a  suggestion  from  the  data  that  the  observed 
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Figure  11.  Conic  distribution  f0(yi  ,V| )  constructed  from 
F(*,y)  with  o  =  8/3  and  a  velocity  scale  of  v0  =  IbGkm/s. 
The  representation  is  identical  to  that  of  Fig.  10  (from  Ref.  [31]). 


distribution  here  is  far  from  an  “ideal”  conic.  Indeed,  examination  of  the 
detailed  distributions  at  this  time  show  a  contamination  with  a  quite  obviously 
unrelated,  energetic  0+  population. 

The  spectra  obtained  from  the  PWI  DC  data  were  averaged  over  three 
spins  to  provide  a  time  resolution  comparable  to  the  panicle  data;  low  time 
resolution  (32  s)  PWI  LFC  data  were  also  used  for  comparison.  Either  way, 
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Figure  12.  Time  series  of  v0  and  o  on  day  81288  based  on 
DE-l  EICS  observations.  In  each  panel,  the  dotted  and  dashed 
lines  correspond  to  the  right  and  left  halves  of  the  conic,  and  the 
solid  points  with  error  bars  correspond  to  their  average.  Squares 
and  triangles  distinguish  the  two  ways  in  which  pairs  of  spins 
may  be  composed  to  construct  a  complete  energy  scan;  these  two 
ways  are  undifferentiated  for  the  lines  (from  Ref.  [31]). 

a  v0  and  o  time  series  is  then  directly  obtained  from  Eqs.  (25)  and  16. 
Because  the  polarization  is  undetermined,  we  have  presented  v0  values 
corresponding  to  a  range  of  q1/3:  1.0,  0.8,  0.6,  and  0.4  for  the  DC  data,  and 
0.8  and  0.  j  for  the  LFC  data.  These  last  two  values  correspond  to  between 
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Figure  13.  Time  series  of  v0  and  o  on  day  81288  based  on 
PWI  observations.  In  each  panel,  the  lines  correspond  to  power- 
law  fits  of  the  DC  data  and  the  open  points  to  power-law  fits  of 
the  LFC  data.  The  multiplicity  of  v0  values  comes  from  the  use 
of  different  values  of  rj1/3:  1.0  (dotted),  0.8  (solid,  squares),  0.6 
(solid,  triangles),  0.4  (dotted;  from  Ref.  [31]). 

13%  and  3%  of  the  total  electric  field  spectral  density  being  left  hand  polar¬ 
ized.  We  note  that  the  general  trend  of  the  event  is  consistent  with  Fig.  12, 
and  that  the  DC  and  LFC  data  are  in  general  agreement.  The  DC  spectra  at 
this  timescale  (18  s)  are  somewhat  more  dynamic  than  the  LFC  data,  which 
have  a  lower  dme  resolution  (32  s).  There  appears  to  be  some  trend  for  the 
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LFC  data  to  provide  somewhat  more  energetic  values  of  v0\  on  the  other 
hand,  at  these  altitudes,  the  lowest  LFC  frequency  band  lies  at  1.78  Hz, 
requiring  an  extrapolation  down  to  the  local  oxygen  cyclotron  frequency.  At 
these  time  scales,  the  spectral  slope  appears  to  be  fairly  constant  over  the 
displayed  interval:  1  <  a  <2;  and  it  is  consistent  between  DC  and  LFC. 

In  comparing  panicles  (Fig.  12)  with  waves  (Fig.  13),  we  note  that  the 
v0  presentations  are  in  reasonable  agreement,  presenting  the  same  general 
trend  of  a  climax  for  the  event  near  UT  74800  s.  A  value  of  q  in  the  ran^c 
0.2  <  q  <  0.5  would  provide  reasonably  good  agreement  between  the  two 
sets  of  observations.  The  situation  for  o  is  perhaps  not  as  conclusive.  How¬ 
ever,  as  we  remarked  above,  its  determination  from  the  EICS  data  is  not 
likely  to  be  well  constrained,  particularly  for  o  >  2.  To  some  extent  the 
larger  values  of  a  near  UT  74800  s  may  be  an  aliasing  of  the  increases  in  v0 
at  that  time. 

To  briefly  summarize  the  results  obtained  for  this  event  [31],  there 
appears  to  be  a  rather  clear  correlation  between  the  va  parameter  determined 
separately  from  EICS  and  PWI.  The  situation  for  o  is  not  as  simple.  On  the 
one  hand,  the  agreement  appears  to  be  much  worse.  On  the  other  hand,  the 
disagreement  is  largely  due  to  the  variation  of  o  as  determined  from  the  ion 
data,  and  as  we  have  pointed  out  above,  there  are  a  number  of  reasons  not  to 
expect  a  to  be  accurately  computed  from  the  ion  data.  Thus,  given  the  agree¬ 
ment  for  va  we  should  be  inclined  to  forgive  the  a  variation  of  the  ions  in  an 
assessment  of  the  validity  of  the  theory  as  formulated. 

Of  course,  some  of  the  assumptions  which  we  have  made  in  order  to 
constrain  the  theory  and  simplify  the  comparison  of  the  two  instruments 
might  be  the  source  of  some  of  the  discrepancy.  For  example,  implicit  in  the 
assumption  that  lE^I2  -  q|£x|2  is  the  assumption  that  the  frequency  depen¬ 
dence  of  IFf.lf.is  the  same  power  law  as  \EX |2.  More  generally,  we  could  take 
|£L|2(/C1)  -  q|£J2(/cl )  but  assign  different  spectral  indices  aL  and  a*  to 
their  frequency  dependence.  Then  for  the  same  values  of  q,  it  is  easy  to 
reconcile  >  a,  with  the  observations,  and  it  is  rather  difficult  to  have 
aL  <  a*.  Since  the  o  parameter  relevant  to  the  conic  is  aL  +  1/3,  we 
would  then  expect  the  PWI  data  to  provide  a  lower  bound  for  the  value  of  a 
computed  from  the  EICS  data,  and  this  is  entirely  consistent  with  the  a  plots 
in  Figs.  12  and  13. 

Why  this  might  be  the  case  (i.e.,  aL  *  a*)  is  another  question.  In  the 
theory,  the  wave  power  at  higher  frequencies  is  used  to  provide  heating  at 
lower  altitudes.  Thus  in  a  model  for  the  event  where  the  wave  power  is  gen¬ 
erated  at  high  altitudes  and  is  propagating  downward,  there  would  be  some 
reflection  of  the  waves  which  would  tend  to  make  the  high  altitude  observa¬ 
tions  somewhat  more  energetic  than  those  at  lower  altitudes.  There  is  indeed 
some  suggestion  (Fig.  9  of  Ref.  [26])  that  the  higher  frequencies  are  less 
energetic  it  lower  altitudes  that  at  higher  ones.  Thus,  even  if  the  turbulence 
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was  unpolarized,  since  we  construct  a  model  of  the  turbulence  on  the  flux 
tube  based  on  the  high  altitude  observations,  we  must  have  aL  >  ax 

Another  significant  assumption  which  affects  the  analysis  is  the 
assumption  of  a  quasi-steady  evolution  process.  As  indicated  above,  this 
assumption  was  made  mostly  in  Older  to  simplify  the  analysis;  but  at  the 
same  time,  observations  by  a  single  satellite  do  not  provide  sufficient  infor¬ 
mation  to  include  these  effects.  Qualitatively,  we  still  expect  to  see  conics; 
quantitatively,  we  would  expect  some  corrections.  For  the  two  events  we 
have  presented,  the  wave  activity  could  be  characterized  by  a  temporal  turn¬ 
on,  turn-off  sequence.  Taking  causality  into  consideration,  we  would  then 
expect  to  see  the  peak  of  wave  activity  precede  that  of  the  conic  energy. 
Indeed,  an  allowance  of  -  20  s  on  day  81288  substantially  improves  the 
correlation  of  the  v0  plots  appearing  in  Figs.  12  and  13. 

V.  SUMMARY 

In  this  paper  we  have  discussed  the  ion  cyclotron  resonance  heating 
mechanism  as  it  applies  to  the  understanding  of  conic  events.  We  have 
approached  this  understanding  from  a  number  of  levels,  starting  with  the 
basic  physical  ideas  and  working  our  way  up  through  kinetic  theory  and  ulti¬ 
mately  on  to  a  comparison  with  the  observations.  Thus  we  have  seen  how  to 
estimate  the  zeroth  order  effects  with  a  simple  ion  heating  rate;  we  have  dis¬ 
cussed  how  to  formulate  the  problem  exactly  with  a  kinetic  equation  for  the 
distribution  /;  we  have  considered  a  number  of  analytic  and  numeric 
approaches  for  determining  / ;  we  have  seen  how  the  auroral  conics  are  best 
viewed  in  the  context  of  a  similarity  regime  determined  by  parameters  a  and 
v0;  and  we  have  seen  that  these  results  are  consistent  with  the  auroral  obser¬ 
vations  in  as  much  detail  as  the  observations  permit.  For  the  cusp  conics,  the 
work  of  comparison  is  still  in  progress,  but  the  results  to  date  [32-34]  are 
quite  promising. 

Of  course,  the  reader  is  reminded  that  one  should  not  conclude  from 
this  work  that  all  conics  are  ICRH  conics — we  have  not  considered  here  any 
evidence  which  detracts  from  the  viability  of  other  conic  generation  mechan¬ 
isms  given  the  proper  ionospheric/magnetospheric  setting.  Indeed,  many  of 
the  low-altitude  rocket  observations  of  energetic  conics  [56-59]  are  probably 
not  ICRH  conics  (or  at  least  not  as  we  have  discussed  them  here).  However, 
future  conics  discovered  in  the  presence  of  an  intense,  low-frequency,  elec¬ 
tromagnetic  turbulence  should  be  seriously  considered  as  ICRH  conics. 

There  are,  of  course,  many  questions  which  arise  in  the  course  of  this 
work.  These  include  the  application  of  ICRH  to  other  domains  in  space  plas¬ 
mas,  majority  species  heating,  the  ultimate  origins  of  the  low-frequency  tur¬ 
bulence,  possible  correlations  of  conic  generation  with  substorm  activity,  and 
other  questions  [31].  There  are  also  observational  uncertainties  due  to  the 
limitations  of  the  current  generation  of  instruments/satellites  which  hopefully 
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will  be  cured  in  future  missions.  (See  for  example  the  discussion  in  this 
volume  on  the  FAST  mission  [60] . )  As  is  the  rule,  one  can  expect  that  more 
questions  will  continue  to  be  raised  and  that  the  conic  story  is  far  from  over. 


REFERENCES 

[1]  R.D.  Sharp,  R.G.  Johnson,  and  E.G.  Shelley,  Observation  of  an  ionos¬ 
pheric  acceleration  mechanism  producing  energetic  (keV)  ions  primarily 
normal  to  the  geomagnetic  field  direction,  J.  Geophys.  Res.,  82,  pp. 
3324-3328  (1977). 

[2]  D.M.  Klumpar,  A  digest  and  comprehensive  bibliography  on  transverse 
auroral  ion  acceleration,  in  Ion  Acceleration  in  the  Magnetosphere  and 
Ionosphere ,  Geophysical  Monograph  Series,  Vol.  38,  Tom  Chang,  M.K. 
Hudson,  J.R.  Jasperse,  R.G.  Johnson,  P.M.  Kintner,  M.  Schulz  and  G.B. 
Crew,  eds.,  pp.  389-398  (American  Geophysical  Union,  Washington 
D.C.,  1986). 

[3]  J.L.  Burch,  Energetic  panicles  and  currents:  results  from  Dynamics 
Explorer,  Rev.  Geophys.,  26,  pp.  215-228  (1988). 

[4]  E.G.  Shelley,  W.K.  Peterson,  A.G.  Ghielmetti,  and  J.  Geiss,  The  polar 
ionosphere  as  a  source  of  energetic  magnetospheric  plasma,  Geophys. 
Res.  Lett.,  9,  pp.  941-944  (1982). 

[5]  A.W.  Yau,  B.A.  Whalen,  W.K.  Peterson,  and  E.G.  Shelley,  Distribution 
of  upflowing  ionospheric  ions  in  the  high-altitude  polar  cap  and  auroral 
ionosphere,  J.  Geophys.  Res.,  89,  pp.  5507-5522  (1984). 

[6]  E.G.  Shelley,  Magnetospheric  energetic  ions  of  terrestrial  origin  from 
the  Eanh’s  ionosphere,  Adv.  Space  Res.,  6,  pp.  121-132  (1987). 

[7]  C.R.  Chappell,  The  terrestrial  plasma  source:  a  new  perspective  in 
solar-terrestrial  processes  from  Dynamics  Explorer,  Rev.  Geophys.,  26, 
pp.  229-248  (1988). 

[8]  A.W.  Yau,  W.K.  Peterson,  and  E.G.  Shelley,  Quantitative  parametriza- 
tion  of  energetic  ionospheric  ion  outflow,  in  Modeling  Magnetospheric 


35. 


Plasma ,  Geophysical  Monograph  Series,  Vol.  44,  T.E.  Moore  and  J.H. 
Waite,  Jr.,  eds.,  pp.  211-217  (American  Geophysical  Union,  Washing¬ 
ton,  D.C.,  1988). 

[9]  R.  L.  Lysak,  in  Ion  Acceleration  in  the  Magnetosphere  and  Ionosphere , 
Geophysical  Monograph  Series,  Vol.  38,  T.  Chang,  M.  K.  Hudson,  J. 
R.  Jasperse,  R.  G.  Johnson,  P.  M.  Kintner,  M.  Schulz,  and  G.  B.  Crew, 
eds.,  pp.  261-270  (American  Geophysical  Union,  Washington,  D.C., 
1986). 

[10]  T.  Chang,  G.B.  Crew,  and  J.M.  Retterer,  Electromagnetic  tornadoes  in 
space:  ion  conics  along  auroral  field  lines  generated  by  lower  hybrid 
waves  and  electromagnetic  turbulence  in  the  ion  cyclotron  range  of  fre¬ 
quencies,  Computer  Physics  Communications,  49,  pp.  61-74  (1988). 

[11]  E.  Ungstrup,  D.M.  Klumpar,  and  W.J.  Heikkila,  Heating  of  ions  to 
superthermal  energies  in  the  topside  ionosphere  by  electrostatic  ion 
cyclotron  waves,  J.  Geophys.  Res.,  84,  pp.  4289-4296  (1979). 

[12]  R.L.  Lysak,  M.K.  Hudson,  and  M.  Temerin,  Ion  heating  by  strong  elec¬ 
trostatic  ion  cyclotron  turbulence,  J.  Geophys.  Res.,  85,  pp.  678-686 
(1980). 

[13]  K.  Papadopoulos,  J.D.  Gaffey,  and  P.J.  Palmadesso,  Stochastic  accelera¬ 
tion  of  large  M/Q  ions  by  hydrogen  cyclotron  waves  in  the  magneto¬ 
sphere,  Geophys.  Res.  Lett.,  7,  pp.  1014-1017  (1980). 

[14J  P.B.  Dusenbury  and  L.R.  Lyons,  Generation  of  ion  conic  distributions 
by  upgoing  ionospheric  electrons,  J.  Geophys.  Res.,  86,  pp.  7627-7638 
(1981). 

[15]  N.  Singh,  R.W.  Schunk,  and  J.J.  Sojka,  Energization  of  ionospheric 
ions  by  electrostatic  hydrogen  waves.  C?ophys.  Res.  Lett.,  8,  p.  1249 
(1981). 

[16]  M.  Ashour-Abdalla  and  H.  Okuda,  Turbulent  heating  of  heavy  ions  on 
auroral  field  lines,  J.  Geophys.  Res.,  89,  pp.  2235-2250  (1984),  and 
references  therein. 

[17]  T.  Chang  and  B.  Coppi,  Lower  hybrid  acceleration  and  ion  evolution  in 
the  suprauroral  region,  Geophys.  Res.  Lett.,  8,  pp.  1253-1256  (1981). 

[18]  G.B.  Crew  and  T.  Chang,  Asymptotic  theory  of  ion  conic  distributions, 
Phys.  Fluids,  28,  pp.  2382-2394  (1985). 

[19]  I.  Roth  and  M.K.  Hudson,  Lower  hybrid  heating  of  ionospheric  ions 
due  to  ion  ring  distributions  in  the  cusp,  J.  Geophys.  Res.,  90,  pp. 
4191-4203  (1985). 

[20]  J.M.  Retterer,  T.  Chang,  and  J.R.  Jasperse,  Ion  acceleration  by  lower 
hybrid  waves  in  the  suprauroral  region,  J.  Geophys.  Res.,  91,  pp. 
1609-1618  (1986). 

[21]  J.M.  Retterer,  T.  Chang,  and  J.R.  Jasperse,  Particle  acceleration  by 
intense  auroral  VLF  turbulence,  in  Physics  of  Space  Plasmas  (1989), 
SPI  Conference  Proceedings  and  Reprint  Series,  No.  9,  T.  Chang,  G.  B. 
Crew,  and  J.  R.  Jasperse,  eds.,  p.  119  (Scientific  Publishers,  Inc., 


36. 


Cambridge,  Mass.,  1990),  and  references  therein. 

[22]  W.  Lennartsson,  On  the  consequences  of  the  interaction  between  the 
auroral  plasma  and  the  geomagnetic  field.  Planet.  Space  Sci.,  28,  p.  135 
(1980). 

[23]  J.  E.  Borovsky,  The  production  of  ion  conics  by  oblique  double  layers, 
J.  Geophys.  Res.,  89,  p.  2251  (1984). 

[24]  M.  E.  Greenspan,  Effects  of  oblique  double  layers  on  upgoing  ion  pitch 
angle  and  gyrophase,  J.  Geophys.  Res.,  89,  p.  2842  (1984). 

[25]  J.D.  Winningham  and  J.  Burch,  Observation  of  large  scale  ion  conic 
generation  with  DE-1,  in  Physics  of  Space  Plasmas  (1982-4),  SPI 
Conference  Proceedings  and  Reprint  Series,  No.  5,  J.  Belcher,  H. 
Bridge,  T.  Chang,  B.  Coppi,  and  J.R.  Jasperse,  eds.,  pp.  137-158 
(Scientific  Publishers,  Inc.,  Cambridge,  MA,  1984). 

[26]  D.A.  Gumett,  R.L.  Huff,  J.D.  Menietti,  J.L.  Burch,  J.D.  Winningham, 
and  S.D.  Shawhan,  Correlated  low-frequency  electric  and  magnetic 
noise  along  the  auroral  field  lines,  J.  Geophys.  Res.,  89,  pp.  8971-8986 
(1984). 

[27]  T.  Chang,  G.B.  Crew,  N.  Hershkowitz,  J.R.  Jasperse,  J.M.  Retterer,  and 
J.D.  Winningham,  Transverse  acceleration  of  oxygen  ions  by  elec¬ 
tromagnetic  ion  cyclotron  resonance  with  broad  band  left-hand  polar¬ 
ized  waves,  Geophys.  Res.  Lett.,  13,  pp.  636-639  (1986). 

[28]  R.  Gendrin,  Wave  Particle  interactions  as  an  energy  transfer  mechanism 
between  different  particle  species.  Space  Sci.  Rev.,  34,  pp.  271-287 
(1983). 

l29]  J.M.  Retterer,  T.  Chang,  G.B.  Crew,  J.R.  Jasperse,  and  J.D.  Winning¬ 
ham,  Monte  Carlo  modeling  of  ionospheric  oxygen  acceleration  by 
cyclotron  resonance  with  broad-band  electromagnetic  turbulence,  Phys. 
Rev.  Lett.,  59,  pp.  148-151  (1987). 

[30]  D.M.  Klumpar,  W.K.  Peterson,  and  E.G.  Shelley,  Direct  evidence  for 
two- stage  (bimoda!)  acceleration  of  ionospheric  ions,  J.  Geophys.  Res., 
89,  pp.  10779-10787  (1984). 

[31]  G.  B.  Crew,  Tom  Chang,  J.  M.  Retterer,  W.  K.  Peterson,  D.  A.  Gumett, 
and  R.  L.  Huff,  Ion  cyclotron  resonance  heated  conics:  theory  and 
observations,  J.  Geophys.  Res..  95  (1990),  in  press. 

[32]  M.  Andr6,  H.  Koskinen,  L.  Matson,  and  R.  Erlandson,  Local  transverse 
ion  energization  in  and  near  the  polar  cusp,  Geophys.  Res.  Lett ,  15,  pp. 
107-110  (1988). 

[33]  W.K.  Peterson,  M.  Andrd,  G.B.  Crew,  A.M.  Persoon,  M.J.  Engebretson, 
C.J.  Pollock,  and  M.  Temerin,  Heating  of  thermal  ions  near  the  equator- 
ward  boundary  of  the  mid-altitude  polar  cleft,  in  Electromagnetic  Cou¬ 
pling  in  the  Polar  Clefts  and  Caps,  NATO  A  SI  Series  C:  Mathematical 
and  Physical  Sciences,  Vol.  278,  P.E.  Sandhold  and  A.  Egeland,  eds., 
pp.  103-113  (Kluwer  Academic  Publishers,  Dordrecht,  1989). 


37. 


[34]  M.  Andr£,  G.B.  Crew,  W.K.  Peterson,  A.M.  Persoon,  CJ.  Pollock,  and 
M.J.  Engebretson,  Heating  of  ion  conics  in  the  cusp/cleft,  in  Physics  of 
Space  Plasmas  (1989),  SPI  Conference  Proceedings  and  Reprint  Series, 
No.  9,  T.  Chang,  G.  B.  Crew,  and  J.  R.  Jasperse,  eds.,  p.  203  (Scientific 
Publishers,  Inc.,  Cambridge,  Mass.,  1990). 

[35]  Jay  R.  Johnson,  Tom  Chang,  G.  B.  Crew,  and  Mats  Andrd,  Equatorially 
generated  ULF  waves  as  a  source  for  the  turbulence  associate  J  with  ion 
conics,  Geophys.  Res.  Lett.,  16,  pp.  1469-1462  (1989). 

[36]  Jay  R.  Johnson,  Tom  Chang,  G.  B.  Crew,  and  Mats  Andrd,  Equatorially 
generated  ULF  waves  as  a  source  for  the  turbulence  associated  with 
ion  conics,  in  Physics  of  Space  Plasmas  (1989),  SPI  Conference 
Proceedings  and  Reprint  Series,  No.  9,  Tom  Chang,  G.  B.  Crew  and  J. 
R.  Jasperse,  eds.,  p.  433  (Scientific  Publishers,  Inc.,  Cambridge,  MA, 
1990). 

[37]  C.F.  Kennel  and  F.  Engelmann,  Velocity  space  diffusion  from  weak 
plasma  turbulence  in  a  magnetic  field,  Fhys.  Fluids,  9,  p.  2377  (1966). 

[38]  D.G.  Swanson,  Plasma  Waves  (Academic  Press,  Boston,  1989),  p.  296. 

[39]  R.Z.  Sagdeev  and  A.A.  Galeev,  Nonlinear  plasma  theory  (W.A.  Benja¬ 
min,  Inc.,  New  York,  1969),  pp.  54-55. 

[40]  G.B.  Crew  and  T.  Chang,  Asymptotic  theory  of  ion  conic  distributions, 
Phys.  Fluids,  28,  pp.  2382-2394  (1985). 

[41]  G.B.  Crew  and  T.  Chang,  Kinetic  treatment  of  oxygen  conic  formation 
in  the  central  plasma  sheet  by  broadband  waves,  in  Modeling  Magne- 
tospheric  Plasma,  Geophysical  Monograph  Series,  Vol.  44,  T.  Moore 
and  J.H.  Waite,  eds.,  pp.  159-163  (American  Geophysical  Union, 
Washington,  D.C.,  1988). 

[42]  G.  B.  Crew  and  T.  Chang,  Path  integral  formulation  of  ion  heating, 
Phys.  Fluids,  31,  pp.  3425-3439  (1988). 

[43]  G.B.  Crew,  M.  Andrd,  and  W.K.  Peterson,  Ion  cyclotron  resonance 
heating  in  the  cusp/cleft.  Bull.  Am.  Phys.  Soc.,  34,  p.  2009  (1989). 

[44]  M.  Andrd,  G.  B.  Crew,  W.  K.  Peterson,  A.  M.  Persoon,  C.  J.  Pollock, 
and  M.  J.  Engebretson,  Ion  heating  by  broadband  low-frequency  waves 
in  the  cusp/cleft,  in  preparation  for  J.  Geophys.  Res. 

[45]  J.M.  Retterer,  T.  Chang,  and  J.R.  Jasperse,  Ion  acceleration  in  the 
suprauroral  region:  a  Monte  Carlo  model,  Geophys.  Res.  Lett.,  10,  pp 
583-586  (1983). 

[46]  H.E.  Coffey,  Geomagnetic  and  solar  data,  J.  Geophys.  Res.,  87,  p.  926 
(1982). 

[47]  H.E.  Coffey,  Geomagnetic  and  solar  data,  J.  Geophys.  Res.,  87,  p.  1733 
(1982). 

[48]  H.E.  Coffey,  Geomagnetic  and  solar  data,  J.  Geophys.  Res.,  87,  p.  5310 
(1982). 

[49]  J.L.  Horwitz,  Velocity  filter  mechanism  for  ion  bowl  distributions 
(bimodal  conics),  J.  Geophys.  Res.,  91,  pp.  4513-4523  (1986). 


[50]  J.B.  Cladis,  Parallel  acceleration  and  transport  of  ions  from  polar  iono¬ 
sphere  to  plasma  sheet,  Geophys.  Res.  Lett.,  13,  pp.  893-896  (1986). 

[51]  M.  Temerin,  Evidence  for  a  large  bulk  ion  conic  heating  region  ,  Geo¬ 
phys.  Res.  Lett.,  13,  pp.  1059-1062  (1986). 

[52]  J.M.  Retterer,  T.  Chang,  G.B.  Crew,  J.R.  Jasperse,  and  J.D.  Winning- 
ham,  Monte  Carlo  modeling  of  oxygen  ion  conic  acceleration  by  cyclo¬ 
tron  resonance  with  broadband  electromagnetic  turbulence,  in  Physics  of 
Space  Plasmas  (1985-7),  SPI  Conference  Proceedings  and  Reprint 
Series,  No.  6,  Tom  Chang,  J.  Belcher,  J.  R.  Jasperse  and  G.  B.  Crew, 
eds.,  pp.  97-111  (Scientific  Publishers,  Inc.,  Cambridge,  MA,  1987). 

[53]  S.D.  Shawhan,  D.A.  Gumett,  D.L.  Odem,  R.A.  Helliwell,  and  C.G. 
Park,  The  plasma  wave  and  quasi-static  electric  field  instrument  (PWI) 
for  Dynamics  Explorer-A,  Space  Sci.  Instrum.,  5,  pp.  535-550  (1981). 

[54]  P.F.  Fougere,  On  the  accuracy  of  spectrum  analysis  of  red  noise 
processes  using  maximum  entropy  and  periodogram  methods:  simula¬ 
tion  studies  and  application  to  geophysical  data,  J.  Geophys.  Res.,  90, 
pp.  4355-4366  (1985). 

[55]  E.G.  Shelley,  D.A.  Simpson,  T.C.  Sanders,  E.  Hertzberg,  H.  Balsiger, 
and  A.  Ghielmetti,  The  energetic  ion  composition  spectrometer  (EICS) 
for  the  Dynamics  Explorer-A,  Space  Sci.  Instr.,  5,  pp.  443-454  (1981). 

[56]  B.A.  Whalen,  W.  Bernstein,  and  D.W.  Daly,  Low  altitude  acceleration 
of  ionospheric  ions,  Geophys.  Res.  Lett.,  5,  p.  55  (1978). 

[57]  A.W.  Yau,  B.A.  Whalen,  A.G.  McNamara,  P.J.  Kellogg,  and  W.  Bern¬ 
stein,  Particle  and  wave  observations  of  low-altitude  ionospheric  ion 
acceleration  events,  J.  Geophys.  Res.,  88,  pp.  341-355  (1983). 

[58]  P.M.  Kintner,  J.  LaBelle,  W.  Scales,  A.W.  Yau,  and  B.A.  Whalen, 
Observations  of  plasma  waves  within  regions  of  perpendicular  ion 
acceleration,  Geophys.  Res.  Lett.,  13,  pp.  1113-1116  (1986). 

[59]  P.M.  Kintner,  W.  Scales,  J.  Vago,  R.  Amoldy,  G.  Garbe,  and  T.  Moore, 
Simultaneous  observations  of  electrostatic  oxygen  cyclotron  waves  and 
ion  conics,  Geophys.  Res.  Lett.,  16,  pp.  739-742  (1989). 

[60]  M.A.  Temerin,  C.W.  Carlson,  C.A.  Cattell,  R.E.  Ergun,  J.P.  McFadden, 
F.S.  Mozer,  D.M.  Klumpar,  W.K.  Peterson,  E.G.  Shelley,  and  R.C. 
Elphic,  Wave-particle  interactions  on  the  FAST  satellite,  in  Physics  of 
Space  Plasmas  (1989),  SPI  Conference  Proceedings  and  Reprint  Series, 
No.  9,  T.  Chang,  G.  B.  Crew,  and  J.  R.  Jasperse,  eds.,  p.  343  (Scientific 
Publishers,  Inc.,  Cambridge,  Mass.,  1990). 


39. 


3.  PARTICLE  ACCELERATION  BY 
INTENSE  AURORAL  VLF  TURBULENCE 

John  M.  Rettercr1 

Air  Force  Geophysics  Laboratory,  Bedford,  VIA  01731 
Tom  Chang 

Center  for  Space  Research,  MIT,  Cambridge,  MA  02139 
J.  R.  Jasperse1 

Air  Force  Geophysics  Laboratory,  Bedford,  MA  01731 
ABSTRACT 

Broadband  turbulence  in  the  lower-hybrid  to  plasma  frequency  range  is 
found  in  a  variety  of  forms  in  the  suprauroral  region,  most  notably  as  auroral 
hiss  and  VLF  saucers.  When  the  turbulence  is  intense,  it  is  observed  to  be 
associated  with  ion  conics  (ions  heated  transverse  to  the  geomagnetic  field) 
and  “counter-streaming”  electron  fluxes  (heated  in  both  directions  parallel  to 
the  field).  This  tutorial  will  begin  with  a  review  of  the  dispersion  and  propa¬ 
gation  characteristics  of  whistler  resonance-cone  waves,  which  comprise  the 
turbulence,  and  go  on  to  discuss  the  theories  for  the  excitation  of  the  tur¬ 
bulence.  Plasma  simulation  and  mesoscale  (Monte  Carlo)  simulation  tech¬ 
niques  will  be  used  to  illustrate  the  interaction  of  the  ambient  plasma  with  the 
turbulence.  These  calculations  will  demonstrate  how  this  interaction  results  in 
transverse  heating  of  the  ions  and  parallel  heating  of  the  electrons  of  the 
plasma,  leading  to  the  formation  of  the  observed  heated  and  accelerated  parti¬ 
cle  fluxes. 


I.  INTRODUCTION 

Auroral  hiss  is  a  plasma-wave  phenomenon  commonly  observed  in  the 
acceleration  region  above  the  Earth’s  auroral  zone,  associated  with  inverted-V 
electron  precipitation  [1].  The  frequencies  of  these  waves,  extending  from 
the  lower-hybrid  resonance  up  to  the  plasma  frequency,  fall  in  the  VLF  range. 
Their  dispersion  characteristics  have  led  them  to  be  identified  with  waves  on 
the  whistler  resonance-cone  plasma  dispersion  surface,  which  are  excited  by  a 
plasma  microinstability  caused  by  the  accelerated  auroral  electrons.  The  elec¬ 
tric  field  intensities  of  these  waves  are  occasionally  large  12],  and  this  tur¬ 
bulence  has  been  observed  to  be  correlated  [3]  with  the  energy  flux  of 
transversely  accelerated  ions  observed  in  the  topside  ionosphere  (see  Ref.  [4] 
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for  a  bibliography). 

Above  the  ionosphere,  these  accelerated  ions  are  observed  in  a  form 
called  an  ion  conic,  because  their  pitch-angle  distribution  peaks  at  an  oblique 
angle.  A  mechanism  by  which  wave-panicle  interaction  of  ions  with  the 
intense  turbulence  near  the  lower-hybrid  frequency  leads  to  ion  conic  forma¬ 
tion  was  suggested  by  Chang  and  Coppi  [5]:  first,  the  turbulence  is  excited  in 
the  ambient  plasma  of  the  suprauroral  region  by  the  accelerated  auroral  elec¬ 
trons.  Wave-panicle  interaction  of  the  ambient  ion  population  with  the  tur¬ 
bulence  near  the  lower-hybrid  frequency  then  leads  to  ion  acceleration  nearly 
perpendicular  to  the  field  line,  which  is  followed  by  the  adiabatic  folding  of 
velocities  as  the  ions  mirror  and  travel  up  the  field  line,  to  create  the  conic 
velocity  distribution. 

Plasma  simulations  [6]  have  shown  this  mechanism  to  be  effective  not 
only  for  the  transverse  acceleration  of  ions,  but  also  for  the  acceleration  of 
the  ambient  electrons  in  the  directions  parallel  to  the  geomagnetic  field.  The 
electron  velocity  distribution  is  found  to  have  enhanced  fluxes  of  energetic 
electrons  in  both  directions  with  respect  to  the  geomagnetic  field,  a  form 
which  is  indeed  observed  in  space  in  conjunction  with  ion  conics:  counter¬ 
streaming  electrons  [7], 

The  wave-panicle  interaction  processes  involved  in  the  generation  of 
the  VLF  turbulence  and  the  acceleration  of  the  ambient  plasma  are  found  to 
occur  on  spatial  and  temporal  scales  large  compared  to  those  which  can  be 
practically  studied  using  self-consistent  panicle  plasma  simulation.  To  model 
ion  conics,  a  larger-scale  simulation  model  has  been  employed  [8].  This 
model  of  global  wave-panicle  interaction  employs  a  Monte  Carlo  technique  to 
follow  the  trajectories  of  many  ions  as  they  undergo  the  influences  of  wave- 
panicle  interaction  and  large-scale  electric  and  magnetic  fields  in  the  suprau¬ 
roral  region.  It  has  enjoyed  considerable  success  in  modeling  the  ion  velocity 
distribution  of  ion  conics  formed  in  other  regions  by  interaction  with  tur¬ 
bulence  in  the  ion  cyclotron  frequency  range  [9;  10).  (See  the  review  by 
Crew  and  Chang  (11],  for  a  discussion  of  this  other  form  of  ion  conic.) 
Because  the  wave-panicle  interaction  in  this  other  case  involves  a  simple  fre¬ 
quency  resonance,  an  empirical  estimate  of  the  rate  of  panicle  heating  using 
the  observed  frequency  spectrum  of  the  turbulence  directly  is  found  to  be 
satisfactory.  In  the  case  of  interaction  with  VLF  turbulence,  on  the  other 
hand,  the  resonance  requires  the  matching  of  velocities,  and  considerable 
uncertainty  exists  in  using  observed  wave  amplitudes  to  estimate  the  wave- 
panicle  interaction  rate  because  the  lack  of  wave-vector  measurements 
prevents  us  from  knowing  the  phase  velocities  of  the  waves.  This  problem 
has  been  addressed  using  models  for  the  micro-physics  of  the  local  wave- 
particle  interaction,  based  on  a  combination  of  observation  and  simulation 
results. 

In  t  :s  review,  we  will  first  touch  on  the  wave  dispersion  and  propaga¬ 
tion  vharacteristics  and  discuss  the  excitation  mechanism  of  auroral  hiss.  The 
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local  model  for  wave-particle  interactions  will  be  described  and  illustrated 
using  simulation  results  and  a  simple  nonlinear  theory,  which  explains  how 
the  counter-streaming  electrons  and  hot  ion  tails  are  produced.  Finally,  the 
application  of  the  local  results  to  the  meso-scale  model  for  ion  conic  forma¬ 
tion  in  the  suprauroral  region  will  be  discussed,  using  two  examples  to  illus¬ 
trate  the  method  and  the  results  achieved  by  the  theory. 

II.  AURORAL  HISS  AND  LOWER-HYBRID  WAVES 

The  plasma  waves  that  are  observed  in  conjunction  with  ion  conics  and 
counter-streaming  electrons  are  an  intense  form  of  an  ionospheric  plasma 
wave  phenomenon  known  as  auroral  hiss.  We  begin  here  with  a  brief  review 
of  the  literature  of  observations  of  auroral  hiss  and  our  theoretical  understand¬ 
ing  of  the  phenomenon. 

From  their  first  flights,  polar-orbitting  satellites  carrying  VLF  radio 
receivers  have  measured  a  broad-band  emission  at  high  latitudes,  ranging  in 
frequency  from  the  lower-hybrid-resonance  frequency  up  through  the  whistler 
band  [12;  13].  Correlation  with  panicle  spectrometer  data  [14]  demonstrated 
the  association  between  auroral-zone  VLF  hiss  and  intense  fluxes  of  precipi¬ 
tating  electrons  in  inverted-V  events  and  suggested  that  the  generation  of  the 
VLF  modes  occured  through  Landau  resonance  with  the  accelerated  auroral 
electrons.  Evaluation  of  the  Poynting  flux  of  these  waves  showed  that  the 
auroral  hiss  was  indeed  propagating  towards  the  Eanh  along  with  the  auroral 
electrons,  and  distinguished  auroral  hiss  from  a  similar  wave  phenomenon, 
known  as  VLF  saucers  [15],  whose  Poynting  flux  at  the  same  altitude  was 
observed  to  flow  away  from  the  Earth.  VLF  saucers  are  located  equatorward 
of  the  auroral  hiss,  and  were  not  found  to  be  correlated  with  energetic  elec¬ 
tron  precipitation,  implying  that  their  source  of  free  energy  resided  in  some 
unobserved  current  carrier,  possibly  cold  ionospheric  electrons  in  the  auroral 
return  current  region.  VLF  emissions  have  also  been  observed  in  conjunction 
with  -  50  eV  electron  beams  flowing  away  from  the  Earth  in  the  cusp  region 
[16].  The  characteristic  funnel-shaped  dispersion  of  auroral  hiss  on  a 
frequency-latitude  spectrogram — that  is,  higher-frequency  waves  being  more 
spread  out  latitudinally  than  lower-frequency  waves — can  be  explained  in 
terms  of  the  propagation  characteristics  of  whistler  waves  propagating  near 
the  resonance  cone  [I;  15],  which  we  will  briefly  describe  below. 

A.  Dispersion  Theory 

To  understand  the  plasma  waves  that  make  up  auroral  hiss,  we  can 
begin  with  die  linear  theory  of  a  stable,  cold,  homogeneous  plasma,  and  then 
make  the  necessary  generalizations.  In  a  cold  magnetized  plasma,  the  linear 
dispersion  of  electromagnetic  waves  is  described  by  the  following  dispersion 
relation,  using  the  notation  of  Stix  [17]: 
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where  n  -ic-co  .s  the  index  of  refraction,  6  is  the  angle  of  the  wave  vector  k 
with  respect  to  the  background  magnetic  field  B ,  and 
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where  copj  is  the  plasma  frequency  and  Qf  the  gyro  frequency  of  species  s . 
We  isolate  the  relevant  root  of  this  equation  for  whistler  resonance-cone 
waves  by  considering  the  resonance  limit,  in  which  n  approaches  infinity.  In 
this  case,  tan2  0  =  -P  /S ,  and  with  the  ordering  appropriate  for  VLF  waves, 
£2?«co2«Q2,  we  find  the  root 


cu2  =  co^  s 


CQp2,  +  co2,cos29 
1  4-  co2  sin29/Q2 


(3) 


where  we  have  assumed  that  the  plasma  consists  of  electrons  and  one  ion 
species,  with  plasma  frequencies  respectively  <ape  and  cup( .  We  see  that  the 
mode  has  frequencies  ranging  from  the  lower  hybrid  resonance  frequency 
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when  k  is  perpendicular  to  B,  up  to  the  plasma  frequency  o)p  (where 
tDp2  =  o)/*+cop2 )  when  k  is  parallel  to  B .  (If  we  analysed  the  dispersion  relation 
slightly  more  generally,  we  would  find  that  the  upper  limit  is  the  smaller  of 
the  plasma  or  electron  gyro  frequencies,  but  the  ordering  we  have  selected  is 
the  appropriate  one  for  the  auroral  acceleration  region,  at  altitudes  of  the 
order  of  1  Rp>t  where  cop</Qe  is  small.) 

These  waves  propagate  near  the  whistler  resonance  cone,  where  to  first 
approximation  the  waves  are  electrostatic  in  nature  and  the  frequency  of  a 
mode  depends  only  on  the  angle  of  propagation  with  respect  to  the  magnetic 
field.  When  to  is  independent  of  It  I,  the  group  velocity,  which  defines  the 
direction  and  speed  with  which  wave  energy  propagates,  is  perpendicular  to 
k.  Waves  with  frequencies  near  co^*,  whose  wave  vectors  are  approxi¬ 
mately  perpendicular  to  B ,  have  group  velocities  approximately  parallel  to  B , 
while  waves  near  the  plasma  frequency  cop,  have  group  velocities  approxi¬ 
mately  perpendicular  to  B.  Thus,  in  the  auroral  region  where  B  is  nearly 
vertical,  higher  frequency  waves  on  the  whistler  resonance  cone  propagate 
further  in.  latitude  than  lower-frequency  waves  from  the  same  source. 


explaining  the  funnel  or  saucer  shape  dispersion  characteristics  of  VLF  hiss 
observed  by  satellite.  By  tracing  the  paths  of  the  low  and  high  frequency 
portions  of  the  observed  spectrum  back  to  where  they  intersect,  the  distance 
from  the  source  of  the  hiss  or  saucer  can  be  inferred  Li;  15],  although  the 
direction  along  B ,  whether  above  or  below  the  spacecraft,  cannot  be  deter¬ 
mined  from  the  dispersion  alone.  The  filled- in  funnel  spectrum  observed  for 
auroral  hiss  presumably  implies  a  source  distributed  along  the  geomagnetic 
field  line,  while  the  saucer  spectrum,  which  is  only  outlined  by  the  dispersion 
curve,  implies  a  source  which  is  confined  in  extent  along  B  [15],  although  it 
can  be  extended  in  longitude  [18]. 

The  first  generalization  we  must  consider  in  preparation  for  our  study 
of  wave-particle  interaction  with  these  waves  is  the  introduction  of  kinetic 
effects  and  the  consequences  of  finite-temperature  plasma  particle  species. 
The  dielectric  tensor  for  electromagnetic  waves  in  a  Maxwellian  plasma  can 
be  found  in  many  texts  [19]  and  we  will  not  reproduce  it  here.  Instead,  we 
can  begin  with  the  facts  that  the  waves  of  interest  are  essentially  electrostatic 
and  within  the  frequency  range  Qi«ax<Qe  to  write  a  simpler  dispersion 
relation.  Because  wave  frequencies  are  larger  than  the  ion  gyrofrequency,  to 
the  first  approximation  we  may  use  the  unmagnetized  susceptibility  [19,  p.56] 
for  the  ions.  For  the  electrons,  we  can  use  the  magnetized  susceptibility  [19, 
p.73],  retaining  from  the  sum  over  cyclotron  resonances  only  the  zeroth  order 
term,  making  the  strongly  magnetized  approximation  for  the  electrons  because 
the  frequencies  of  interest  are  much  smaller  than  the  electron  gyrofequency. 
The  dispersion  relation  for  electrostatic  waves  in  a  Maxwellian  plasma  is  then 


e(*,o»  =  l  +  £x,(*,ci>)  =  0 
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where  the  sum  of  susceptibilities  is  over  ail  species  in  the  plasma,  or 
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for  the  simple  ion-electron  plasma.  In  this  expression,  W  is  the  plasma 
dispersion  function  as  presented  by  Ichimaru  [19,  p.56],  is  the  component 
of  k  along  B ,  and  v„  is  the  thermal  velocity,  yjTslms.  By  introducing  kinetic 
effects,  we  allow  the  possibility  of  damping  or  excitation  of  the  waves  by 
resonant  wave-particle  interaction.  If  we  consider  both  kvti  and  to  be 
small  compared  to  the  frequencies  of  interest,  we  can  solve  the  dispersion 
relation,  Eq.(6),  using  a  perturbative  approach.  To  zeroth  order  in  these 
terms,  we  recover  the  solution,  Eq.(3),  except  with  the  slight  error  that  sin2© 
is  missing  from  the  second  term  in  the  denominator.  To  first  order  in  (kvti  )2 
and  (Jk»vM)2,  co2  contains  positive  additional  terms  proportional  to  these  fac¬ 
tors,  much  like  the  thermal  correction  of  the  Langmuir  mode  dispersion.  The 
imaginary  part  of  the  frequency  can  be  calculated  using  the  following 
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approximation,  valid  for  small  imaginary  frequency 

-  -Ime 
^  3Ree/8o) 

where  e  is  evaluated  using  the  real  frequency.  We  will  not  write  out  this 
result  until  the  next  section,  where  we  include  the  auroral  electron  beam 
population  and  the  opportunity  for  positive  growth  rates  and  instability. 

The  introduction  of  the  effects  of  magnetized  ions  changes  the  disper¬ 
sion  at  short  wavelengths  for  nearly  perpendicular  propagation.  When  the 
perpendicular  wavelength  becomes  comparable  to  the  ion  gyroradius. 
Pi  =  v(1  /Q, ,  the  frequency  plateau  near  the  lower-hybrid  resonance  frequency 
which  existed  at  longer  wavelength  begins  to  break  up  into  the  high-harmonic 
ion  Bernstein  modes  around  co^ ,  separated  approximately  by  the  ion  gyroffe- 
quency.  The  upper  envelope  of  the  group  of  Bernstein  modes  is  still  roughly 
described  by  the  unmagnetized  lower-hybrid  frequency  including  the  thermal 
correction.  This  phenomenon  is  illustrated  well  by  the  plasma  dispersion  sur¬ 
faces  presented  by  Andrd  [20],  calculated  numerically  using  the  dispersion 
code  WHAMP  [21].  The  damping  effect  of  magnetized  ions  was  demon¬ 
strated  in  calculations  by  Gomey  et  a/. [22]  who  showed  that  damping  of 
auroral  hiss  by  energetic  ions,  i.e.,  ion  conics  themselves,  can  produce  observ¬ 
able  absorption  lines  in  the  hiss  spectrum  at  multiples  of  the  ion  gyrofre- 
quency. 


B.  Resonance-Cone  Plasma  Wave  Instability 

With  the  association  of  auroral  hiss  with  energetic  electron  precipitation 
in  inverted-V  events  came  the  recognition  [14]  that  a  kinetic  instability  of  the 
plasma  containing  both  the  accelerated  auroral  electrons  and  a  ambient  popu¬ 
lation  of  electron  plus  ion  species  could  excite  the  waves  through  the  Landau 
resonance  with  the  auroral  electrons,  where  ubt  is  the  velocity  of  the 

auroral  electrons  along  the  geomagnetic  field.  We  first  consider  simple 
modeling  using  Maxwellians  for  all  plasma  species,  and  then  discuss  more 
detailed  modeling  that  can  be  performed. 

The  effect  of  the  DC  electric  field  within  the  inverted-V  potential  struc¬ 
ture  in  the  direction  of  the  geomagnetic  field  is  to  accelerate  the  plasma  sheet 
electron  population  of  the  magnetosphere  along  the  magnetic  field  lines 
toward  the  Earth  [2;23].  The  resulting  plasma  below  the  potential  drop  is  a 
mixture  of  the  hot,  drifting  auroral  electron  component  and  the  much  cooler 
ambient  plasma  of  ionospheric  origin.  (In  addition,  electron  components  of 
intermediate  energies,  due  to  atmospheric  backscatter  and  energy  degradation 
of  the  auroral  electrons,  wave-particle  interaction  processes,  etc.,  are  also 
present;  we  will  mention  their  effects  later.)  This  situation  gives  rise  to  a 
beam-pl  ma  instability  of  the  bump-on-tail  form  which  can  excite  the 
v  rstler  resonance-cone  waves  of  auroral  hiss. 
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The  simple,  electrostatic  linear  theory  will  be  our  first  introduction  to 
this  instability  [5;24].  In  addition  to  the  ambient  electron  and  ion  species, 
modeled  as  Maxweilians  of  equal  temperature  (7,  -  Tt~  1  eV),  the  unstable 
plasma  contains  the  energetic  electrons  accelerated  through  the  field-aligned 
potential  drop.  This  species  is  modeled  as  a  warm  Maxwellian  distribution, 
drifting  parallel  to  the  magnetic  field  with  velocity  ubl  such  that  the  drift 
kinetic  energy  is  >  1  keV  and  its  temperature  is  a  few  hundred  eV.  The  den¬ 
sity  of  the  beam  population  is  small — typically,  less  than  10~2  times  the 
ambient  density.  To  incorporate  the  beam  population  into  the  dielectric  func¬ 
tion  of  the  electron-ion  plasma,  we  add  its  susceptibility 


to  Eq.(6).  If  the  beam  density  is  small  enough,  its  effect  on  the  real  fre¬ 
quency  of  the  resonance-cone  mode  is  small.  When  we  examine  the  ima¬ 
ginary  component  of  the  frequency,  describing  the  growth  or  damping  of  the 
wave,  we  obtain 
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where  to  is  the  real  frequency  of  the  wave,  Eq.(3).  We  see  that  y,  the  ima¬ 
ginary  part  of  the  frequency,  can  be  positive,  implying  exponential  growth  of 
the  waves,  when  the  first  term,  due  to  the  auroral  beam  population,  is  positive 
and  is  larger  than  the  damping  contributions  from  the  ambient  plasma  elec¬ 
trons  and  ions,  described  by  the  second  and  third  terms.  Typically,  this  can 
happen  for  phase  velocities  ubt-r\vtb  <u>ikt<ubv  where  q  is  a  small  factor. 

The  plasma  of  the  auroral  ionosphere,  however,  is  not  a  cold  Maxwel¬ 
lian,  and  the  presence  of  a  warm  secondary  population  of  elections  can 
strongly  affect  the  instability  by  acting  to  damp  the  waves  [25].  The  study  of 
the  stability  of  the  auroral  plasma  must  be  done  using  an  empirical 
specification  of  the  electron  velocity  distribution.  A  useful  technique 
employed  by  Lotko  and  Maggs  [25]  is  to  fit  an  observed  distribution  with  a 
sum  of  Maxweilians  of  different  densities,  temperatures,  and  drifts;  this 
allows  us  to  use  the  well-understood  and  accurately  tabulated  Maxwellian 
plasma  dispersion  function  to  analyze  a  more  general  plasma.  An  additional 
refinement  is  to  calculate  the  solution  of  the  complete  electromagnetic  disper¬ 
sion  relation.  This  must  necessarily  be  done  numerically,  for  which  we  use 
the  WHAMP  code  of  K.  Ronnmark,  kindly  provided  to  us  by  M.  Andrf.  Fig. 


Figure  l.  Linear  calculation  of  resonance-cone  frequency  and 
growth  rate.  The  top  panel  presents  the  real  part  of  the  frequency 
as  a  function  of  k i  and  *1  using  contours,  while  the  lower  panel 
prt--:nts  'he  imaginary  part  or  growth  rate. 


1  presents  the  results  of  the  instability  calculation  at  an  altitude  of  a  few 
thousand  kilometers,  using  a  representative  combination  of  plasma  popula¬ 
tions:  a  cool  ionospheric  electron  component  of  density  130  cm-3  and  a  tem¬ 
perature  of  1  eV;  a  combination  of  three  warm  electron  components  to 
represent  an  approximate  power  law  flux,  between  10  eV  and  1000 

eV,  with  a  total  density  of  about  1.6;  the  auroral  electron  component  of  den¬ 
sity  0.05,  temperature  100  eV,  drifting  with  an  energy  of  4000  eV;  and  finally 
a  cool  proton  component  of  temperature  1  eV.  The  top  panel  gives  the  real 
pan  of  the  frequency  plotted  using  contours  of  constant  frequency  as  a  func¬ 
tion  of  fcnPi  and  k1pt .  The  contours  are  labeled  with  the  corresponding  fre¬ 
quencies,  expressed  as  a  multiple  of  the  ion  gyrofrequency.  The  straight-line 
nature  of  the  contours  in  the  center  of  the  plot  shows  the  essential  depen¬ 
dence  of  the  frequency  on  the  angle  of  propagation,  as  discussed  in  the  previ¬ 
ous  subsection.  Near  parallel  propagation  we  see  the  plasma- frequency  pla¬ 
teau  (cop=180Q,),  with  a  slight  frequency  perturbation  due  to  the  instability. 
Near  perpendicular  propagation,  we  find  the  lower-hybrid  plateau 
^lhr  =4.2Q, ),  modified  by  the  thermal  correction  near  the  border  of  the 
figure.  In  the  lower  panel,  the  contours  of  positive  growth  rate  as  a  fraction 
of  the  ion  gyrofrequency  are  spaced  logarithmically,  two  per  decade  from 
0.001  for  the  outermost  contour  up  to  2.0  at  the  peak.  The  growth  rate  con¬ 
tours  approximately  follow  lines  of  constant  c o/kw,  due  to  the  Landau  reso¬ 
nance  that  excites  the  instability.  (We  see  this  by  noting  the  horizontal  con¬ 
tours  in  the  two  frequency  plateau  regions,  and  the  vertical  contours  in  the 
intermediate  region.)  From  the  larger  growth  rates  at  higher  frequency,  one 
might  infer  that  the  spectrum  of  auroral  hiss  should  be  dominated  by  emis¬ 
sions  near  the  plasma  frequency.  That  this  is  not  the  case  is  due  to  the  fact 
that  the  instability  is  a  convective  one,  which  in  the  auroral-arc  geometry  has 
great  consequences  for  the  form  of  the  auroral  hiss  spectrum. 

C.  Convective  Saturation  and  the  Auroral  Hiss  Spectrum 

The  geometry  of  an  auroral  arc  implies  that  the  source  of  free  energy 
for  the  excitation  of  auroral  hiss  is  narrowly  confined  in  latitude,  but  can  be 
extended  along  the  geomagnetic  field  line,  following  the  trajectories  of  the 
accelerated  auroral  electrons.  In  addition  to  the  structure  of  the  electron 
beam,  the  background  plasma  in  the  vicinity  of  auroral  arcs  is  often  struc¬ 
tured,  too;  characteristically,  a  field-aligned  density  depletion  or  cavity  is 
observed,  oriented  along  the  auroral  field  lines  (26).  This  structure  affects  the 
propagation  of  the  resonance-cone  waves  excited  within  the  arc,  and  can 
affect  the  resonant  absorption  of  the  waves  by  the  ambient  plasma. 

The  consequences  of  the  geometry  of  the  precipitation  region  and  the 
propagation  of  the  auroral-hiss  waves  through  the  inhomogeneous  ionosphere 
were  pointed  out  by  Maggs  [27;  28).  Waves  near  the  electron  plasma  fre¬ 
quency,  with  wave  vectors  nearly  aligned  with  the  geomagnetic  field,  have 
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group  velocities  nearly  perpendicular  to  the  field  line  and  so  can  propagate 
quickly  out  of  the  auroral  arc,  limiting  their  amplitude.  The  group  velocities 
of  waves  near  the  lower  hybrid  frequency,  on  the  other  hand,  are  directed 
nearly  along  the  magnetic  field  line,  allowing  these  waves  to  stay  in  reso¬ 
nance  longer  and  possibly  grow  to  amplitudes  at  which  nonlinear  effects  are 
important  for  saturation.  Using  a  ray-tracing  calculation,  Maggs  showed  that 
many  of  the  features  of  the  auroral  hiss  spectrum  could  be  explained  by  the 
convective  beam  amplification  of  incoherent  whistler  radiation  as  it  pro¬ 
pagates  through  the  inhomogeneous  ionosphere. 

The  basic  assumption  of  the  ray -trajectory  method  is  that  the 
wavelength  must  be  short  compared  to  the  distance  over  which  the  refractive 
index  changes  appreciably  [17],  In  a  time-independent  medium  where  the 
geometrical-optic  assumption  above  applies,  the  path  along  which  the  energy 
of  a  stable  wave  propagates  is  described  by  the  set  of  equations 

dx  _  deJdk  _  3co  _  dk_  Be/Bx 

dt  BeJB  a)  Bk  g  dt  BeJBw 

[17],  Here  the  first  equation  describes  the  path  through  space-time  whose 
tangent  is  the  group  velocity  of  the  wave,  while  the  second  equation  describes 
the  change  in  the  wave  vector  which  is  necessary  so  that  the  dispersion  rela¬ 
tion  for  a  locally  homogeneous  plasma,  e(* ,  to,  x )  =  0,  is  satisfied  everywhere 
along  the  ray  path  for  a  wave  whose  frequency  stays  constant. 

Let  us  consider  the  propagation  of  a  resonance-cone  wave  excited  at 
high  altitude  propagating  toward  the  Earth  in  the  auroral  density  cavity,  with 
density  gradients  along  the  field  line  towards  the  Earth  and  transverse  (in  the 
latitudinal  direction)  in  both  directions  away  from  the  density  minimum  at  the 
center  of  the  cavity.  As  the  wave  propagates  toward  the  Earth,  the  parallel 
component  of  its  wave  vector  falls  to  keep  the  wave  on  the  resonance -cone  of 
the  cold-plasma  dispersion  relation,  until  the  local  lower-hybrid  resonance  fre¬ 
quency  has  risen  to  meet  the  wave’s  frequency  and  k N  falls  to  zero.  At  this 
point  the  wave  will  be  reflected,  although  the  description  of  this  process 
requires  a  full-wave  treatment,  because  the  assumptions  of  ray-tracing  will 
have  failed  before  this  point,  when  the  wave’s  wavelength  became  so  long 
that  the  geometrical  optics  approximation  was  violated.  Propagation  in  the 
transverse  direction  is  marked  by  a  tendency  for  the  resonance-cone  waves  to 
be  refracted  towards  high  density,  or  out  of  the  depleted  density  cavity;  as  the 
waves  propagate  insversely  into  increasing  density,  k{  increases  until  the 
resonance  is  reac;  d  where  co  equals  At  this  point  the  wave  will  be 

absorbed.  Again,  because  refractive  characteristics  of  the  medium  are  chang¬ 
ing  rapidly  at  this  point,  a  full-wave  treatment  of  propagation  is  needed  to 
describe  the  situation  accurately.  The  phenomena  of  absorption  at  the  lower- 
hybrid  resonance  remain  poorly  understood;  in  addition  to  the  need  for  a 
kinetic  d  'crintion,  the  multiple  Bemstein-mode  branches  nearby  in  the 
disocrciuu  relation  make  the  description  of  possible  linear  mode  conversion 
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processes  a  nightmare  [29], 

The  transverse  structure  of  the  inhomogeneous  electron  beam  in  the 
auroral  arc  can  cause  the  dielectric  function  to  change  rapidly  in  space,  too, 
modifying  the  growth  of  the  resonance-cone  waves.  One  way  is  through  the 
“de-tuning"  of  the  Landau  resonance  caused  by  the  shear  in  the  magnetic- 
field  geometry  produced  by  the  current  generated  by  the  auroral  electron 
beam.  A  WKB  analysis  of  the  wave  equation  and  us  turning  points  [ 30 1 , 
however,  shows  for  the  appropriate  scales  and  parameters  of  the  auroral  arc 
situation  that  the  excitation  of  resonance-cone  waves  there  is  little  affected  by 
non-local  effects. 

The  waves  excited  by  the  auroral  beam  are  subject  to  the  influences  of 
both  the  parallel  and  perpendicular  density  gradients,  to  differing  degrees 
depending  on  their  propagation  direcdons  and  thus  their  frequencies.  Fig.  2 
illustrates  four  ray  trajectories  for  waves  of  different  frequencies.  These  four 
rays  start  with  kL  oriented  along  the  arc,  (in,  say,  the  East- West  direcdon),  but 
due  to  the  transverse  density  gradient  acquire  a  component  of  kL  perpendicu¬ 
lar  to  the  arc  (the  North-South  direction),  and  begin  propagadng  out  of  it. 
The  highest  frequency  wave  propagates  transversely  out  of  the  auroral  arc 
soonest  because  its  group  velocity  is  directed  more  nearly  perpendicular  to  B . 
The  ends  of  the  trajectories  are  depicted  with  dashed  lines  to  remind  us  that 
the  geometrical-optic  assumption  of  the  calculation  begins  to  fail  in  that  por¬ 
tion  of  the  trajectory.  As  a  wave  propagates  along  its  path,  its  phase  velocity 
changes,  so  it  can  fall  into  or  out  of  resonance  with  the  auroral  beam,  as  well 
as  be  convected  out  of  the  auroral  arc.  When  the  amplitude  of  the  auroral 
hiss  is  weak  and  does  not  perturb  the  plasma,  Maggs  [27]  showed  how  to  use 
ray-tracing  calculations  to  predict  the  spectral  density  of  the  auroral  hiss  at  a 
point  in  space  for  a  given  frequency.  The  calculation  is  performed  by 
integrating  along  a  ray  path  the  whistler  noise  emitted  at  each  point  and 
amplified  in  its  propagation  along  the  path,  and  then  summing  over  all  ray 
paths  passing  through  the  observation  point.  The  formula  is 


P(u,khsobj)  =  £  F  exp 

J 

(11) 


where  E  is  the  incoherent  emission  rate  and  F  is  a  geometric  factor  resulting 
from  the  divergence  of  ray  paths.  In  the  exponential  amplification  factor,  the 
growth  rate  y  is  calculated  using  a  detailed  model  of  the  auroral  electron  flux, 
while  the  group  velocity  vg  and  the  path  s  are  calculated  using  a  model  for 
the  plasma  in  the  auroral  region.  In  a  series  of  studies,  Maggs  and  Maggs 
and  Lotko  [31]  have  shown  that  the  features  of  the  observed  auroral  hiss 
spectrum  can  be  explained  by  this  model.  The  work  culminated  in  a  study 
(32]  in  which  a  detailed  fit  to  an  observed  hiss  spectrum  was  made  using 
simultaneous  rocket  data  for  the  electron  flux.  Because  the  plasma  density 
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Figure  2.  One  view  of  a  three-dimensional  perspective  plot  of 
the  trajectories  of  four  resonance-cone  waves,  with  frequencies  of 
0.1,  0.2,  0.4,  and  0.7  of  the  plasma  frequency  at  the  altitude  at 
which  they  were  launched.  The  highest  frequency  mode  is  the 
one  which  propagates  farthest  from  the  arc.  The  scale  of  the  plot 
and  the  orientation  of  the  plasma  density  gradient  are  shown  in 
the  inset. 

profile  and  composition  cannot  be  measured  along  the  entire  length  of  any  ray 
path,  the  plasma  profile  could  not  be  specified  by  observation.  It  was  found, 
however,  that  a  reasonable  choice  for  the  plasma  profile  produced  a  theoreti¬ 
cal  hiss  spectrum  which  agreed  well  with  the  observed  spectrum,  predicting 
both  the  frequency  of  the  spectral  peak  at  a  few  times  c o^,  as  well  as  the 
peak  inter. rity  a"d  spectral  shape. 
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It  is  expected  that  when  the  intensity  of  the  hiss  is  high  enough  [331, 
quasilinear  modification  of  the  electron  beam  through  wave-particle  interac¬ 
tion  with  the  hiss  will  become  important.  One  way  to  address  this  problem  is 
to  formulate  a  set  of  equations  describing  the  evolution  of  the  parameters  of  a 
beam  population  subject  to  wave-particle  interaction  with  the  waves  the  beam 
excites  as  well  as  the  other  influences  on  the  beam  [34],  Another  approach  is 
to  folio'-/  the  evolution  of  a  beam-plasma  system  self-consistently  using  a 
plasma  simulation.  We  will  follow  the  latter  approach  in  the  next  section  as 
we  begin  to  explore  the  consequences  of  wave-particle  interaction  with 
intense  auroral  hiss. 

IH.  LOCAL  VELOCITY  DIFFUSION  THEORY 


Let  us  first  examine  the  wave-particle  interaction  process  in  the  simpler 
situation  of  a  homogeneous  medium.  With  a  basic  understanding  of  the  phy¬ 
sics  involved  in  the  particle  heating  process,  we  will  then  turn  in  the  next  sec¬ 
tion  to  the  treatment  of  ion  heating  in  the  more  complicated  geometry  of  the 
suprauroral  region,  in  which  auroral  hiss  turbulence,  the  inhomogeneous 
geomagrutic  field,  and  possible  parallel  electric  fields  all  play  a  role  in  deter¬ 
mining  the  observed  panicle  fluxes.  In  the  presence  of  a  turbulent  electric 
field  spectrum,  the  resonant  (irreversible)  effect  on  the  ion  velocity  distribu¬ 
tion  function  /  is  adequately  described  by  a  diffusion  equation  [35], 


=  _a_ 

dt  dv 


(12) 


The  velocity  diffusion  tensor  D  depends  on  the  spectral  density  of  the 
electric-field  turbulence  as  a  function  of  frequency  and  wave-vector.  If  we 
specialize  to  the  study  of  ion  interactions,  and  use  the  unmagnetized  trajectory 
approximation,  we  have 

[-—■  f-iLL  l£  l2(*,o))  it  5(0)-* -v  )  (13) 

J  2n  J  fbtri3  lc2 


D  = 


JL 

m 


where  q  and  m  are  the  ion  charge  and  mass,  w  is  the  frequency,  k  is  the 
wave  vector,  an,-4  l£  1 2(k  ,co)  is  the  spectral  density  of  the  magnitude  of  the 
electric  field,  which  is  assumed  to  be  electrostatic  in  nature.  Observations  of 
auroral  hiss  tell  us  something  about  the  electric -field  spectral  density  which 
allows  us  to  simplify  this  picture.  The  largest  amplitudes  of  the  fields  are 
found  at  frequencies  which  are  only  a  few  times  larger  than  the  lower-hybrid 
resonance  frequency,  .  The  whistler  resonance-cone  dispersion  relation 
tells  us  that  the  wave  vectors  of  these  waves  are  directed  nearly  perpendicular 
to  the  geomagnetic  field.  Thus,  the  largest  contribution  to  the  integral  over 
the  spectral  density  for  D  comes  for  the  perpendicular  component 
Di=kLD  -kL/k2.  This  explains  the  essentially  transverse  heating  of  the  ions 


observed  in  ion  conic  events,  and  means  that  we  can  write  the  local  evolution 
equation  for  the  ion  velocity  distribution  keeping  only  the  perpendicular  con¬ 
tribution 


a/  =  j__a_ 

9 1  vA  9v| 


l°l  8v, 


(14) 


where 


D ,  = 
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~  l£l2(/k,co)  rc  5(to-*  v  ).  (15) 
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Note  the  delta  function  in  the  integTand,  imposing  a  resonance  condition  that 
states  that  the  only  waves  which  will  cause  ions  of  velocity  v  to  diffuse  in 
velocity  are  those  waves  for  which  the  projection  of  the  phase  velocity  in  the 
direction  of  v  equals  the  ion  speed. 

If  we  attempted  to  use  an  auroral-hiss  spectrum  predicted  using  local 
linear  excitation,  the  resonance  requirements  would  pose  a  great  difficulty. 
The  linear  excitation  of  auroral  hiss  is  achieved  through  a  Landau  resonance 
with  the  accelerated  auroral  electrons,  whose  velocities  (corresponding  to 
energies  of  1  keV  or  more)  serve  to  give  us  an  estimate  of  the  parallel  phase 
velocities  of  the  waves.  By  exciting  oblique  modes,  whose  wave  vectors  per¬ 
pendicular  to  B  are  much  larger  than  their  parallel  components,  we  reduce  the 
transverse  resonant  velocity  by  a  factor  of  the  order  of  the  square  root  of  the 
mass  ratio,  (o/i  =  (<a/kt){k^k)  -  wb||  cos8,  but  that  is  not  enough.  If  we  con¬ 
sider  the  phase  velocities  of  linearly  excited  resonance-cone  waves  of  fre¬ 
quency  about  twice  (o^hr,  we  find  that  they  propagate  at  an  angle  given  by 
cos28  =  3me  lml .  With  an  electron  beam  of  energy  I  keV,  this  corresponds  to 
an  ion  energy  of  3  keV.  There  are  very  few  ions  of  this  energy  in  the  undis¬ 
turbed  ionosphere,  where  the  temperature  is  of  the  order  of  l  eV.  This 
difficulty  is  resolved  when  we  recognize  that  a  number  of  processes  which  act 
on  the  waves  after  their  excitation  tend  to  reduce  their  phase  velocities.  Pro¬ 
pagation  of  the  resonance-cone  waves  through  a  transverse  density  gradient 
toward  the  lower-hybrid  resonance  can  reduce  the  phase  velocities  of  the 
waves,  as  will  linear  mode  conversion  at  the  resonance,  and  there  is  also  evi¬ 
dence  (described  below)  for  nonlinear  mode-mode  coupling  processes  which 
can  effectively  aid  the  ion  acceleration  process  by  generating  short 
wavelength  modes.  To  study  the  nonlinear  evolution  of  the  electron  beam 
and  the  heating  of  the  ambient  plasma  in  a  self-consistent  way,  we  performed 
plasma  simulations  of  the  phenomena. 
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A.  Plasma  Simulations 


To  study  the  generation  of  the  turbulence  and  the  resulting  ion  and 
electron  acceleration,  we  simulated  the  phenomena  using  an  electrostatic  par¬ 
ticle  model  of  the  plasma  [36].  The  suprauroral  situation  beneath  a  field- 
aligned  potential  drop  was  modeled  by  allowing  a  weak,  energetic  (1  keV), 
warm  (100  eV)  electron  beam  traveling  along  the  magnetic  field  to  destabilize 
a  cool  (1  eV)  electron-ion  plasma.  All  species  are  represented  initially  as 
Maxwellian  velocity  distributions. 

We  begin  with  one-dimensional  simulations  [6],  where  we  can  study 
waves  propagating  at  only  one  angle  with  respect  to  the  magnetic  field.  To 
study  the  waves  near  the  peak  of  the  auroral  hiss  spectrum,  we  choose  the 
angle  of  our  simulation  direction  to  be  nearly  perpendicular  to  B , 
kB  -  cos©  -  yjme  lm, ,  allowing  waves  near  the  lower  hybrid  resonance  fre¬ 
quency  to  be  excited.  (We  of  course  cannot  choose  k  to  be  exactly  perpen¬ 
dicular  to  B ,  or  no  coupling  to  the  electron  beam  traveling  along  B  would  be 
possible.) 

As  in  our  analysis  above,  we  treat  the  ions  as  unmagnetized  in  the 
simulation,  although  runs  with  magnetized  ions  showed  basically  the  same 
results.  For  the  electrons  we  employed  the  guiding-center  drift  approximation 
[37].  In  a  one-dimensional  model,  panicle  cross-field  drifts  are  out  of  the 
plane  of  k  and  B  ;  in  this  case  an  electron’s  guiding-center  velocity  in  the  k 
direction  is  simply  the  projection  in  that  direction  of  its  velocity  along  B ,  if 
the  small  polarization  drift  is  neglected.  With  this  approximation,  the 
acceleration  of  a  guiding  center  in  the  k  direction  is  qe  E  (x  ,t )  cos2*)  tme , 
implying  that  the  electron  guiding  center  can  be  treated  as  a  particle  moving 
in  the  k  direction  with  an  effective  inertia  mt  /cos20.  At  the  angles  of  propa¬ 
gation  of  lower  hydrid  waves,  the  effective  inertia  of  the  electron  guiding 
center  is  comparable  to  the  ion  mass,  implying  that  electrons  and  ions  will 
respond  in  similar  ways  to  the  waves,  and  incidently  permitting  realistic  mass 
ratios  to  be  used  in  the  simulation. 

As  we  follow  the  evolution  of  the  plasma  in  the  simulation,  we  first 
find  the  exponential  growth  of  the  waves  linearly  excited  by  the  instability  of 
the  electron  beam.  Then  as  this  instability  saturates,  the  beam  velocity  distri¬ 
bution  broadens  and  flattens  as  energy  and  momentum  are  transferred  from  it 
to  the  other  species  and  the  waves.  The  evolution  of  the  fluctuation  spectum 
is  illustrated  in  Fig.  3.  In  this  figure  we  have  plotted  I E  l2(fc)  as  a  function  of 
k  at  four  times.  The  first  snapshot,  in  the  upper  left,  is  taken  during  the  stage 
of  linear  instability  and  shows  the  excitation  of  the  waves  in  the  narrow  range 
of  wave  vectors  of  the  linearly  unstable  modes.  The  next  snapshot,  in  the 
upper  right,  shows  the  spectrum  just  prior  to  saturation,  with  a  higher  ampli¬ 
tude  and  slightly  broader  range  of  k .  From  this  concentration  in  a  relatively 
small  region  of  k  space,  we  subsequently  see  the  wave  energy  redistribute 
itself  throughout  k  space  in  the  following  two  snapshots  on  the  bottom  of  the 
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Figure  3.  The  wave  spectrum  in  a  1-d  simulation.  In  each  panel 
the  spectral  density  of  the  electric  field  is  plotted  as  a  function  of 
wave  vector  at  a  fixed  time,  in  sequence  first  from  left  to  right, 
then  from  top  to  bottom.  The  wave  vectors  are  expressed  as  a 
multiple  of  the  fundamental  wave  vector  of  the  simulation  length. 

figure,  through  the  action  of  mode-coupling  processes.  The  results  at  the  end 
is  an  approximately  steady-state  spectrum  which  at  smaller  k  take  the 
Rayleigh-Jeans  form,  l£  l2(ifc)  =  constant. 

The  evolution  of  the  particle  velocity  distributions  is  illustrated  in  Figs. 
4  and  5.  Fig.  4  shows  two  snapshots  of  the  electron  distribution  of  parallel 
velocities,  one  showing  the  initial  velocity  distribution  and  the  other  showing 
the  distribution  following  saturation  of  the  instability.  The  ambient  (cool) 
electrons  and  the  beam  electrons  are  shown  as  separate  distributions,  although 
of  course  for  stability  analyses,  etc.,  it  is  the  sum  of  the  two  for  the  total 
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Figure  4.-  Electron  parallel- velocity  distribution  in  a  1-d  simula¬ 
tion.  The  dashed  curves  give  the  initial  distributions  of  the  am¬ 
bient  and  beam  electron  populations,  while  the  solid  curves  give 
the  distributions  following  the  end  of  the  simulation. 

distribution  which  is  relevant.  The  figure  shows  the  plateau  formation  on  the 
unstable  beam  distribution,  but  more  interestingly,  we  see  that  the  ambient 
electrons  have  been  heated  in  both  directions  with  respect  to  the  magnetic 
field.  The  acceleration  of  electrons  in  the  direction  of  the  beam  is  not  unex¬ 
pected,  because  if  produced  by  a  resonant  interaction  with  waves,  then  this 
heating  was  produced  by  waves  traveling  in  the  same  direction  as  the  beam, 
just  as  the  waves  that  are  linearly  excited  by  the  beam  are.  It  is  more  surpris¬ 
ing  that  nearly  as  many  electrons  were  accelerated  in  the  opposite  direction, 
by  waves  that  are  traveling  in  the  direction  opposite  to  the  beam.  Evidence 
of  such  parallel  acceleration  of  electrons  in  the  suprauroral  region  can  be 
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Figure  5.  The  ion  velocity  distribution  in  a  1-d  simulation.  The 
multiple  curves  give  the  ion  velocity  distribution  at  successive 
times  throughout  the  simulation,  showing  the  formation  of  the 
high-velocity  tails. 

found  in  the  satellite  observations  of  counter- streaming  electrons  [7],  those  of 
type  1,  in  which  there  is  no  indication  of  acceleration  by  a  steady  parallel 
electric  held.  An  analysis  of  the  magnitude  of  the  phase  velocities  needed  to 
resonantly  interact  with  the  ambient  electrons,  like  the  earlier  one  for  the 
ambient  ions,  shows  that  the  linearly  excited  waves  would  have  phase  veloci¬ 
ties  too  high  to  affect  the  ambient  electrons.  From  these  two  observations  we 
infer  that  the  nonlinear  process  we  noted  in  the  evolution  of  the  wave  spec¬ 
trum  is  also  aiding  the  electron  acceleration.  Fig.  5  shows  us  that  the  ion 
velocity  distributions  undergoes  the  formation  of  similar  high-energy  tails. 
Because  ti  '  ions  are  unmagnetized,  they  are  accelerated  in  the  direction  of 
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the  simulation  electric  fields,  nearly  perpendicular  to  B .  In  this  simulation, 
where  cos29  =  me  lmi ,  the  ambient  electrons  and  ions  are  heated  by  equal 
amounts. 

In  the  interest  of  simplicity,  this  first  set  of  simulations  [6]  were  one¬ 
dimensional  simulations,  in  which  the  waves  could  propagate  in  only  one 
direction  with  respect  to  the  ambient  magnetic  field.  A  number  of  interesting 
phenomena,  however,  can  be  realized  only  in  higher  dimension,  where  a 
range  of  propagation  angles  can  be  studied  simultaneously.  For  example, 
both  the  real  frequency  and  the  convective  linear  growth  rates  of  waves  on 
the  whistler  resonance-cone  dispersion  surface  are  a  function  of  the  angle  of 
propagation  with  respect  to  the  magnetic  field.  Mode-coupling  processes  can¬ 
not  be  studied  in  complete  generality  in  one  dimension,  because  wave  vectors 
are  there  constrained  to  be  aligned  in  the  same  direction.  Finally,  because 
electrons  and  ions  interact  with  different  efficiencies  to  waves  of  different  fre¬ 
quencies,  an  estimation  of  the  relative  effectiveness  of  electron  and  ion 
acceleration  requires  a  simulation  in  which  waves  of  an  appropriate  range  of 
frequencies  can  be  excited. 

To  address  these  questions,  a  scries  of  two-dimensional  simulations 
were  performed  [38].  Full  electron  dynamics  were  employed  and  the  ions 
were  treated  as  magnetized  particles,  although  to  compress  the  dynamical 
range  of  frequencies  in  the  simulation  a  smaller  than  actual  mass  ratio  was 
necessary.  (For  the  lightest  ions  the  ratio  lme  was  25.)  An  angular  cutoff 
was  used  to  simulate  the  effect  of  the  preferential  convective  saturation  of  the 
parallel  propagating  modes.  It  was  found  that  a  broadband  frequency  spec¬ 
trum,  spanning  the  allowed  range  of  resonance-cone  modes,  was  excited.  The 
nature  of  the  electron  and  ion  acceleration,  however,  was  not  qualitatively 
changed.  Figs.  6  and  7  show  the  electron  and  ion  velocity  distributions  from 
a  two-dimensional  simulation  run.  A  high  energy  tail  in  the  direction  perpen¬ 
dicular  to  the  magnetic  field  is  observed  on  the  ion  velocity  distribution  plot¬ 
ted  as  a  function  of  energy,  while  the  ambient  electrons  are  heated  with  the 
counter-streaming  tails  as  in  the  one-dimensional  simulations.  The  amount  of 
heating  of  the  two  ambient  populations  did  depend  on  the  choice  of  the 
resonance-cone  angles  permitted  in  the  simulation,  but  because  a  range  of 
resonance -cone  modes  were  always  excited  in  the  two-dimensional  simula¬ 
tions,  the  amount  of  heating  was  less  sensitively  dependent  on  the  angle  of 
propagation  than  it  would  have  been  in  a  one-dimensional  simulation. 

In  one  set  of  simulations  a  mixture  of  ions  species  was  included,  using 
both  light  ions  representing  the  hydrogen  of  the  ionosphere,  and  heavy  ions 
representing  oxygen.  It  was  found  that  the  heating  of  the  heavier  species  was 
negligible.  This  can  be  understood  as  a  consequence  of  the  velocity  depen¬ 
dent  resonance  condition  for  interaction  with  the  waves.  The  waves  are 
excited  at  very  high  phase  velocities  by  the  electron  beam,  and  then  are  sub¬ 
ject  to  mode-coupling  processes  that  reduce  their  phase  velocities.  These 
processes  continue  until  the  phase  velocities  of  the  waves  are  sufficiently 
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Figure  6.  Electron  velocity  distribution  in  the  2-d  simulation. 

The  dot-dash  curves  give  the  initial  parallel  velocity  distributions 
of  the  ambient  and  beam  electron  populations,  while  the  dashed 
and  solid  curves  give  the  distributions  at  two  times  near  the  end 
of  the  simulation  run. 

degraded  that  the  cool  ambient  plasma  species  can  begin  to  interact  with 
them.  The  species  of  particles  with  the  largest  velocities  will  be  able  to 
interact  with  the  waves  first.  By  absorbing  the  wave  energy,  this  first  ion 
species  prevents  the  waves  from  undergoing  further  degradation  in  phase 
velocity  and  isolates  the  species  with  smaller  velocities,  namely  species  with 
heavier  masses,  from  interaction  with  the  waves. 
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Figure  7.  Ion  velocity  distribution  in  the  2-d  simulation.  The 
dot-dash  line  gives  the  initial  velocity  distribution  (a  Maxwellian) 
of  the  ion  population  as  a  function  of  the  perpendicular  energy 
per  unit  mass,  while  the  solid  curve  give  the  distributions  near  the 
end  of  the  simulation  run. 

B.  Nonlinear  Theory 

It  is  a  common  observation  that  in  plasma  simulations,  with  their  full 
treatment  of  nonlinear  plasma  phenomena,  many  plasma  waves  are  excited 
which  linear  theory  would  predict  to  be  damped  or  only  weakly  growing. 
Interaction  of  ions  and  electrons  with  such  modes  at  low  phase  velocities  can 
produce  the  energetic  tails  found  on  the  their  velocity  distributions.  One  of 
the  first  nonlinear  processes  to  act  as  waves  begin  to  grow  in  the  plasma 
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involves  the  coupling  of  only  three  (or  four)  waves.  In  weak  turbulence 
theory  [35]  the  coupling  of  three  waves  requires  that  the  waves  satisfy  the 
resonance  conditions 


co0  =  cuj  +  co2  (16) 

*o  =  *t+*2  (17) 

where  (c o0,Jt0)  are  the  frequency  and  wave  vector  of  a  parent  wave,  and 
((Oi Jfc  i)  and  (o>2,ifc  2)  are  the  corresponding  quantities  for  the  two  waves  into 
which  the  parent  wave  decays.  Let  us  suppose  that  the  parent  wave  is  a 
lower  hybrid  wave  linearly  excited  by  the  electron  beam  instability.  If  the 
sideband  wave  (0)i^|)  is  another  lower  hybrid  wave,  then  to  satisfy  the  reso¬ 
nance  conditions  the  beat  wave  must  be  of  very  low  frequency,  of  the 

order  of  the  difference  in  thermal  corrections  to  the  frequencies  of  the  two 
lower  hybrid  waves.  The  beat  mode  may  be  either  a  heavily  damped  acoustic 
quasimode  or  the  near-zero  frequency  oscillating  two-stream  instability  [39]. 
The  phase  velocities  of  the  sideband  modes  are  (if  all  three  wave  vectors  are 
aligned  and  allowance  is  made  for  the  two  sign  possibilities  for  k2 ) 

0)!  _  (OpiO)^  (Op  _  +  COp 

*i  ”  k(ytk2  k(ytk2  ~  k2 

If  /fc2»*o  the0  G>i/*i  can  **  much  smaller  than  the  phase  velocity  of  the 
parent  wave,  ^kQ,  allowing  resonant  interaction  with  the  ions  and  ambient 
electrons.  Note  that  sidebands  can  be  excited  with  phase  velocities  of  both 
signs  with  nearly  equal  magnitudes,  accounting  for  the  symmetry  of  the  hot 
tails  on  the  velocity  distributions. 

A  dispersion  relation  for  this  nonlinear  process  can  be  derived.  Details 
can  be  found  .in  Ref.  [6]  and  references  therein.  From  this  dispersion  rela¬ 
tion.  the  wave-amplitude  threshold  for  the  decay  process  can  be  calculated,  as 
well  as  the  rate  of  the  process  when  the  amplitude  of  the  parent  wave  exceeds 
this  threshold.  Such  an  analysis  of  lower  hybrid  parametric  decay  under 
suprauroral  conditions  has  been  performed  by  Koskinen  [40].  He  reports  that 
decay  through  nonresonant  quasimodes  is  one  of  the  dominant  processes, 
although  other  decay  routes,  such  as  decay  into  a  lower  hybrid  mode  and  an 
ion  Bernstein  mode,  are  possible.  He  found  that  the  threshold  of  wave  ampli¬ 
tude  for  decay  may  be  only  a  few  millivolts  per  meter  or  smaller,  insuring 
that  this  mode  coupling  process  can  be  effective  in  the  suprauroral  region. 

Visualizing  the  nonlinear  evolution  of  the  plasma  wave  spectrum  in 
terms  of  parametric  decay  processes  is  a  good  way  to  describe  the  initial 
stages  of  the  evolution,  but  as  soon  as  the  sideband  modes  have  grown  to 
appreciable  magnitude,  they  too  will  be  susceptible  to  decay,  and  a  more  gen¬ 
eral  approach  is  necessary.  In  addition,  to  verify  our  interpretation  of  the 
simulatio.  results  we  would  like  to  be  able  to  calculate  the  evolution  of  the 
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turbulent  spectrum  and  its  effect  on  the  particle  populations. 

A  nonlinear  kinetic  equation  that  describes  the  evolution  of  the  wave 
spectrum  subject  to  mode  coupling  effects  has  been  derived.  We  refer  to  Ref. 
[6]  for  the  details  of  the  derivation,  but  the  essence  of  the  calculation  can  be 
summarized  by  saying  that  we  begin  with  the  warm  fluid  equations  for  the 
electron  and  ion  plasma  species  (for  simplicity,  since  the  lower  hybrid  waves 
are  supported  in  that  approximation).  Assuming  that  the  electric  field  and 
other  quantities  consist  of  separate  components  oscillating  at  the  lower  hybrid 
and  acoustic  frequencies,  we  retain  only  the  dominant  nonlinear  terms  and 
derive  an  equation  describing  the  time  variation  of  the  amplitude  of  the 
lower-hybrid  field  component.  In  one  dimension,  the  result  is 

f =  °  (W 

where  E(x,t)  is  the  amplitude  of  the  electric  field  oscillating  at  <0^, 
assumed  to  vary  slowly  with  time.  The  spectral  density  at  co ^  is  propor¬ 
tional  to  l£  I2.  Cc  and  Cd  are  coefficients  for  the  coupling  and  dispersion 
processes,  respectively,  expressed  in  terms  of  plasma  frequencies,  tempera¬ 
tures,  charges,  and  masses  of  the  species  of  the  plasma.  It  is  of  the  same 
form  as  the  nonlinear  Schroedinger  equation  derived  to  describe  Langmuir 
turbulence  [41].  Perhaps  the  most  fascinating  results  of  the  nonlinear 
Schroedinger  equation  are  the  solitary  wave  solutions  (for  the  envelope  of  the 
lower  hybrid  waves  in  this  case)  which  it  supports.  In  one  dimension,  and  in 
the  absence  of  dissipation,  these  sharply  peaked  wave  packets  propagate 
without  alteration  because  the  nonlinear  refraction  balances  the  effects  of 
dispersion.  With  Landau  damping  present,  the  formation  of  solitary  waves 
through  the  modulational  instability  of  an  intense  lower  hybrid  wave  is  inter¬ 
rupted  when  the  wave  energy  has  spread  to  wave  vectors  large  enough  to 
become  susceptible  to  damping.  To  be  able  to  describe  the  effects  of  such 
dissipation  on  the  waves,  as  well  as  the  resulting  effects  on  the  particles,  we 
introduce  kinetic  effects  into  the  model,  following  Thomson  et  al.[ 42].  To 
introduce  Landau  damping  into  the  wave  equation,  we  Fourier  transform  spa¬ 
tially  and  add  the  linear  dissipative  term. 

%  =  -  ±Cdk2Ek  +iCc  I  Ek-E^Ek,k-_k.  (20) 

dr  2  *'jr 

where  the  linear  damping  or  growth  rate  is  given  by  the  quasilinear  expres¬ 
sion 

/l%  ft  ^  01 p s  s  |  , 

The  effect  of  the  wave-panicle  interaction  on  the  particles  is  described  by  the 
resonant  quasilinear  diffusion  equation 
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(22) 


a/, 

dt 


where  in  one  dimension  the  velocity  diffusion  coefficient  is 


Ds(v) 


I Ek  !2tc5(  03(k) 


-k  v  ) 


(23) 


The  set  of  equations,  Eqs.(20)  and  (22),  with  Eqs.(21)  and  (23), 
comprise  our  theoretical  model,  a  generalization  of  the  quasilinear  equations 
to  include  the  effects  of  mode  coupling.  We  use  this  set  of  equations  now  to 
study  the  saturation  of  the  bump-on-tail  instability  that  excites  the  intense 
long-wavelength  lower  hybrid  modes.  The  calculation  begins  with  initial  con¬ 
ditions  corresponding  to  those  with  which  we  started  the  one-dimensional  par¬ 
ticle  simulations,  and  the  equations  are  integrated  forward  in  time  to  study  the 
evolution  of  the  system.  Fig. 8  compares  the  evolution  of  the  total  lower 
hybrid  wave  energy  in  the  one-d  simulation  and  in  the  theoretical  model.  The 
dashed  curve  in  the  figure  gives  the  evolution  of  the  electrostatic  energy  in 
the  simulation,  while  the  solid  curve  gives  the  result  from  the  theoretical 
model.  We  experimented  with  the  model  by  removing  the  mode-coupling 
term  from  the  equation;  the  dot-dash  curve  gives  the  results  for  the  wave 
energy  in  this  case.  We  see  at  initial  times  the  exponential  growth  in  energy 
associated  with  the  linear  instability.  Once  the  threshold  for  mode-coupling 
processes  has  been  exceeded,  a  turbulent  steady-state  in  the  waves  is  soon 
reached,  while  energy  is  transferred  to  the  particle  populations  (although  their 
energies  are  not  shown  here).  The  level  of  energy  in  the  saturated  state  is 
very  similar  to  the  level  observed  in  the  simulation.  In  the  solution  without 
the  mode-coupling  term,  on  the  other  hand,  the  wave  energy  continues  to 
grow,  although  the  growth  rate  slows  as  the  beam  velocity  distribution  pla¬ 
teaus  at  the  phase  velocities  of  the  fastest  growing  modes.  Without  the 
mode-coupling  term,  the  energy  transferred  to  the  ambient  particle  popula¬ 
tions  is  negligible. 

The  detailed  ion  velocity  distribution  in  the  mode  coupling  calculation 
shows  the  high  velocity  tails  characteristic  of  the  particle  simulation.  Fig.  9 
shows  a  snapshot  of  the  ion  velocity  distribution  following  the  saturation  of 
the  instability,  illustrating  the  formation  of  the  tails.  The  solid  curve  gives 
the  result  from  the  theoretical  model,  while  the  dashed  curve  in  the  figure 
gives  the  comparable  result  from  the  particle  simulation.  We  conclude  from 
these  comparisons  that  our  theoretical  model  with  mode-coupling  does  satis¬ 
factorily  describe  the  phenomena  observed  in  the  particle  simulations.  While 
the  simple  quasilinear  velocity  diffusion  model  is  adequate  to  describe  the 
effect  of  the  turbulence  on  the  particles,  a  full  nonlinear  treatment  is  neces¬ 
sary  to  describe  the  evolution  of  the  turbulence. 
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Figure  8.  The  energetics  of  the  lower-hybrid  waves  in  the 
theoretical  model.  The  dashed  curve  gives  as  a  function  of  time 
the  electrostatic  energy  of  the  waves  in  a  1-d  particle  simulation 
run,  while  the  solid  curve  gives  the  prediction  of  this  quantity 
from  our  theoretical  model.  The  dot-dash  curve  gives  the  theoret¬ 
ical  result  when  the  mode-coupling  term  is  removed  from  the 
model.  Time  is  in  units  of  . 

C.  Diffusion  Models 

Although  the  simulation  results  can  serve  as  a  qualitative  guide  to  the 
phenomena  observed  in  space,  the  situation  in  the  auroral  region  is  more 
complicated  than  we  can  practically  study  using  simulations.  Important 
phenomena  occur  on  a  wider  range  of  spatial  and  temporal  scales  than  self- 
consistent  numerical  models  can  accomodate  simultaneously.  For  success  in 
modeling  the  observations  we  must  hope  to  extract  from  the  simulations  a 
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Figure  9.  The  ion  velocity  distribution  in  the  theoretical  model. 

The  dashed  curve  gives  the  ion  velocity  distribution  at  the  end  of 
a  particle  simulation  run,  while  the  solid  curve  gives  the  theoreti¬ 
cal  prediction  for  the  ion  distribution. 

description  of  the  micro-physics  of  the  wave-particle  interaction,  to  incor¬ 
porate  into  a  model  of  the  larger  scale  variations  in  the  suprauroral  region. 
For  the  heating  of  the  ambient  particle  species,  the  model  is  the  velocity 
diffusion  equation,  to  which  the  simulations  provide  a  model  for  the  aspects 
of  the  spectral  density  of  the  electric  fields  which  cannot  be  measured 
directly. 

Because  of  the  variety  of  circumstances  to  be  found  in  the  suprauroral 
region,  w  ?nusi  also  expect  to  employ  wave  observations  in  our  specification 


of  particle  heating  rates  to  insure  the  applicability  of  our  models.  Satellites 
rarely  measure  wavevector  information  for  the  electric  field  spectral  density: 
the  observed  spectral  density  of  any  component  of  the  electric  field  is  a  func¬ 
tion  only  of  frequency  (which  unfortunately  includes  the  possibility  of  alias¬ 
ing  of  spatial  structure  if  the  spacecraft  is  moving  through  an  inhomogeneous 
medium).  Thus,  if  the  satellite  antenna  system  were  equally  responsive  to  all 
wavelengths,  the  spacecraft  receiver  would  measure 


|£|2(*fco)  =  l£  l2(Ci>)  . 

J  (2k)3 


One  can  obtain  a  rough  order  of  magnitude  of  the  transverse  velocity 
diffusion  rate  by  using  the  observed  value  of  I E  l2(co)  at  some  frequency  near 
the  lower  hybrid  frequency  in  the  order-of-magnitude  evaluation  of  Eq.(15) 

l2 

DL~  -9-  I£12(oj)  (25) 


and  this  is  frequently  all  that  can  be  done  to  estimate  the  diffusion  rate.  One 
cannot,  however,  describe  the  velocity  dependence  of  Di  using  a  description 
this  crude,  and  the  high-energy  tail  structure  of  the  particle  fluxes  seen  both 
in  space  and  the  simulations  requires  a  velocity  dependent  diffusion  rate. 
This  is  the  point  where  our  simulation  results  can  suggest  the  form  for  a 
model  of  the  unobservable  wave  vector  specta  of  the  resonance-cone  waves. 

To  use  the  frequency-spectrum  observations  to  estimate  the  velocity 
diffusion  rate,  let  us  first  condense  all  our  ignorance  of  the  spectral  density 
into  one  function,  the  distribution  of  spectral  density  with  respect  to  the  wave 
vector,  given  a  specified  frequency,  S(k  I  co) 

lfl2(*,co)  =  lf!2(co)  S(k  I  co)  (26) 


where 


S(k  I  co)  a  l 


For  frequencies  near  the  lower  hybrid  frequency,  we  can  write  the  most 
important  functional  dependence  of  5  on 

S(k  !oj)  =Si(*i)5,(*fl)  .  (28) 

Although  in  a  turbulent  plasma  we  cannot  assume  in  general  that  the  frequen¬ 
cies  and  wavevectors  of  the  modes  in  the  plasma  satisfy  a  dispersion  relation 
perfectly,  for  this  application  this  assumption  does  not  cause  a  great  error, 
and  we  write 


5,,(*)  =  2nd(k  -k#) 


where 
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*yo  2  kt(kL,(o)  =  ±A(co)  kL 


(30) 


-  ,  ,  (U2-<Oujd 

A2(CO)  =  (1+0)2  /Q2) - __  (3l) 

C0*-(ir 

Let  us  substitute  these  expressions  into  the  formula  for  the  transverse  velocity 
diffusion  rate 


(*r*i)2 


SL(kL)  5(o>-ti  vi-^l0v„)(32) 


For  waves  near  the  lower  hybrid  frequency,  klt0  is  much  less  than  k±,  and  the 
resonance  condition  is  essentially  one  on  the  perpendicular  velocity.  If  we 
assume  that  SL(ki)  is  isotropic  in  the  plane  perpendicular  to  B ,  we  can  do  one 
more  integration  of  this  expression,  and  we  get 


0,  =  -r 
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To  proceed  further,  we  need  a  model  for  Crew  and  Chang  [43] 

suggested  the  following  form,  to  describe  the  general  situation  where  Sj^  peaks 
near  a  wave  vector  k0 


SL(kL) 


(kL/kor 

i +(kL/k0y 


(34) 


In  this  expression  the  free  parameters  are  the  wave  vector  k0  and  the  two 
spectral  indices  p  and  v.  For  the  normalizing  integral,  Eq.(27)  to  be  defined, 
we  need  v>p+3,  and  for  D ^  to  be  nonsingular,  p>0.  This  form,  it  turns  out, 
fits  the  body  of  the  wave  vector  spectra  measured  in  the  one-dimensional 
simulations  quite  well,  and  we  will  use  it  to  parametrize  the  wave  vector 
spectra  to  calculate  the  velocity  diffusion  rate.  Crew  and  Chang  [43]  next 
showed  that  to  a  good  approximation,  the  diffusion  coeficicnt  one  obtains 
using  this  form  has  the  following  velocity  dependence 


DL(VL)  oc 


(Vj/Vj^-3 
1  +n(M.v)(vi/v0  )*■** 


(35) 


Here  q  is  a  numerical  constant,  dependent  on  p  and  v,  and  v0  is  a  charac¬ 
teristic  velocity  near  which  the  diffusion  rate  is  a  maximum;  if  the  frequency 
spectrum  peaks  near  a  frequency  o0,  then  k0  =0)o/vo.  We  expect  that  va  is 
of  the  order  of  the  perpendicular  phase  velocity  of  the  waves  excited  by  the 
auroral  electron  beam,  or  smaller  if  mode-coupling  processes  have  been 
effective  n  transferring  wave  energy  to  shorter  wavelength,  lower  phase- 
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velocity  waves.  At  Vj  goes  to  zero,  D[  goes  to  zero  as  viv~M'  3 ;  at  velocities 
much  larger  than  v0 ,  Dy  falls  as  vf3 . 

The  peaked  nature  of  D |  reflects  the  way  in  which  wave-particle 
interaction  leads  to  the  formation  of  high-energy  tails  on  the  velocity  distribu¬ 
tions,  rather  than,  say,  simple  bulk  heating.  Because  DL  falls  to  zero  at  small 
velocity,  the  bulk  of  the  velocity  distribution  is  unaffected  by  interaction  with 
the  waves.  At  higher  velocity  begins  to  become  appreciable,  and  the 
diffusion  process  begins  to  affect  particles  in  the  distribution.  Because  there 
are  initially  more  particles  at  lower  velocity  than  at  higher  velocity,  the  net 
flux  in  velocity  space  is  toward  higher  velocities,  leading  to  the  formation  of 
high-energy  tails.  Crew  and  Chang  [43]  demonstrated  this  by  developing  an 
asymptotic  method  for  analytically  solving  the  velocity  diffusion  equation 
with  the  peaked  diffusion  coefficient  and  presenting  the  solutions. 

IV.  ION  CONIC  FORMATION 


We  must  now  consider  how  our  local  description  of  particle  heating  by 
wave-particle  interaction  with  intense  auroral  hiss  can  be  applied  in  the 
suprauroral  region,  where  ion  conics  and  counter-streaming  electrons  are 
observed  in  the  vicinity  of  auroral  arcs.  The  basic  point  of  view  for  our  stu¬ 
dies  of  this  problem  is  the  test-particle  approach,  in  which  the  motion  and 
heating  of  particles  are  studied  in  specified  global  models  of  the  geomagnetic 
field,  the  steady  electric  fields,  and  the  spectral  density  of  the  turbulent 
fluctuating  electric  fields.  The  justification  for  this  approach  lies  in  the  expec¬ 
tation  that  the  perturbation  caused  by  these  particles  will  be  small  because  the 
number  of  ionospheric  ions  and  electrons  which  are  affected  by  this  tur¬ 
bulence  will  be  small,  being  those  particles  which  can  meet  the  velocity  reso¬ 
nance  conditions  for  interaction.  We  will  first  outline  the  method  of  the 
model,  and  then  go  on  to  describe  several  applications. 

In  the  absence  of  wave-particle  interactions,  the  evolution  along  a 
goemagnetic  field  line  of  the  ion  velocity  distribution,  / ,  would  be  described 
by  the  Liouville  equation  [44].  If  we  use  the  guiding  center  approximation, 
the  equation  is  the  following: 
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where  s  is  a  coordinate  denoting  position  along  a  geomagnetic  field  line,  B  is 
the  magnitude  of  the  geomagnetic  field,  and  q  and  m  are  the  charge  and 
mass  of  an  ion.  We  ignore  possible  cross-field  drifts  caused  by  a  perpendicu¬ 
lar  electric  field,  Eit  and  only  retain  the  effects  of  a  possible  parallel  field,  £,. 
The  ion  motion  is  then  constrained  to  lie  along  a  field  line.  The  ions  move 
adiabatically,  conserving  their  magnetic  moment  and  energy  in  the  static 
geomagnetic  field  and  parallel  electric  field.  The  component  of  the 
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gravitational  field  that  lies  along  B  can  easily  be  included  together  with  the 
parallel  electric  field,  although  its  effect  is  less  important  for  the  light  hydro¬ 
gen  ions  of  this  study  than  it  would  be  for  oxygen  conics. 

In  the  presence  of  a  transverse  heating  process,  the  effect  on  the  ions  of 
their  interaction  with  the  turbulence  is  adequately  described  by  a  diffusion 
equation. 


df_  _  1  3 

dt  vA 
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(37) 


where  the  left-hand  side  of  the  equation  contains  the  terms  of  the  Liouville 
equation  and  Dl  is  the  quasilinear  velocity  diffusion  rate  perpendicular  to  the 
geomagnetic  field  obtained  in  the  previous  section. 


A.  The  Mean-Panicle  Equations 


Before  solving  the  spatially  dependent  kinetic  equations  for  the  ion  dis¬ 
tribution  function,  it  is  useful  to  have  an  estimate  of  the  relevant  ranges  of 
energy  and  pitch  angle  that  will  be  important  in  the  problem  at  hand.  For 
such  an  analysis,  a  simple  model  described  by  the  mean-particle  equations 
can  be  derived  from  the  terms  in  the  kinetic  equation.  This  pair  of  ordinary 
differential  equations  describe  the  evolution  of  the  perpendicular  and  parallel 
energies  of  a  typical  particle  within  the  ion  distribution.  The  equations  arc 
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where  W^Vvnv^  and  are  the  "typical”  parallel  and  perpendicular  ener¬ 
gies,  respectively.  The  perpendicular  heating  rate  due  to  wave-particle 
interaction,  Wlturb,  is  estimated  from  the  right-hand  side  of  the  kinetic  equa¬ 
tion.  For  example,  if  we  assume  that  the  velocity  distribution  is  Maxwellian 
and  that  the  diffusion  rate  DL  is  independent  of  velocity,  then  the  appropriate 
moment  of  Eq.(14)  gives  the  heating  rate  2mDi .  Chang  and  Coppi  [5]  used 
these  equations  to  demonstrate  that  lower  hybrid  waves  could  be  an  effective 
heating  mechanism  for  the  production  of  ion  conics. 

B.  The  Monte  Carlo  Method 


To  predict  the  form  of  the  distribution  function,  it  is  necessary  to  return 
to  the  kinetic  equation  and  solve  it  in  detail.  Even  in  the  simplest  situation 
where  a  steady  state  solution  is  sought,  the  problem  is  a  three  dimensional 
one,  with  e  independent  variables  position  along  the  field  line,  perpendicular 


velocity,  and  parallel  velocity.  In  addition,  the  specified  fields  in  a  problem 
will  not  necessarily  be  found  in  forms  that  make  the  problem  amenable  to 
analysis.  Test-particle  simulation  techniques  offer  a  flexibile  and  robust  way 
of  attacking  such  complicated  problems.  We  will  solve  the  kinetic  equation 
for  the  ion  velocity  distribution  using  a  Monte  Carlo  model  which  was 
developed  to  investigate  problems  of  this  kind  [8].  From  an  initial  distribu¬ 
tion  in  velocity  and  space,  the  calculation  of  the  evolution  of  the  distribution 
proceeds  by  following  the  trajectories  of  a  large  number  of  ions  with  time. 
Between  the  velocity  perturbations  caused  by  interaction  with  the  waves,  it  is 
assumed  that  the  ions  travel  in  the  static  geomagnetic  field  with  constant 
energy  and  magnetic  moment. 

The  wave-panicle  interactions  are  taken  into  account  by  perturbing  the 
ion  velocities  with  random  increments  Av  in  each  small  time  step,  Ar.  To 
see  how  these  perturbations  affect  the  velocity  distribution,  consider  the  situa¬ 
tion  in  which  the  steady  fields  are  absent  and  the  perturbations  act  alone.  If 
PAj(Av)d3Av  is  the  probability  of  a  change  of  velocity  Av  during  a  time 
interval  A r,  then  the  effect  of  velocity  perturbations  over  A r  on  the  distribu¬ 
tion  /(v,r)  is  the  convolution 

/(v,r+Ar)  =  J/(v-Av,r)/>Ai(Av)d3Av  (40) 

If  we  perform  the  Fokker-Planck  expansion  of  Eq.(40)  for  the  limit  that  A I 
goes  to  zero,  we  recover  the  local  velocity  diffusion  equation,  Eq.(12),  if  we 
choose  P^(Av)  to  be  the  gaussian  distribution  with  moments 

<Av>  =  D Ar  (41) 

and 

<Av  Av  >  =  2  D  Ar  .  (42) 

In  the  calculation,  this  probability  distribution  is  sampled  using  random 
numbers.  To  simulate  transverse  velocity  diffusion,  the  two  components  of 
Av  in  the  plane  perpendicular  to  B  are  chosen  to  be  gaussian  random  vari¬ 
ables  such  that  <Avx>  =  (v^/vj  A t  dDjdv^,  <Avy>  =  (vy/v{)A t  dD^/dv^,  and 
<Av/>  =  <Av,2>  =  ID  fit .  The  vector  Av  is  added  to  the  ion’s  transverse 
velocity  assuming  a  random  orientation  with  respect  to  the  gyrophase. 

Because  we  are  studying  the  heating  and  acceleration  of  ions  of  ionos¬ 
pheric  origin,  the  initial  velocity  distribution  of  the  ions  is  chosen  to  be  a 
Maxwellian  distribution  at  a  characteristic  topside  ionospheric  temperature,  1 
eV,  and  for  simplicity  we  will  generally  start  the  particles  all  at  one  altitude. 
As  ions  pass  the  observation  points,  statistics  on  their  velocities  are  accumu¬ 
lated  to  calculate  the  steady  state  ion  velocity  distribution  at  these  points. 
The  velocity  distribution  of  the  ions  at  their  initial  location  s0  can  be  written 
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(43) 


/(Vi,vii40)  =  -jr  Z  5(vi2-vi)  5(v,rvllw) 

where  viW  and  v|(^  are  the  initial  velocity  components  of  an  ion,  chosen  so 
that  /  equals  the  desired  initial  distribution,  and  the  density  of  the  N  ions  is 
normalized  to  n0.  After  the  ions  have  followed  their  trajectories  to  the  point 
s ,  the  velocity  distribution  at  that  point  can  be  calculated.  It  is 

f(vL,v«,s)  =  ~  I  -7s-  “  SO^-v  ?)  8(VN-VW)  (44) 

XN  ,=1  v»  Bo 

where  and  are  the  components  of  the  ion  velocity  at  s .  The  weighting 
factor  under  the  sum  is  a  consequence  of  the  invariance  of  the  parallel  parti¬ 
cle  flux  n(s  )<v  «(s  )>/fi  (5)  in  a  steady  state[43].  Physically,  the  weighting 
factor  expresses  two  effects:  v  ilio  /v  w  describes  the  reduction  in  phase-space 
density  that  occurs  when  particles  get  spread  out  along  a  field  line  when  v, 
changes,  while  Bs/B0  describes  the  dilution  due  to  the  increasing  flux-tube 
area  as  B  decreases  with  increasing  s .  The  detailed  velocity  distribution, 
Eq.(44),  contains  more  information  than  we  really  need  or  want;  the  details 
are  condensed  by  binning  the  distribution  on  a  32  by  32  grid  in  v^-v,  space 
for  plotting  purposes.  Binning  according  to  the  magnitude  of  vL  introduces  an 
additional  weighting  factor  of  vx,  which  is  divided  out  of  the  bin  occupation 
numbers  to  give  the  phase-space  densities  which  we  will  be  presenting. 
Velocity  moments  can  also  be  calculated  from  / . 

C.  High  Altitude  Ion  Conic 

The  first  application  of  the  Monte  Carlo  model  presented  here  will  be 
an  ion  conic  that  was  observed  at  high  altitude  by  the  S3-3  satellite  on  day 
280,  1976,  orbit  718.  This  conic  has  been  described  by  Temerin  et  al.  [2] 
and  by  Gomey  et  al.  [22].  The  conic  is  seen  in  conjunction  with  VLF  hiss, 
as  well  as  with  evidence  of  parallel  acceleration  of  electrons  at  energies  of 
1.6-5  keV  above  the  satellite,  which  was  at  an  altitude  of  5700  km.  The 
conic  is  folded  up  from  perpendicular  velocity  by  about  40°,  implying  that  the 
region  of  peak  transverse  heating  was  at  an  altitude  of  about  3600  km;  from 
the  width  of  the  conic  it  can  be  estimated  that  the  region  of  enhanced 
transverse  heating  was  about  700  km  in  extent  along  the  field  line.  This 
means  that  virtually  all  of  the  transverse  heating  of  the  conic  occurred  at  alti¬ 
tudes  below  the  spacecraft,  and  that  the  observations  at  the  satellite  of  the 
amplitude  of  the  auroral  hiss  will  be  of  little  help  in  estimating  the  rate  of  ion 
heating.  In  this  case  what  we  can  do  by  modeling  the  conic  is  to  infer  some 
of  the  characteristics  of  the  turbulence  which  would  be  required  to  explain  the 
remote  heating,  and  then  judge  how  well  these  characteristics  fit  our  expecta¬ 
tions  for  ;  sm. 
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To  model  the  heating  process,  then,  we  must  specify  a  model  for  the 
velocity  diffusion  rate  as  a  function  of  perpendicular  velocity  and  position 
along  the  field  line.  For  the  velocity  dependence  we  choose  the  rational  form 
given  above,  Eq.(35),  while  for  the  spatial  variation,  we  choose  a  power  law 
in  geocentric  distance.  We  do  not  make  any  claims  about  the  uniqueness  of 
this  form;  rather,  one  should  conclude  simply  that  it  is  a  form  that  can  pro¬ 
duce  conics  with  the  observed  shape.  Altogether,  the  diffusion  rate  is  written 
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where  SE0  is  a  reference  value  of  the  electric-field  spectral  density,  RE  is  the 
Earth’s  radius,  a  is  an  altitude,  a0  is  the  reference  altitude,  a  is  the  spatial 
power-law  index,  v,^  is  a  perpendicular  velocity  below  which  is  essen¬ 
tially  constant  with  a  small  value,  v0  is  the  reference  phase  velocity,  and  (3  is 
the  velocity  power-law  index  ((3=v-p  from  the  last  section). 

The  fitting  of  the  parameters  of  the  diffusion  rate  was  done  by  trial  and 
error,  by  comparing  the  phase-space  density  at  the  pitch  angle  of  maximum  / 
of  the  calculated  distribution  function  to  the  same  quantity  from  the  observed 
conic,  read  from  the  published  plots  of  Gomey  et  al.  [22].  The  resulting  fit  is 
shown  in  Fig.  10.  The  two  solid  curves  are  the  phase  space  densities  from 
the  observed  conic,  the  two  traces  being  from  the  contours  plotted  for  ±vA. 
For  this  calculation,  we  began  with  a  population  of  hydrogen  ions  traveling 
up  the  field  line  at  2800  km,  initially  with  a  Maxwellian  velocity  distribution 
of  temperature  1  eV.  The  plotted  distribution  function  was  tabulated  at  the 
altitude  of  the  satellite,  5700  km,  after  the  ions  crossed  a  region  containing 
auroral  hiss  turbulence  with  the  parameters;  5Ey=(20mV/m)2/kHz,  a0=  2S00 
km,  a=  9,  vmin=  35  km/s  (corresponding  to  an  energy  of  6.4  eV),  v0=  310 
km/s  (corresponding  to  500  eV),  and  (3=  3.  Note  that  this  value  of  Sg o 
represents  a  fairly  modest  value  of  the  total  electric  field  in  the  turbulence, 
whose  values  have  been  observed  to  range  up  to  50  mV/m  or  more.  The 
phase  velocity  of  the  peak  of  diffusion,  vQ,  corresponds  to  an  energy  of  500 
eV,  which  is  certainly  of  the  right  order  of  magnitude  if  the  turbulence  was 
generated  by  electron  beams  of  a  few  keV  energy  and  then  subjected  to 
degradation  of  phase  velocity  due  to  propagation  and  mode  coupling  effects. 
The  full  velocity-space  plot  of  the  calculated  ion  conic  distribution  is  shown 
in  Fig.  11.  Note  the  characteristic  folded  shape.  It  is  interesting  to  note  that 
although  the  model  for  the  diffusion  rate  decreases  monotonically  with  alti¬ 
tude,  the  conic  shows  its  peak  fiux  not  at  a  pitch  angle  corresponding  to  the 
minimum  altitude  of  the  simulation,  but  to  an  intermediate  altitude.  This  alti¬ 
tude  corresponds  to  the  point  where  the  ions  are  at  velocities  at  which  they 
experience  the  largest  transverse  heating,  that  is,  . 
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Figure  10.  The  phase-space  density  of  the  S3-3  conic  from  orbit 
718.  The  solid  curves  give  two  measurements  of  /  as  a  function 
of  energy  (in  eV)  at  the  pitch  angle  of  maximum  /,  while  the 
dot-dashed  curve  gives  the  fitted  theoretical  model. 

This  class  of  ion  conics  is  distinguished  in  form  from  the  second  class 
(second  only  in  order  of  discovery),  whose  formation  is  reviewed  by  Crew 
and  Chang  [i  1],  by  the  behavior  of  /  at  small  vA.  In  these  conics,  as  vA  goes 
to  zero,  the  parallel  velocities  of  the  ions  is  also  found  to  go  to  zero,  while  in 
the  other  conic  form,  the  ions  of  small  are  found  with  still  sizeable  values 
of  V|.  (Strictly  speaking  then,  the  other  type  of  transversely  accelerated  ion 
event  is  not  a  conic,  and  is  sometimes  called  an  ion  bowl,  for  the  shape  of  its 
velocity  distribution  in  vj-Vi  space.)  The  explanation  is  that,  due  to  the 
velocity-dependent  resonance  condition  that  low-velocity  ions  cannot  satisfy, 
the  bulk  of  the  ions  at  low  velocities  travel  up  the  field  lines  unaffected  by 
the  wave-particle  interaction  process  and  serve  to  anchor  the  velocity  distribu¬ 
tion  to  the  origin  in  velocity  space.  In  the  other  form  of  conic,  as  Crew  and 
Chang  [111  describe,  the  resonance  condition  is  independent  of  velocity  and 
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Figure  11.  The  calculated  ion  velocity  distribution  for  the  S3-3 
conic  of  orbit  718.  The  velocity  distribution  as  a  function  of  per¬ 
pendicular  (the  horizontal  axis)  and  parallel  velocities  (in  km/sec) 
is  presented  as  a  contour  plot.  The  bolder  contours  are  a  decade 
apart. 

all  ions  in  the  distribution  cam  participate  in  the  wave -particle  interaction, 
which  again  causes  velocity  diffusion  in  the  direction  perpendicular  to  the 
magnetic  field.  The  adiabatic  folding  effect  of  the  inhomogeneous  geomag¬ 
netic  field  acting  at  the  same  time  serves  as  a  one-way  street  for  the  flow  of 
energy,  from  perpendicular  degrees  of  freedom  to  the  parallel  direction. 
Thus,  while  the  parallel  energy  of  an  ion  is  always  increasing,  the  perpendic¬ 
ular  energy  is  being  continually  scattered  both  up  and  down,  and  at  any  time 
some  ions  of  a  given  parallel  energy  will  be  found  with  small  perpendicular 
energy. 
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D.  Low  Altitude  Ion  Conic 


Our  second  example  will  be  the  ion  conic  measured  during  the  sound¬ 
ing  rocket  flight  code-named  MARIE.  This  rocket  campaign  over  Greenland 
has  provided  us  with  what  is  probably  the  most  complete  and  detailed  set  of 
observations  within  the  regions  of  low-altitude  transversely  accelerated  ions 
yet  available  [3,45;  46].  The  rocket  passed  through  an  auroral  arc  on  its 
upleg  and  downleg  flights;  at  each  crossing  a  strong  correlation  was  seen 
between  the  intensity  of  the  auroral  hiss  electric  fields  at  frequencies  between 
3  and  16  kHz  and  the  flux  of  ions  flowing  up  the  field  line  at  energies  of  50 
to  370  eV,  with  pitch  angles  near  ninety  degrees.  The  intense  fluxes  of  both 
waves  and  particles  were  observed  at  altitudes  between  500  and  800  kilome¬ 
ters.  A  mass  spectrometer  on  board  determined  that  the  ions  were  primarily 
hydrogen,  consistent  with  our  expectation  that  the  lighter  ions  are  preferen¬ 
tially  heated  by  auroral  hiss  turbulence.  The  noted[46]  lack  of  direct  correla¬ 
tion  with  precipitating  electron  fluxes  probably  results  because  the  trajectories 
of  the  resonance -cone  waves  do  not  follow  the  geomagnetic  field  lines  like 
the  trajectories  of  the  precipitating  electrons  do.  There  is  some  suggestion  of 
parallel  heating  of  the  ambient  electron  distribution. 

To  model  this  event,  we  first  recognize  that  hydrogen  is  generally  a 
minority  species  of  the  plasma  at  these  altitudes,  where  the  major  constituent 
is  expected  to  be  oxygen.  Thus  the  hydrogen  ions  observed  here  have  prob¬ 
ably  been  drawn  down  from  higher  altitudes.  This  affects  the  residence  time 
of  the  ions  in  the  vicinity  of  the  intense  turbulence,  and  thus  the  energies  the 
ions  achieve.  Ior.s  traveling  down  the  field  line  enter  the  region  of  intense 
turbulence  and  begin  to  be  transversely  heated  by  interaction  wiui  the  waves. 
At  the  same  time,  the  ions’s  parallel  motion  slows  down  (and  eventually  rev¬ 
erses)  as  they  mirror  in  the  geometry  of  the  geomagnetic  field.  Thus,  the  ions 
spend  more  time  in  the  region  of  wave-particle  interaction  than  they  would  if 
they  were  already  traveling  up  the  field  line, increasing  the  effectiveness  of  the 
heating  process.  (Another  mechanism  to  note  that  enhances  residence  time  is 
effective  in  regions  where  the  steady  parallel  electric  field  is  directed  to 
accelerate  ions  downwards.  This  electric  field  causes  the  ions  to  bounce 
repeatedly  between  the  magnetic  mirror  below  and  the  electrostatic  mirror 
above,  until  the  ions  have  gained  enough  energy  to  pass  the  electrostatic  bar¬ 
rier  above.  This  “pressure-cooker”  effect,  studied  by  Gomey  et  a/.[47],  may 
be  effective  in  the  polar  cusp  region,  where  beams  of  electrons  accelerated 
out  of  the  ionosphere  are  observed.) 

For  our  simulation  of  the  MARIE  event,  then,  we  begin  with  a  hydro¬ 
gen  ion  species  with  a  temperature  of  1  eV  at  the  upper  end  of  the  simulated 
region,  from  500  to  800  km  altitude.  The  peak  spectral  density  of  the  electric 
field  observed  in  the  first  pass  through  the  arc,  from  about  05:29  to  5:31  UT, 
is  around  (7mV/m)2  per  kHz,  averaged  over  the  range  from  3  to  16  kHz. 
Fig.  12  sh  vs  the  results  of  the  simulation  as  two  contour  plots  of  the  ion 
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Figure  12.  Simulation  of  the  MARIE  ion  conic  event.  Contour 
plots  of  the  ion  velocity  distribution  as  a  function  of  perpendicu¬ 
lar  and  parallel  velocities  (in  km/sec)  at  700  km  altitude  (bottom) 
and  760  km  (top).  The  contours  are  a  half  decade  apart. 

phase-space  density  as  a  function  of  vL  and  at  the  altitudes  700  and  760 
km.  We  see  in  both  cases  the  characteristic  near-ninety  degree  pitch  angle  of 
the  conic,  despite  the  fact  that  the  headng  has  occurred  over  the  course  of 
several  hundred  kilometers.  Fig.  13  shows  the  phase-space  density  integrated 
over  pitch  angle  at  two  altitudes.  Yau  et  al.  [46]  reported  that  the  charac¬ 
teristic  energy  of  the  conic  was  about  22  eV,  with  a  tail  emerging  from  the 
body  of  the  distribution  at  around  150  to  200  eV,  extending  out  to  a  few  hun¬ 
dred  eV.  We  see  here  that  in  the  simulation,  the  observed  level  of  turbulence 
was  able  to  produce  the  heating  required  to  allow  the  ions  to  reach  these 
energies. 
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Figure  13.  The  calculated  phase-space  density  of  the  MARIE 
conic,  integrated  over  pitch  angle.  The  velocity  distribution  at 
the  altitudes  of  700  and  760  km  are  plotted  as  a  function  of  ener¬ 
gy  (in  eV). 

The  exact  shape  of  the  tail  distribution  could  have  been  fitted  by  adjust¬ 
ing  the  spadal  and  velocity  dependences  of  the  diffusion  coefficient,  which 
for  this  run  was  assumed  to  be  independent  of  both  variables.  From  our 
experience  with  the  S3-3  conic,  we  expect  that  the  phase  velocity  of  the  peak 
spectral  density  would  be  the  velocity  corresponding  to  an  energy  above  the 
point  where  the  tails  emerge  from  the  body  of  the  velocity  distribution, 
around  200  eV.  This  energy  corresponds  to  a  velocity  of  about  200  krrv/s. 
On  the  MARIE  rocket,  an  electric  field  interferometer  was  flown  to  attempt  to 
measure  wavelengths  of  the  observed  plasma  waves  [48].  When  the  rocket 
passed  through  the  ion-acceleration  region,  the  electric  field  interferometer 
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detected  a  feature  with  measureable  coherence  which  at  a  frequency  of  9  kHz 
had  a  wavelength  of  about  20  m.  This  represents  a  phase  velocity  of  180 
km/s,  in  agreement  with  the  phase  velocity  of  the  peak  spectral  density 
inferred  from  the  ion  velocity-distribution  measurement. 

To  summarize  the  points  on  which  the  MARIE  in-situ  plasma  and  wave 
measurements  agree  well  with  the  predictions  of  the  model  of  ion  interaction 
with  intense  auroral  hiss,  we  have  1)  a  strong  correlation  between  VLF  wave 
intensity  and  energetic  ion  flux;  2)  preferential  acceleration  of  hydrogen  over 
oxygen;  3)  enhancements  of  field-aligned  electron  fluxes;  4)  the  frequency 
spectrum  of  the  VLF  turbulence  around  the  lower  hybrid  frequency  shows 
structure  spaced  at  the  hydrogen  gyro  frequency,  suggesting  that  ion  damping 
of  the  waves  is  responsible  for  the  acceleration;  5)  the  amplitude  of  the  VLF 
turbulence  is  several  tens  of  mV/m,  strong  enough  to  account  for  the  ion 
energies;  and  finally,  6)  the  form  of  the  ion  velocity  distribution  has  the  high 
energy  tails  that  are  characteristic  of  the  acceleration  observed  in  the  simula¬ 
tion,  and  the  measured  phase  velocities  of  the  waves  match  the  features  of  the 
ion  distribution.  Many  of  these  points  have  been  observed  on  other  missions 
and  satellites,  but  never  perhaps  as  convincingly  as  observed  on  MARIE. 


V  DISCUSSION 

One  of  the  longlasting  controversies  in  the  study  of  ion  conics  near 
auroral  arcs  has  been  the  question  of  identifying  the  mechanism  for  the 
transverse  ion  heating,  whether  through  wave-panicle  interaction  with  the 
lower-hybrid  waves  reviewed  here  or  with  electrostatic  ion  cyclotron  waves 
[49]  or  through  other  mechanisms.  Searches  for  satellite  observations  of  con¬ 
ics  of  pitch  angle  near  ninety  degrees  (that  were  presumably  in  the  process  of 
being  accelerated)  were  inconclusive  [50],  but  suffered  because  fast  velocities 
and  slow  spin  rates  make  satellites  a  poor  platform  from  which  to  observe 
transverse  ion  acceleration  near  auroral  arcs.  The  advantages  of  sounding 
rockets  for  these  observations  were  demonstrated  in  the  MARIE  campaign(3), 
in  which  the  simultaneous  observation  of  transverse  ion  acceleration  and 
intense  lower-hybrid  waves  provided  strong  evidence  for  the  lower-hybrid 
heating  mechanism. 

To  model  the  observations  of  ion  conics  and  counter-streaming  elec¬ 
trons  near  auroral  arcs  requires  the  consideration  of  the  global  evolution  of 
both  the  particle  populations  and  the  wave  spectra  over  ranges  of  distance  and 
time  that  stretch  our  present  models.  The  complete  and  self-consistent 
analysis  of  the  fields  and  panic  es  near  auroral  arcs  remains  still  only  a  dis¬ 
tant  goal  of  space  physics.  The  origins  of  the  field-aligned  potential  drops 
that  accelerate  the  auroral  electrons  are  still  poorly  understood;  our  under¬ 
standing  of  the  spectra  of  the  plasma  waves  in  the  arc  environment  is  only 
beginning  to  develop.  For  the  calculation  of  ion  conics,  in  particular,  we  are 
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still  hampered  by  the  difficuhy  of  measuring  the  relevant  quantities  for  our 
theories,  and  thus  considerable  uncertainty  remains  in  the  empirical  estimation 
of  the  rate  of  transverse  ion  acceleration  because  of  the  lack  of  information 
on  the  wave-vector  spectrum  of  the  VLF  turbulence.  Efforts  to  measure  the 
wavelengths  of  the  plasma  waves  [48]  are  most  welcome,  although  such 
results  are  still  generally  limited  to  the  characteristic  or  dominant  wavelength 
rather  than  the  spectrum  of  wavelengths  that  would  be  needed  to  calculate 
wave-particle  interaction  rates.  Such  measurements  certainly  will  help  eluci¬ 
date  and  discriminate  among  the  variety  of  mechanisms  that  intervene 
between  the  linear  excitation  of  plasma  waves  in  the  auroral  arc  environment 
and  the  heating  and  acceleration  of  ionospheric  ions  and  electrons.  Using  the 
models  we  have  developed,  it  may  be  possible  to  use  the  form  of  the  ion 
conic  velocity  distribution  to  infer  constraints  on  the  spectral  shape  and  mag¬ 
nitude  of  the  VLF  turbulence,  a  great  benefit  when  the  acceleration  region  is 
localized  and  remote  from  the  platform  observing  the  accelerated  particles. 
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ABSTRACT 

Low  frequency  turbulence  present  on  closed  field  lines  in  the  central 
plasma  sheet  has  been  used  to  explain  ion  heating  and  conic  formation  with 
remarkable  success.  However,  the  source  of  the  turbulence  has  yet  to  be 
established,  and  there  are  no  obvious  local  sources  which  could  power  such  a 
broadband  spectrum.  Alternatively,  observations  reveal  that  ion  distributions 
in  the  equatorial  region  are  often  anisotropic,  and  such  distributions  excite 
waves  both  above  and  below  the  proton  gyrofrequency.  As  these  waves  pro¬ 
pagate  to  lower  altitudes  where  the  magnetic  field  is  stronger,  their  left-hand 
circularly  polarized  component  resonates  with  heavy  ions.  The  presence  of  a 
parallel  gradient  in  the  magnetic  field  complicates  the  details  of  wave  propa¬ 
gation,  and  as  a  result,  downcoming  right-hand  circularly  polarized  waves, 
which  acquire  a  left-hand  circularly  polarized  component  at  the  crossover  fre¬ 
quency,  may  tunnel  through  a  “stop  zone”  to  altitudes  where  they  resonate 
with  the  ions  and  thus  contribute  to  the  observed  ion  heating. 

I.  INTRODUCTION 

The  outflow  of  ionospheric  ions  in  the  form  of  ion  conics  into  the  cen¬ 
tral  plasma  sheet  (CPS)  is  often  accompanied  by  a  broadband  electric  field 
spectrum.  This  turbulence  which  exists  along  closed  field  lines  that  emanate 
from  the  auroral  region  can,  in  fact,  account  for  the  ion  heating  and  conic  for¬ 
mation  via  resonant  cyclotron  heating  in  the  presence  of  the  mirror  geometry 
of  the  Earth’s  magnetic  field  [1-4].  A  typical  electric  field  spectral  density  is 
shown  in  Fig.  I.  The  spectrum  is  peaked  at  low  frequencies  and  tends  to 
obey  a  power  law  over  the  range  of  frequencies  which  resonantly  match  the 
gyrofrequency  of  the  ionospheric  ions.  Such  spectra  are  commonly  observed 
up  to  several  RE  in  altitude.  Although  the  turbulence  and  conic  form  are  well 
correlated  for  the  events  studied  to  date  [5],  the  origin  and  nature  of  the  tur¬ 
bulence  remains  uncertain. 
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Figure  1.  Typical  electric  field  spectral  density  in  the  central 
plasma  sheet  (1=2.0  RE  ILAT  =  60°),  from  Ref.  (IJ. 

A  number  of  possible  sources  for  the  turbulence  exist  throughout  the 
magnetospheric  environment.  Of  course,  the  most  obvious  possibility  for  a 
source  would  be  one  that  is  local  in  nature;  however,  all  of  the  events  scrutin¬ 
ized  to  date  fail  to  point  out  any  local  sources  and  instead  suggest  that  a  more 
global  approach  is  required  [1,  6],  Velocity  shears,  for  example,  observed  in 
the  cusp  region  might  excite  the  Kelvin-Helmholtz  instability  and  generate 
Alfvdn  waves  which  might  propagate  upwards  into  the  heating  region. 
Alternatively,  instabilities  in  the  magnetotail  associated  with  reconnection 
might  generate  waves  which  propagate  earthward  into  the  central  plasma 
sheet.  In  this  paper,  however,  we  will  consider,  in  particular,  the  possibility 
of  a  nonlocal  source  emanating  from  the  equatorial  region. 

n.  EQUATORIALLY  GENERATED  WAVES 

Anisotropic  ion  distributions  are  often  observed  in  the  nightside  equa¬ 
torial  region  at  geosynchronous  orbit  [7].  As  ions  are  injected  earthward 
from  the  magnetotail,  they  drift  adiabaticaliy  about  the  Earth  along  drift 
shells.  Because  the  shell  on  which  a  particle  is  constrained  to  move  depends 
on  the  equatorial  pitch  angle,  initially  isotropic  distribution  functions  tend  to 
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develop  thermal  anisotropies  (T±>  T „)  and  acquire  anomalous  loss  cones 
(enlarged  over  the  3°  atmosphenc  loss  cone)  [8,  9],  In  addition,  the  observa¬ 
tions  reveal  that  these  distributions  often  contain  a  large  energetic  population 
with  90°  pitch  angle  [7|,  Moreover,  the  proton  gyrofrequency  near  geosyn¬ 
chronous  orbit  in  the  equatorial  region  (/'£■)  matches  the  oxygen  cyclotron 
frequency  (fc0 .)  over  the  range  of  altitudes  where  heating  takes  place.  Below 
we  consider  how  anisotropic  ion  distributions  can  generate  waves  with  fre¬ 
quencies  near  //J-  in  the  equatorial  region  and  how  these  waves  may  reach 
the  CPS  at  altitudes  ranging  from  the  upper  ionosphere  to  several  RE. 

A.  Waves  with  /  <  //*- 

Ion  distributions  with  loss-cones  or  temperature  anisotropies  can  excite 
waves  below  the  proton  gyrofrequency.  Indeed,  this  instability  has  been  stu¬ 
died  comprehensively  in  the  equatorial  region  in  connection  with  Pci  (0.2  < 
/  <  5  Hz)  wave  generation  and  helium  heating  [10-12].  Waves  generated 
with  frequencies  below  the  helium  gyrofrequency  (typically  /  <  1  Hz)  are 
guided  along  field  lines  to  the  ground  [13,  14],  On  the  other  hand,  it  is 
believed  that  waves  generated  above  the  helium  gyrofrequency  but  below  the 
hydrogen  cyclotron  frequency  ( 1  <  /  <5  Hz),  rather  than  penetrating  directly 
to  the  ground,  travel  between  geomagnetic  conjugate  points  near  the  magnetic 
equator  which  correspond  to  the  ion-ion  hybrid  frequency.  Although  the 
growth  rates  for  these  waves  are  smaller  than  those  for  the  waves  near  the 
helium  gyrofrequency,  the  power  in  the  waves  may  be  comparable  because 
they  amplify  as  they  bounce  back  and  forth  through  the  equatorial  plane  [13]. 
Furthermore,  when  the  concentration  of  helium  is  low  enough,  these  waves 
tunnel  through  the  reflection  layer  to  lower  altitudes  in  the  right  hand  circu¬ 
larly  polarized  (RHCP)  mode  (see  Fig.  6  of  Ref.  [14]).  In  fact,  ground  based 
instruments  often  detect  both  types  of  waves — an  observation  which  strongly 
suggests  that  they  contribute  to  the  auroral  turbulence  [14], 

B.  Waves  with  /  > 

Propagating  waves  generated  below  may  explain  a  significant  part 
of  the  wave  power  observed  at  auroral  altitudes,  but  much  of  the  wave  power 
has  /  >  5  Hz  (Fig.  1).  Moreover,  the  initial  energization  of  ionospheric  oxy¬ 
gen  to  magnetospheric  energies  occurs  at  low  altitudes  where  /<0.  is  larger 
than  Pci  frequencies.  Therefore,  it  is  important  to  consider  the  generation 
and  propagation  of  waves  with  /  >  f  ^  in  the  equatorial  region. 

The  measurements  taken  from  GEOS  1  and  2  indicate  that  ion  distribu¬ 
tions  with  loss-cone  character  ( dF/dvL  >  0)  exist  at  geosynchronous  orbit  in 
the  equatorial  region  and  can  generate  waves  at  multiples  of  f’} J.  [7].  Indeed 
Perraut  v  al.  17]  modeled  the  hot  ion  population  with  a  ring  distribution  and 
exami"i’'_  the  solutions  of  the  linear  dispersion  relation  with  kt  =  0.  They 
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found  that  indeed  waves  are  generated  at  multiples  of  f’#,  but  the  waves 
have  a  group  velocity  perpendicular  to  the  magnetic  field  which  implies  that 
they  are  confined  to  the  equatorial  region. 

On  the  other  hand,  we  have  solved  the  dispersion  relation  of  linear 
waves  using  the  somewhat  more  realistic  model  used  by  the  WHAMP  code 
[15,  16].  In  particular,  we  have  included  a  thermal  spread  in  the  hot  ion  dis¬ 
tribution  function  which  we  have  taken  to  be  of  the  form 
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with  the  parameters  listed  in  Table  I.  The  distribution  consists  of  a  cold  back¬ 
ground  plasma  with  a  small  energetic  proton  population  which  peaks  in  the 
perpendicular  direction  at  7  keV. 


Table  I. 

Parameters  Used  in  Eq.  (1) 

Species 

T  (keV) 

n  (cm-3) 

<*i 

02 

A 

H+ 

0.001 

65.0 

1.0 

0.0 

H+ 

1.0 

1.6 

1.0 

0.0 

e~ 

0.00 1 

66.6 

1.0 

0.0 

We  find  that  for  k h  =  0  electrostatic  waves  are  generated  at  multiples  of 
//d.  with  group  velocity  perpendicular  to  the  magnetic  field;  however,  for 
k >  10^*  a  broad  band  of  unstable  electromagnetic  modes  arises.  In  Fig.  2 
we  plot  the  dispersion  surface  of  the  real  frequency  in  a  two  dimensional 
wave  vector  space,  and  in  addition,  we  indicate  the  growth  rates  for  the  par¬ 
ticular  mode  with  shading.  This  figure  illustrates  that  for  larger  kt  the  har¬ 
monic  spectrum  at  *t  =  0  becomes  more  broadband.  In  addition,  these  broad¬ 
band  mtxies  with  kn  *  0  also  have  vgllvgL~  1  as  is  illustrated  in  Fig.  3.  In 
this  figure,  we  have  taken  a  slice  at  fixed  fcj  along  one  band  in  Fig.  2.  We 
have  plotted  one  curve  which  corresponds  to  maximum  growth  at  =  0  and 
another  curve  which  is  slightly  offset  in  kL  from  the  first  curve.  This  figure 
illustrates  that  although  the  maximum  growth  may  be  in  the  modes  which  are 
electrostatic  and  propagate  perpendicular  to  the  magnetic  field,  electromag¬ 
netic  modes,  which  excite  a  much  larger  region  of  k  -space,  have  comparable 
growth  rates  and  also  have  1.  Thus,  the  electromagnetic  waves  will 

initially  propagate  out  of  the  equatorial  region.  Increasing  the  temperature 
and  density  of  the  energetic  population  increases  the  growth  rates  as  one 
would  expect;  however,  for  sufficiently  warm  background  electrons  or  ions. 
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Figure  2.  Dispersion  relation  (co  vs.  logio logio(lfcip,)  where 
p,  is  the  ion  gyroradius)  using  the  distribution  in  Eq.  (1)  with  the 
parameters  listed  in  Table  I.  Shading  indicates  the  growth  rates 
(white  for  growing,  dark  for  damping)  for  waves  on  the  surface. 

Notice  that  electrostatic  waves  are  generated  at  multiples  of  /fir 
which  are  largest  at  k\  =0  whereas  the  waves  at  larger  kt  tend  to 
be  electromagnetic  and  occupy  a  much  larger  region  of  k- space. 

the  in  *  0  modes  become  damped. 

If  we  compare  the  power  in  the  equatorial  and  auroral  waves  using  the 
observed  wave  amplitudes,  we  may  determine  whether  the  free  energy  in  the 
equatorial  region  is  a  viable  source  for  the  turbulence.  With  the  assumptions 
that  the  auroral  waves  are  essentially  Alfv^n  waves  and  that  the  WHAMP 
results  accurately  reflect  the  relationship  between  the  observed  magnetic  field 
and  the  electric  field  in  the  equatorial  region,  we  find  that  the  Poynting  fluxes 
are  comparable.  Although  detailed  ray  tracing  calculations  will  be  required  to 
confirm  the  notion  that  the  auroral  turbulence  originates  in  the  equatorial 
plane,  this  simple  comparison  elicits  some  confidence. 


—  log10(k1pi)=-0.804  log  io(kjDi)  =  - 0.800 


Figure  3.  Growth  rates  (a)  and  group  velocity  ratio  /v?i  (b) 
for  a  plasma  with  parameters  listed  in  Table  I.  The  curves  are 
taken  at  constant  fcj),  where  p,  is  the  ion  gyroradius.  The 
dashed  curve  corresponds  to  maximum  growth  for  =0  modes 
near  the  fourth  harmonic  of  the  proton  gyrofrequency  with 
logiof/tj), )  =  -0.804.  The  solid  curve  corresponds  to  modes  with 
logiofijp, )  »  -0.800  which  are  slightly  offset  from  maximum 
growth  at  k%  =0.  It  is  to  be  noted  that  these  modes  have  a  slight¬ 
ly  smaller  growth  rate,  but  they  also  have  v  |(  ivg[  >  1. 

01.  WAVE  PROPAGATION 

In  the  previous  section  we  identified  an  instability  in  the  equatorial 
region  that  generates  electromagnetic  waves  which  initially  propagate  earth¬ 
ward.  Such  an  instability  would  be  a  likely  source  for  the  auroral  turbulence 
observed  at  lower  altitudes.  In  order  to  determine  whether  the  auroral  region 
is  accessible  to  these  waves,  we  now  consider  the  characteristics  of  elec¬ 
tromagnetic  wave  propagation.  In  a  cold  plasma  only  two  electromagnetic 
modes  exist  at  a  given  frequency.  Because  only  the  left  hand  circularly 
polarized  (LHCP)  component  of  these  waves  resonates  with  ions,  the  waves 
must  contain  a  non-negligibie  fraction  of  LHCP  at  the  local  oxygen  cyclotron 
frequency  in  order  for  heating  to  occur.  In  a  cold  plasma,  however,  the 


88. 


branch  which  is  RHCP  for  kv  =  0  contains  no  LHCP  at  the  ion  gyrofrequency 
even  for  *  0.  The  waves  that  we  are  considering  are  generated  with 
/  >  //it  and  thus,  are  most  likely  RHCP  because  the  LHCP  branch  is  cut  off 
above  /cK-.  In  order  to  heat  ions,  these  waves  must  acquire  some  fraction  of 
LHCP  as  they  propagate  to  lower  altitudes.  Such  a  transmutadon  of  wave 
polarization  may  occur  via  mode  conversion  and  tunneling. 

In  order  to  understand  how  RHCP  waves  might  gain  LHCP  as  they 
propagate  to  auroral  altitudes,  we  consider  linear  theory  for  an  idealized 
model  consisting  of  a  two  ion  component  plasma  as  depicted  by  the  disper¬ 
sion  relation  in  Fig.  4.  Because  wavelengths  are  very  long  and  the  density 
does  not  vary  appreciably  [17],  we  assume  that  the  plasma  is  cold  with  no 
variation  in  density.  Furthermore,  we  assume  that  the  magnetic  field  varies 
slowly  along  itself  and  that  the  wavelength  X<c LB  =  [d\n(B)/dz]~l.  This 


Figure  4.  Dispersion  relation  for  a  cold  two  ion  component  plas¬ 
ma  for  nonparallel  propagation.  The  various  branches  are  labeled 
according  to  their  predominate  polarization  with  RHCP  dashed 
and  LHCP  solid.  Coupling  occurs  at  the  crossover  frequency 
(cocr)  where  the  polarization  of  branches  I  and  II  switches  abrupt¬ 
ly.  Tunneling  occurs  across  the  shaded  "gap”  between  the 
cutoff  frequency  (0)C(?)  and  the  heavy  ion  gyrofrequency  (flcl) 
from  branch  I  to  ID. 
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assumption  reduces  the  problem  to  one  dimension  and  fixes 

In  the  ensuing  analysis  we  consider  a  particular  scenario  whereby 
waves  generated  with  /  >  fc\ propagate  earthward.  The  wave  initially  pro¬ 
pagates  downward  in  the  RHCP  mode  (I  in  Fig.  4).  At  the  crossover  fre¬ 
quency  (0)cr),  where  linear  theory  predicts  that  the  predominate  polarization 
of  modes  I  and  II  reverses,  mode  coupling  may  occur  [18].  As  the  modes 
couple,  some  fraction  of  the  wave  power  will  be  transferred  to  mode  II  and 
will  remain  predominately  RHCP.  For  an  idealized  cold  plasma,  mode  II 
contains  no  LHCP  at  the  heavy  ion  gyrofrequency  and  thus,  this  transferred 
wave  power  will  not  heat  ions.  Nevertheless,  the  rest  of  the  wave  power 
remains  in  mode  I,  which  becomes  predominately  LHCP  for  co  <  cocr .  As  the 
wave  propagates  along  mode  I  to  lower  altitudes,  it  encounters  the  cutoff  fre¬ 
quency  (coco)  where  the  wave  is  reflected.  However,  some  of  the  wave 
power  may  tunnel  through  the  evanescent  layer  to  mode  III  where  the  waves 
can  resonate  with  the  ions.  Below,  we  calculate  the  coupling  (C),  transmis¬ 
sion  (T),  and  reflection  (/?)  coefficients  where  C  is  the  fraction  of  the  wave 
power  that  remains  in  mode  I  as  the  wave  passes  through  0)cr,  T  is  the  frac¬ 
tion  of  wave  power  that  tunnels  from  mode  I  to  mode  III  at  (Oco ,  and  R  is  the 
fraction  of  wave  power  that  is  reflected  back  along  mode  I  at  coco . 

To  proceed  with  the  analysis.  Maxwell’s  equations  combined  with  the 
continuity  and  momentum  transfer  equations  yield  coupled  differential  equa¬ 
tions  for  the  electric  field.  Because  the  waves  that  we  are  considering  have 
£ ii  =  0  [19]  the  equations  for  the  electric  field  components  take  the  simple 
form 


where 

£t  =  Ex  ±  iEy  e  =  k?c2/2<s>2  (3) 

and 


_  y 


(4) 


where  o)pj  and  Q,  are  the  plasma  and  cyclotron  frequencies  for  each  species 

s. 
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A.  Coupling  at  the  Crossover  Frequency 

Because  the  functions  R  and  L  vary  slowly  in  altitude,  we  may  apply 
the  WKB  analysis.  Therefore,  given  the  boundary  conditions  (downcoming 
RHCP  wave),  the  zeroth  order  (in  e  «*:  1)  electric  fields  take  the  form 

l 

E^=E0{R(z))-Viexp(-ij[R{s)}u2ds)  ;  £'0)=0  (5) 

where  h  is  some  suitably  chosen  reference  altitude.  The  parameter  e  couples 
the  two  parallel  propagating  modes  so  that  we  may  solve  for  the  electric  field 
perturbatively.  To  first  order  in  e  we  obtain  the  equation  for  the  coupled  elec¬ 
tric  field  £+. 


d2£}1) 


- f-  +£(z)£j1)  =  — £  £(0) 

dz1 

(6) 

We  may  solve  this  equation  via  variation  of  parameters  to  obtain 

+  A_<j>_ 

(7) 

where 

Z 

At(z)  =  ~jds  [L(s)R(s))-Vi  x 

11  H 

(8) 

z 

?rv(-i  j([R(<!)}'r2±[L{q)}'r2)dq  ) 

h 

and 

Z 

<|>±(z )  =  £° [L (z  )]"l/4  exp(  ±ij[L  {s )] [ads  ) 

h 

(9) 

The  amplitude  of  the  downgoing  coupled  wave  A_  involves  a 
over  an  oscillating  function  which  becomes  stationary  at  the 
quency 

phase  integral 
crossover  fre- 

n  1 

(10) 

defined  by  D=(R  -L)/ 2  =  0.  Because  the  contribution  to  the  overall 
integral  in  (8)  is  exponentially  small  except  at  the  stationary  point,  we  may 
expand  the  integral  in  the  integrand  of  (8)  about  the  stationary  point.  We 
may  then  extend  the  limits  of  the  overall  integration  to  be  ±<*»  which  intro¬ 
duces  only  exponentially  small  errors.  Then  we  find  that  the  coupling 
coefficirm  is,  up  to  an  unimportant  phase,  given  by 
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A  _(-°°) 


1 


2  [R  (coc,)]l/2ji. 


jds  exp^ 


-i  D'( cocr)  s2  \ 
2  [/?  (cocr  )]1/2  ) 


(It) 


where 


D\  to  =oicr) 


2c  2  ^tQpl^el't‘(dp2/Qc2)  (nl~*~n2)2 

(Qc2!-Qc22)2  rtln2 


(12) 


We  may  then  rotate  the  path  of  integration  to  evaluate  the  trivial  Gaussian 
integral  in  (11).  By  comparing  the  Poynting’s  flux  of  the  incident  wave 
above  the  coupling  region  (mode  I,  RHCP)  to  that  of  the  coupled  wave 
below  the  coupling  region  (mode  I,  LHCP)  we  may  obtain  the  coupling 
coefficient 


4 below 

TC  e2 

if 

J  above 

2  'fRD' 

(13) 


Thus,  although  the  dispersion  curves  appear  to  be  continuous  and  separate  for 
modes  I  and  II,  the  presence  of  a  gradient  does  introduce  coupling.  At  small 
angles  of  propagation  most  of  the  wave  power  is  transferred  from  mode  I  to 
II  in  which  case  the  waves  remain  RHCP,  while  at  larger  angles  of  propaga¬ 
tion  much  of  the  wave  power  remains  in  mode  I,  which  is  predominately 
LHCP  below  cocr . 


B.  Tunneling  Through  the  Evanescent  Layer 


Once  the  downcoming  RHCP  wave  has  gained  some  LHCP  at  the 
crossover  frequency,  the  wave  encounters  the  cutoff-resonance  pair  of  fre¬ 
quencies.  At  this  point,  intermode  coupling  is  small  so  that  we  neglect  e  in 
the  ensuing  calculations.  We  are  now  concerned  with  the  lower  component 
of  equations  (l)  which  may  be  written 


dZE+  _  CD2  flc^c2  M  ~(0co  ^ 

dz 2  ”  2  wco  (o)-Ocl)(a>-Oc2)  * 


(14) 


where 


®co 


n1Qc2+/l  2^cl 


n{+n: 


=  s 

species 


4  nn,  ms 
B2 


(15) 


If  we  linearize  about  Oci  we  find  that  this  equation  is  of  the  Budden  fonr 
[20].  The  transmission  and  reflection  coefficients  are  then  given  by 

T  =  exp(-n(})  /?  =(1— exp(— reP))2  (16) 


where 
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(17) 


Qc  i  ( 1  -m2lm  i)1'’2 

(3  -  — —Lg _ - _ - _ 

H  vA  B  (1  +n2m2ln  xm  j) 

As  is  characteristic  of  Budden  Tunneling  R+T  *  1,  and  the  missing  wave  flux 
is  absorbed  by  the  medium.  The  absorption  typically  occurs  because  wave 
power  is  converted  to  a  thermal  branch  which  heats  the  ions  [20], 

C.  Numerical  Results 

We  consider  an  0+-H+  plasma  in  the  auroral  region  at  2  RE  along  a 
field  line  with  invariant  latitude  60°.  If  we  take  parameters  typical  of  this 
region:  n  -  100cm-3,  B  =  0.05  G,  /cH-  =  80  Hz,  J  0.  =  5  Hz,  LB  =  1  RE,  we 
find  that  for  sufficiently  low  concentrations  of  O*  a  few  percent  LHCP  is 
obtained  at  the  heavy  ion  gyrofrequency  via  mode  conversion  and  tunneling 
as  indicated  in  Table  II.  These  results  are  consistent  with  ion  heating  theory, 
which  required  only  a  small  fraction  of  the  waves  to  be  LHCP  [2]. 


Table  II.  Transmission  and  Coupling  Coefficients 

nO’/nH-  tan-1(*i^ii) 

C 

T 

R 

0.05  18° 

10% 

0 

100% 

0.01  27° 

10% 

18% 

67% 

IV.  DISCUSSION 

In  this  paper  we  have  identified  a  plausible  source  region  for  the  tur¬ 
bulence  associated  with  ion  conics.  It  is  evident  that  equatorially  generated 
low  frequency  (/  <  f‘£*)  waves  do  propagate  to  low  altitudes  where  they  may 
contain  sufficient  LHCP  to  heat  the  ions.  In  addition,  we  have  determined 
that  waves  at  higher  frequencies  may  also  be  generated  in  the  equatorial 
region,  where  they  initially  propagate  earthward.  Careful  ray  tracing  will 
demonstrate  whether  these  waves  do  indeed  propagate  to  auroral  altitudes. 

The  high  frequency  waves  (/  >  f^r)  are  important  because  they  may 
be  responsible  for  extracting  oxygen  from  the  ionosphere  where  /t0.  >  f,] J.. 
Although  the  wave  power  at  high  frequencies  provides  a  resonant  mechanism 
whereby  a  small  fraction  of  the  ionospheric  ions  may  be  accelerated  to  mag- 
netospheric  energies,  at  low  altitudes  where  the  oxygen  cyclotron  frequency  is 
increased,  the  power  available  to  heat  the  ions  is  diminished  by  several  orders 
of  magnitude.  On  the  other  hand,  the  spectral  density  seems  to  be  peaked  at 
the  lowest  frequencies.  This  fact  elicits  the  question  of  whether  at  some  alti¬ 
tude,  non-resonant  heating,  resulting  from  small  fluctuations  in  the  fields, 
might  be*,  me  important  [21,  22]. 
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We  have  also  demonstrated  that  the  effects  of  mode  conversion  and 
tunneling  become  significant  in  the  presence  of  parallel  magnetic  field  gra¬ 
dients.  In  the  special  case  that  we  considered,  we  found  that  at  the  crossover 
frequency,  power  in  the  RHCP  branch  may  be  transferred  to  the  LHCP 
branch,  and  at  lower  altitudes,  the  waves  tunnel  through  an  evanescent  layer 
and  resonate  with  the  ions.  Of  course,  a  similar  analysis  would  apply  to 
waves  generated  in  the  LHCP  mode.  Furthermore,  the  addition  of  another 
ion  species  would  complicate  wave  propagation  further  and  introduce  another 
crossover  frequency  and  resonance. 

Consideration  should  also  be  given  to  other  effects  which  might  effect 
wave  propagation  and  heating.  Although  the  inclusion  of  finite  temperature 
should  not  affect  the  crossover  coupling,  it  will  affect  the  resonance-cutoff 
tunneling  by  introducing  higher  ordei  derivatives  which  become  important  at 
the  singularity  of  the  differential  equation  (7).  However,  careful  considera¬ 
tion  of  the  effects  of  finite  temperature  typically  demonstrates  that  the  missing 
wave  flux  goes  into  a  thermal  branch  which  ends  up  heating  the  plasma.  In 
addition,  nonlinear  coupling  of  the  turbulence  may  allow  more  heating  to  take 
place  near  the  gyrofrequency. 

In  conclusion,  the  notion  that  the  CPS  turbulence  arises  in  the  equa¬ 
torial  region  appears  to  be  a  very  promising  possibility,  and  the  details  of 
wave  propagation  from  this  region  to  auroral  altitudes  will  be  the  subject  of 
further  investigation. 
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ABSTRACT 

The  one  dimensional,  nonlinear  ion  hole  model  of  double  layer  forma¬ 
tion  on  auroral  field  lines  is  modified  to  include  three  dimensional  effects.  In 
this  expanded  model,  the  hole  potential  drives  electrostatic  ion  cyclotron 
waves  in  the  background  plasma,  much  as  in  a  discrete  particle  model.  The 
field  aligned  dynamics  (trapping)  which  result  in  hole/double  layer  formation 
are  retained  as  in  the  one  dimensional  model.  The  three  dimensional 
hole/double  layer  structure  has  parallel  and  perpendicular  scale  lengths  in 
agreement  with  Viking  satellite  data. 

I.  INTRODUCTION 

Nonlinear  ion  hole/clump  instability  [1,  2]  has  been  proposed  as  the 
cause  of  double  layers  observed  in  the  auroral  particle  acceleration  region  [3]. 
An  ion  hole  is  a  localized  depletion  in  the  ion  phase  space  density.  Because 
of  the  neutralizing  charge  density  of  the  plasma  electrons,  ions  become 
trapped  in  the  hole’s  negative  potential.  Such  a  self-trapped,  virialized  phase 
space  structure  is  a  BGK  mode. 

A  novel  feature  of  the  instability  is  that  the  holes  can  form  and  grow  in 
linearly  stable  plasma  where  (for  example),  for  weak  currents  and  Te  =  Ti ,  ion 
acoustic  waves  can  not  grow.  An  ion  hole  grows  as  electrons  arc  resonantly 
reflected  by  the  ion  hole's  negative  potential.  In  a  drifting  plasma,  the  elec¬ 
tron  reflection ‘gives  the  ion  hole  a  double  layer  potential  structure,  and  causes 
the  hole  to  accelerate  through  the  ion  distribution  function.  As  a  result,  the 
plasma,  in  its  fully  developed  turbulent  state,  is  composed  of  an  intermittent 
distribution  of  colliding,  growing  hole  fluctuations  (clumps).  Though  the 
hole/clump  model  of  Refs.  [1-3]  is  one  dimensional,  its  predictions  agree 
well  with  double  layer  observations  from  the  S3-3  and  Viking  satellites.  This 
agreement  results  from  the  fact  that,  in  a  strongly  magnetized  plasma,  the 
magnetic  field  restricts  the  particle  trapping  to  one  dimension  (the  component 
of  velocity  along  the  magnetic  field). 

However,  the  satellite  data  indicates  that  the  double  layers  have  three 
dimensional  structure,  the  perpendicular  scale  of  the  double  layer  being  some¬ 
what  larger  than  the  field  aligned  dimension  14-71  In  addition,  the  data 
reveals  ihe  presense  of  electrostatic  ion  cyclotron  (EIC)  waves.  We  show  here 
how  the  model  of  Ref.  [3]  can  be  extended  to  three  dimensions.  A  fully  3D 
model  cf  touble  layer  formation  in  a  magnetic  field  requires  a  fully  nonlinear 
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treatment  of  hole  dynamics  in  three  dimensions.  As  in  the  nonlinear  hole 
dynamics  in  ID,  a  fully  three  dimensional  hole  model  requires  an  extensive 
investigation  of  hole  dynamics  from  first  principles.  The  resulting  model 
becomes  significantly  more  complex  than  in  the  ID  case.  For  example,  in 
three  dimensional  geometry  with  a  magnetic  field,  it  is  possible  for  groups 
(“clumps”)  of  ions  to  become  correlated  along  their  gyro  orbits.  Each  such 
clump  of  ions  (“gyro-clump”)  will  gyrate  ballistically  around  a  magnetic 
field  line.  Conceptually,  we  would  expect  such  nonlinearly  enhanced  particle 
correlations  to  effect  the  perpendicular  scale  of  the  holes.  In  addition,  ion 
cyclotron  wave  emmision  would  be  nonlinearly  enhanced.  The  mathematical 
evaluation  of  these  collective,  nonlinear  effects  from  the  gyro-kinetic  Vlasov 
equation  is  formidable. 

Here,  we  present  a  more  modest,  but  certainly  not  trivial,  extension  of 
the  ID  model.  In  this,  the  ID  dynamics  of  ion  trapping  along  the  field  lines  is 
fully  nonlinear.  However,  the  perpendicular  dynamics  is  treated  linearly.  We 
show  that  the  essential  hole  dynamics  of  Ref.  [3]  carry  over  to  the  3D  case, 
but  with  the  hole  linear  dielectric  function  generalized  to  the  appropriate  3D 
analogue.  The  model  allows  for  the  existence  of  both  parallel  ion  (velocity 
space)  trapping  along  the  field  lines  (ion  hole  formation  as  in  Refs.  [1-3])  and 
ion  cyclotron  motion  perpendicular  to  the  field  lines  (3D  double  layer  struc¬ 
ture  as  well  as  the  generation  of  EIC  waves).  In  this  way,  we  show  that  the 
double  layers  and  EIC  waves  can  be  integrated  naturally  into  the  same  non¬ 
linear  model.  The  model  demonstates  that  hole/clump  fluctuations  are  non¬ 
linearly  unstable  to  growth  in  a  3D  plasma  with  strong  magnetic  field. 
Specifically,  we  find  that  the  double  layers  can  be  the  source  for  EIC  waves, 
much  as  in  the  Cerenkov  emission  of  waves  by  discrete  panicles.  Of  course, 
as  we  have  just  noted  above,  the  nonlinear  formation  of  gyro-clumps  will 
enhance  this  Cerenkov  emission.  However,  preliminary  investigation  shows 
that  these  enhancements  appear  to  make  quantatative  (of  order  unity)  rather 
than  qualitative  changes  in  the  model  we  present  here. 

II.  THREE  DIMENSIONAL  MODEL 


We  invoke  the  disparity  between  the  ion  dynamics  parallel  and  perpen¬ 
dicular  to  the  magnetic  field.  The  characteristic  time  for  the  parallel  ion 
motion  is  the  trapping  time  (Axj/Av,)  of  an  ion  that  is  trapped  in  the  negative 
potential  ($)  of  an  ion  phase  space  hole.  The  ion  gyro  frequency  (Q,) 
characterizes  the  ion  perpendicular  dynamics.  We  note  that 
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where  v,  and  p(  are  the  ion  thermal  velocity  and  gyro  radius,  and  we  have 
assumed  that  Tt  =  Ti .  With  this  disparity  in  time  scales,  the  Vlasov  equation 
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splits  up  into  a  response,  5/5,  for  parallel  ion  motion,  and  a  response,  8/[, 
for  perpendicular  ion  motion.  The  ion  hole  can  thus  be  thought  of  forming  in 
a  background  plasma  of  gyrating  ions. 

Assuming  that  the  perpendicular  ion  dynamics  are  linear,  the  Fourier 
transform  of  8/[  is  given  by  the  usual  linear  response, 

5/{(k,oi)  =  y — -7—77/0, <Kk’Co)  .  (2) 

cu2-n,2 

Poission’s  equation  is 

-JL-V20  =  Jtfvny/o.^  -  jdvt8ft  -  jdvL8f[  ,  (3) 

where  the  first  term  on  the  RHS  is  the  electron  response.  With  (2),  the 
Fourier  transform  of  (3)  can  be  written  as 


0(k,<o) 


0(k,(o) 
e(k,to)  ’ 


(4) 


where 


(5) 

(6) 


is  the  linear  dielectric  function  for  EIC  waves.  Accordingly,  (4)  describes  the 
shielding  of  the  hole/double  layer  potential  by  the  background  plasma.  The 
double  layer  potential  can  excite  EIC  waves  in  the  background  plasma,  much 
as  in  a  discrete  particle  model.  The  parallel  (trapping)  dynamics  of  the  dou¬ 
ble  layer  is  as  described  in  Ref.  [3],  and  enters  the  3D  model  here  via  (5)  in 
the  numerator  of  (4).  The  perpendicular  dynamics  enters  through  the  dielectric 
function  in  the  denominator  of  (4).  While  the  existence  of  EIC  waves  in  the 
auroral  acceleration  region  may  also  be  do  to  other  dynamical  mechanisims 
than  the  Cerenkov  emission  model  (4),  the  3D  model  presented  here  shows 
that  die  hole/double  layers  can  form  and  grow  in  a  background  plasma  con¬ 
taining  EIC  waves. 

The  hole/double  layers  are  described  by  the  low  frequency 
(d)  =  Jfc|Avi  <  Q,)  portion  of  (4).  For  this  frequency  range,  we  revert  back  to 
perpendicular  position  and,  using  (5)  and  (6),  write  (4)  as 


99. 


0U.,ki,a>) 


(7) 


=  inejdvnSftixi.k^a)) 


From  (7),  the  ratio  of  perpendicular  to  parallel  scale  length  for  a  3D  double 
layer  fluctuation  is 


Ax,,  2  l+(£ii ^-o)~2 

.^7J  '  1+(P,^d)2 


(8) 


This  ratio  is  on  the  order  of  unity  or  less  for  typical  parameters  and  is  con¬ 
sistent  with  the  satellite  measurements  [4-6]. 
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